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Abstract 

In this paper we show that there exist analytic initial data in the stable regime for the Muskat 
problem such that the solution turns to the unstable regime and later breaks down, i.e., no longer 
belongs to C 4 . 
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1 Introduction 

The Muskat equation governs the motion of an interface separating two fluids in a porous medium 
(e.g. oil and water in sand). See [ 10 1 . The same equation governs an interface separating two 
fluids trapped between two closely spaced parallel vertical plates (a "Helle Shaw cell"). See ]9). 

In this paper, we exhibit a solution of the Muskat equation for which the interface is initially 
smooth and stable, but later enters an unstable regime, and still later develops a singularity. 

Let us briefly recall the Muskat equation and its derivation from physical laws. Imagine the 
plane M 2 filled with two incompressible fluids that cannot mix. At time t, let be the region 
occupied by FLUID 1, and let fi 2 (t) be the region occupied by FLUID 2. These two regions are 
separated by their common boundary c?!! 1 (t) = dfl 2 (t) (the "interface"), which we represent as 
a parametrized curve 

{z(a,t) : a S K}, 

where 

z(a,t) = (zi(a,t), z 2 {a,t)) e R 2 . 

We suppose that all points (xi, x 2 ) £ R 2 with large positive x 2 belong to J7 1 (t) and that all 
(xi, x 2 ) with large negative x 2 belong to fl 2 (t). 

Let p t (x, t) and u l (x, t) denote (respectively) the pressure and fluid velocity in region fl l (t), 
for i — 1,2. Then the p l and u l satisfy the following conditions. 

1. divu 1 = in ft l (t), (i = 1, 2), i.e., the fluids are incompressible. 

2. = -Wp l ~ g(0, pi) in VL l {t) (i = 1, 2); this is Darcy's law HU) for fluid flow in a 
porous medium. (Here k is a physical constant associated to the porous medium, /Ltj and pi 
are physical constants associated to FLUID i, and g is the acceleration due to gravity. In 
particular, pi is the density of FLUID i and pi is the viscosity of FLUID i.) 

Moreover, on the interface, 

3. p 1 = p 2 
and 

4. u 1 — u 2 is tangent to the interface. 

Finally, the interface moves with the fluid, i.e., 

5. dtz(a,t) — u 1 (z(a,t),t) + ci(a,t)d a z(a,t). 

Here the function ci(a,t) may be chosen arbitrarily; it affects only the parametrization of 
the interface, which has no intrinsic physical meaning. (We could have used u 2 in 5, in 
place of u 1 .) 

See, e.g., Q where the Muskat problem is interpreted in terms of weak solutions. 
We may study two scenarios, namely 
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ASYMPTOTICALLY FLAT INTERFACES: lim Q ^ ±00 [z(a, t) - (a, 0)] = 
and 

PERIODIC INTERFACES: z(a, t) - (a, 0) is 27r-periodic in a for fixed t. 

We restrict attention to the case of two fluids with the same viscosity but different densities. 
In that case, /i 1 = fi 2 but p\ ^ p 2 in 2. 

Equations 1,...,5 take pressure and gravity into account but neglect surface tension. 

Once we know the interface z(a,t), we can derive explicit formulas for u 1 , u 2 , p 1 and p 2 
from 1,...,4 and elementary potential theory. Substituting these formulas into 5, and making a 
convenient choice of the function C\ (a, t), we obtain the Muskat equation for z(a, t). 

For asymptotically flat interfaces, the Muskat equation takes the form 



d t z,(a,t) = ^^P.V. I 1 J^l' [d a z» (a,t) - fyz M (J3, t)]d(3 ( M =l,2), 



2vr ' \z{a,t)-z(f3,t)f 

(1.1) 

where "P.V." denotes a principal value integral (the integral does not converge absolutely at 
infinity). 

For periodic interfaces, the Muskat equation is 

- / ,,_P2-Pi f sin(zi(a,t) - z 1 (/3,t))[d a z f _ l (a,t) - dpz^^t)} _ 

OtZn\C£, l) — — / — r — — — rr r — — — rzClp [Li — 1, A). 

2tt y R/27rZ cosh(z2{a,t) - z 2 {/3,t)) - cos(zi{a,t) - Zi(/9,*)) 

(1.2) 

Given an initial interface z°(a) = (zi(a), z® (a)) (a € W), we try to solve jl.l) or ( |1.2| i, with 
the initial condition 

z(a,0) = z°(a) (aUaeM). (1.3) 

There is no significant difference between the two flavors (JTTTJ, ( |1.2| i of the Muskat equation. 
Except for the introduction, we will work in the periodic setting. For the moment, we work with 
asymptotically flat interfaces, since the relevant formulas are then a bit simpler. 

Suppose first that the interface is the graph of a function, x 2 = f{xi). Then it is well-known 
that the initial-value problem (JTTTJ, ( |1.3| l is linearly stable for positive time, if and only if the 
heavy fluid lies below the interface and the light fluid lies above it (i.e., p 2 > pi). If instead 
the heavy fluid lies above the light fluid (pi > p 2 ), then the initial-value problem ( |1.1[ ), ( |1.3| l is 
linearly stable for negative time. See [7|. If the interface is not the graph of a function, as in 
Figure [T] then the heavy fluid lies below the light fluid locally near point A, and the reverse holds 
locally near point B. 

In that case, the initial-value problem ( |1.1[ ), jl.3\ is linearly unstable, both for positive time 
(thanks to point B), and for negative time (thanks to point A). 

Therefore, one asks whether an interface may begin in a (presumably) stable configuration as 
in Figure [2] and then "turn over" to a (presumably) highly unstable configuration as in Figure [3] 

For "very small" initial data, Constantin-Pugh Q, Yi Il20l . Siegel-Caflish-Howison 03), 
Cordoba-Gancedo [7] and Escher-Matioc |8|, prove that such a "turnover" cannot occur. It is 
shown in |3 1 that a turnover cannot occur if the initial interface has the form 

{{x 1 ,x 2 ) : x 2 = fixi)} 

with | /' | < 1 everywhere. 

Nevertheless, Castro et al. fl], show that a turnover from Figure|2]to Figure[3]can occur. 
Moreover, the Muskat solution exists and remains smooth for some time after the turnover, despite 
the instability of interfaces that look like Figure [3] 

The reason why the Muskat solution continues to exist in a highly unstable regime is that, 
by the time the turnover occurs, the interface has become real-analytic. In fact, starting from a 
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LIGHT FLUID 




Figure 1 : Non-graph interface. 




Figure 2: Initially stable interface. 

smooth, stable configuration as in Figure |2j the Muskat solution becomes instantly real-analytic. 
For t > 0, the function z(a,t) defining the interface dft 1 ^) = dVL 2 (t) = {z(a,t) : a € K} 
continues analytically to a strip {a E C : |3a| < h(t)}, with h(t) > ct for small t. As £ 
approaches the turnover time, the function h(t) may decrease, but it remains strictly positive up 
to the moment when the turnover occurs (i.e. up to the appearance of a vertical tangent to the 
interface). 

Therefore, the variants IfTTI . |[T2l of the Cauchy-Kowalewski theorem show that a real-analytic 
Muskat solution continues to exist for a short time after the turnover. The interface becomes 
(presumably) more and more unstable as the turnover progresses. For details see Q. 
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LIGHT FLUID 




HEAVY FLUID 



Figure 3: Turnover. 



LIGHT FLUID 




Figure 4: Breakdown. 



In (following the scheme of lfT3Tl for a scenario of two fluids with different viscosities and 
the same density) the construction of global-in-time solutions leads to ill-posedness (in Sobolev 
spaces in the Hadamard sense) for the unstable regime by going backwards in time. In |j 8 1 and 
|fl9l it was proposed, in the unstable regime, a relaxation of the phase distribution that leads 
to a well-posed problem and captures a mixing profile. From a completely different approach, 
a convex integration framework, unstable weak solutions are constructed for porous media in 

D2HHEI. 
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It is natural to believe that at some finite time T, the Muskat solution will break down. How- 
ever, the technique in H], is not strong enough to prove that such a breakdown occurs. 
In this paper, we show that a breakdown occurs, by proving the following result. 



Theorem 1.1 (Main Theorem). There exists a solution z(x,t) of the Muskat equation {1.2\, 
defined for t € [to, t?\, a £ R/2-7rZ, such that the following hold: 

A. At time to, the interface is a graph. 

B. At some time t\ G (to, t^), the interface is no longer a graph. 

C. For each time t G [to, tz), the interface is real analytic. 

D. At time ti, the interface is C 3 smooth but not C 4 smooth. 

Thus, at time ty,, there is no longer a C 4 solution of the Muskat equation. 

Note that (D) includes the assertion that the curve a — > z(a, t%) cannot be made C 4 -smooth 
by reparametrization. It is not enough merely to show that z(a, £2) ^ Cf oc . In our proof of the 
Main Theorem, the failure of C 4 smoothness in (D) will occur at a single point, as shown in 
Figure [4] 

To start to understand why Muskat solutions behave as described above, let us linearize equa- 
tion (JTTTJ about the trivial solution z(a, t) = (a, 0), a stationary horizontal interface. Perturbing 
z{a, t) to z(a, t) + Sz(a, t), and keeping only terms up to first order in Sz in (1.1) , we obtain the 
linear equation 



dt[5zix[a)] = Pl-flpy. r d a [5z»(a,t)]-de[5z,(M} d p (L4) 



with /i = 1, 2. 
The operator 

f^AfW :=--P.V.r fix) ~ fiy) dy (1.5) 

is the square root of the Laplacian in one dimension, i.e., A/(£) = (all £ € M), where ^ 

denotes the Fourier transform. 

The equation (1.4) takes the form 

d t F(x,t) = {Pl 2 ft) AF(x,f) (1.6) 
F(x,0)=f(x) (given). 

The operator A has a well-known connection with harmonic functions (see (TBI ). In fact, let 
u(x + iy) denote the Poisson integral of f(x), defined for y > 0. Similarly, for each fixed t, let 
U(x + iy, t) denote the Poisson integral of F(x, t). Then (1.6) is equivalent to the initial-value 
problem 



2 

U(x + iy, 0) —u(x + iy), 



8 t U(x + iy, t)=- yH \ H2J d y U(x + iy, t) 



which (formally) admits the solution U(x + iy, t) = u(x + i[y- te^-t]). 
In particular setting y = 0, we find that 

F(x,t)=u(x- V Pl ~ p2 ' t). (1.7) 
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This is the formal solution of our linearized Muskat equation ( |1.6| l. 

It is now obvious that the sign of p% — p2 plays a crucial role, and that solutions are real- 
analytic in the stable case. Indeed the harmonic function u(x + iy) is real-analytic in the upper 
half-plane y > but undefined in the lower half-plane y < and unless / is real-analytic. 
Therefore, if (p2 — pi)t > 0, then ( |1.7[ ) yields a real-analytic function of x; more precisely, 
F(x, t) continues analytically to the strip < <,p2 ~ Pl * > t. On the other hand, if (p2 — pi)t < 0, 
then (|1.7|> makes no sense unless the initial datum f(x) continues analytically to the strip \^sx\ < 



Note that the same "Rayleigh-Taylor" condition that predicts stability as time flows forward 
also predicts instability as time flows backward. 

More generally, if we linearize the Muskat equation about any given solution z(a, t) = 
(z\(a, t), Z2(a, t)), we arrive at a linearized Muskat equation of the form 

F t (x,t) = —a(x, t)AF(x, t) + HARMLESS TERMS, 
F(x,0) = f(x), 

where 

/ ,x (P2-P1) d x z 1 (x,t) . „ 

aM = (d xZlM y + (d^tw the R ^ h - T ^ 01 ' functlon - 

We therefore expect that our Muskat solution will be real-analytic for t > provided a(x, t) > 
everywhere, and highly unstable or non-existent for t > if a(x, t) < everywhere. 

If we pass from initial conditions at time zero to solutions at times t < 0, then the roles of 
cr > and a < are reversed. 

Observe that (assuming p2 > pi) the function <r(x, t) is positive if and only if the interface 
{zi(a, t), Z2(a, t)) is the graph of a function. In particular the significance of a turnover is that 
the Rayleigh-Taylor condition cr > for linear stability breaks down. 

In 13, we made a rigorous analysis of the full nonlinear problem and established analytic 
continuation of Muskat solutions to the time-varying strip {\^sx\ < h(t)}. This allowed us to 
construct real-analytic Muskat solutions z°(x, t) defined for t £ [-T, T] that "turn over" at time 
0. Here, we prove our Main Theorem |l.l| by constructing Muskat solutions z(x, t) analytic on a 
carefully chosen time-varying domain il(t) = {\?sx\ < h($lx,t)}. 

We take as initial datum the interface z(x,t) at some small positive time r, and solve the 
Muskat equation backwards in time until we reach time — r 2 . We take our initial datum z(x, r) to 
be a small perturbation of z°(x, t), where z° is the real-analytic "turnover" solution constructed 
in 12. 

Thus, we obtain Muskat solutions z analytic on tt(t) and close to z° in a suitable Sobolev 
norm. 

For all times t g [— r 2 , r), we take the domain f2(i) to be a neighborhood of the real axis, as 
in Figuref5] However, at time t = t, we take to be a domain that pinches to a point as shown 
in Figurelq 

In particular an initial datum z(x, r) may continue analytically to fi(r), yet fail to belong to 
C 4 as a function of a; € ML We take z(x, r) to belong to C 3 but not C 4 . Solving the Muskat 
equation (backwards in time) on 0(f) (t 6 [— r 2 , r] ) for such an initial datum, we obtain a Muskat 
solution as in the statement of our Main Theorem |1.1| In particular, assertions (A) and (B) of the 
Main Theorem hold because z is a small perturbation of z°, which is known to be a graph at time 
to = — t 2 , but not at time t\ = t 2 . (Recall that z° turns over at time zero.) Assertion (C) holds 
because z(x,t) continues analytically to Q(t), with fl(t) as in Figure |5]for t less than t2 = r. 
Finally, assertion (D) holds because we picked our initial datum z(-.t) to be in C 3 but not in C 4 . 
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Figure 5: Domain Q(t) at time t < r. 




Figure 6: Domain f2(t) at time t = r. 



8 



This concludes our summary of the proof of our main theorem. We warn the reader that the 
summary is somewhat oversimplified. The correct version of our proofs appears in the sections 
below. 



2 Preliminary estimates 

This section is devoted to proving crucial estimates which will be used in sections 4 and 5. 

2.1 Basic Properties of A 

Let T = R/2irZ and let / £ C 2 (T). For x € T we define 

Af(x) = - 1 [ \cot ( X ^]{f{x)-f{y))dy. (2.1) 

For x € T, let 

||z|| = distance(x, 27rZ). (2.2) 
Then a carefully justified integration by parts yields 

A/(x) = ~lim / esc 2 ( )(f( x ) _ f(y)) d y (2 .3) 



4-7T e^O 

\\x-y\\>e 

where the limit converges uniformly in x. Consequently, 

f(x)Af(x)dx = i-lim jj /(x) esc 2 ( ^ ) ( / ( ,■ ) /(.,/ ) U'll'l-r 

||X-J/||>E 

L lim // (7(x)-7( y )) C sc 2 f^ /N )(/(x)-/(y))dydx 



||x-y||>£ 



f(x)Af(x)dx = ±- / / esc 2 f H ) |/(x) - /(2/)| 2 dy^. (2.4) 



i.e., 

/ f(x)Af(x)dx = 

More generally, let A G C 2 (T), and suppose that \A"\ < C on T. We write C", C", etc. to 
denote constants depending only on C. 
Then the expression ( |2.3| l gives 

4(:r)7(z)A/(s)dx = ±- lim / / A(a;)7(a0 esc 2 f ^) (/(*) - /(y))dydi 



||x-j/||>s 

= ^lim jj (A(x)7(x)-A(y)J(y))c S c 2 (^-)lJir > -funuh/d: 

||x— y||>e 
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~lim ff (A(x)(f(x)-f(y)) + (A(x)-A(y))f(y))c S c 2 (^)iJU ) -/(,,)),/,„/,• 



\\x-y\\>e 



A(x)csc 2 ( ^-1) \f(x) - f(y)\ 2 dydx. 



J Jx,y£T V 2 

+ r- / 7(y) ( / ( A (x) - Aw)) esc 2 (/(*) - /(?/) k/.r }> ./// (2.5) 



8tt ./j/gt Uxer V - 



The second term on the right in ( |2.5| > is bounded by CH/Hj^a^, by Calderon's theorem on the 
commutator integral (see [16|). Therefore (|2.5|l gives 



A(x)f(x)Af(x)dx 



1 

x,y£T 



i . A{x) esc' I ) |/(x) - f(y)\ z dxdy + Error, (2.6) 



where 



|£rr r|<C'||/||f 2(T) 

provided \A'\ < ConT. 

Remark 2.1 Afofe fnaf the formula \2.6\ holds for complex-valued A and f. 



2.2 Contour Integrals Analogous to A 

Let h be a positive function in C 2 (T) (h > 0), such that \h\, \h'\, \h"\ and \h"'\ < C. In this 
section, C, C" etc. denote constants depending only on C. 

Let r + = {x + ih(x) : x € T} and r_ = {x - £/i(x) :ieT}. Here, both T + and T_ are 
contours in T + iR, oriented so that x increases as x ± i/i(x) moves along T±. 

Let ^(z) be 27r— periodic in z, and analytic in the region between T + and T_, and suppose 
F extends smoothly to the closure of that region. We write f+(x) = F(x + ih(x)) and f-(x) = 
F(x-ih(x)). 

We study the contour integral 



A r+ F(z) = -i J Icot(^)(F'(z)-F>))^, 



(2.7) 



for 2 e r + . 



We write z = x + ih(x) and w = u + ih(u). Hence, in the expression (2.7 1 we have 

+ ( 



F'(z) = (1 + iti {x))- 1 f' + (x) and = (1 + iti (it)) -1 fUu). Also, dw = (1 + ih'(u))du. 



Next, we show the following lemma: 
Lemma 2.2 The function 

\ 1 { (x + ih(x)) - (u + ih(u))\ 1. . (x — u 

A(x,u) = - cot o _ o( 1 + *' 1 cot 



2 V 2 y 2 V wy V 2 

is a C 1 function of (x, w) G T x T w/fn C 1 norm < C". 



10 



Proof: Since A(x, u) is 2tt— periodic in u for fixed x, we may suppose that \x — u\ < tt. In the 
region c < \x — u\ < tt (for small c to be chosen below), we have that cot((a; — u)/2) is C 1 with 
C 1 norm < C', and same holds for cot(((a; + ih(x)) — (u + ih(u)))/2) since 

dist((x + ih{x)) — (u + ih(u)), IttTL) > dist{x — u, 2tt%~) > c. 

In the region \x — u\ < c for small enough c, we have 

\(x + ih(xj) - (u + ih{u))\ < — , 

and therefore 

1 f(x + ih(x)) - (u + ih(u))\ 1 

t: cot — 7 : , : , , = ola;, u), 

where the C 1 norm of u) is bounded a priori in the region we are considering. Since also 
cot((a; — u) /2) /2 differs from l/(x — u) in \x — u\ < c by an error whose C 1 norm is bounded 
a priori, lemma ( |2.2[ i follows easily. 

Substituting the above results into the expression (2.7 1, we learn that 

A r+ F(x + ih(x)) = -- / {(l + i/^uJJ-^cot f ] +A(x,it) 

t J«eT I 2 V 2 / 



x {(l + t/i'ta:))- 1 /;^) - (l + iV(«))- x /+(«)} (l + iA'Cu))* 
If 1 f x — u 



i cot {(1 + ^'(a;))" 1 /;^) - (1 + ^»r7»} 

(1 + u) {(1 + ^'(a;))- 1 /!^) - (1 + ^») _1 /|00} dw 

(1 + to'to)" 1 {-~ j 2 cot (*^) (f' + (x) - /») J 



" 27 J ueT + ~ (1 + ^' W)" 1 ) COt ( V) } ^ (ti)dU 

— (l+tV(ar)) _1 /5-(a!) / (l + iti(u))A(x,u)du+ ~ [ A(x,u)f + (u)du. 
71 JueT n JueT 

On the right hand side, the first term is (1 + i/i'(x)) _1 )A/ + (a;); the expression in curly brackets 
in the second term has C 1 norm bounded a priori, thus allowing us to integrate by parts. We can 
perform a similar argument with the last term. Therefore 

A r+ F(x + ih(x)) = (1 + ih'(x))- 1 Af + (x) 

--{l + ih'ix))- 1 /'^) j (l + ih'(u))A(x,u)du + Error(x), (2.8) 
n JueT 

with 

\Error(x)\ < C||/+|| i2(T) . 
We study the A— integral in the expression ( |2.8| >. By definition, 

A(x, u) = i cot (^^) ~ (1 + ifc'fo))- 1 ^ cat (^^) . 
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with z = x + ih(x) and w — u + ih(u). Moreover, dw = (1 + ih'(u))du. Therefore, 

A(x.u)(l + iti(u))du = 



P. V. I 1 COt ( dw - P.V. f 1 COt ( da. 



- + 

with z = x + ih(x). Since 



wer + 2 \ 2 J Just 2 V 2 



we have that 



1 / x — u . 
/>.K.y -cat (—-)./„ 0. 



A(x,u)(l + ih'(u))du = P.V. I cot ( Z —^- ) dw = 0. (2.9) 



The last equality in (|2.9[) is an exercise in contour integration. 



Therefore, the expressions ( |2.8[ > and (2.9 1 yield the following result, 

A r+ F{x + ih{x)) = (1 + i/i'(.T))- 1 A/ + (x) + Error + {x), (2.10) 

where 

|| J Brror + || L 2 (T) < C||/ + || L 2 (T ). 

Similarly, 

Ar_F(z) = -^j r ^ cot (^) (F'W - F'H)d w , (2.11) 
for z£L satisfies 

A r _F(z) = (l-i/i'(a;))- :L A/_(a;) + Brror_(a;), (2.12) 

where 

||£;rror_|| L 2 (T) < C||/_|| L 2 (T) . 



Next we return to the definition ( 2.7 i, and shift the contour from T + to T_ (note that the integrand 
in the expression ( |2.7) has no poles). 
Thus, 

A r+ FO + i/i(:E)) = --* y i cot ^^Vf'(z) - F'(w))dw 

for z = x + ih{x). 

Another exercise in contour integration shows that 



r_ 2 



I f z — w . 
cot — - — dw = — ^7^, 
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and therefore, 



A r+ F(x + ih(x)) = iF'(x + ih{x)) + 



2n 



cot 



iF'(x + ih(x)) 



1 

4tt 



CSC 



F(w)dv 



F'(w)dw 



(2.13) 



For z = x + ih(x) 6 T + , w = u — ih(u) G T_, we have 

C 



CSC 



< 



C 



< 



C 



Hence, 





f z — w\ 


/ esc 2 


V 2 ) 







z - w)\\ 2 + |9f(z - w)\ 2 \\x -u\\ 2 + (h(x) + h{u)) 2 
, \F(w)\ = |/_(«)|, < Cdu. 

!/-(«)! 



||x-u|| 2 + /i(a;) 2 

F(w)dw 



< C 



- zT \\x~u\\ 2 + h(xy 



du 



<Ch- 1 (x)M[f.}{x), 

where M [/_] denotes the Hardy-Littlewood maximal function of /_. (See Jl6 |, page 9). 
Together with ( |2.13| l, this implies: 



where 

Similarly, 

where 



Ar + F(x + ih(x)) — iF'(x + ih(x)) + Erreur + (x), 

\\h(-)Erreur + (-)\\ LHT) <C\\f-\\ L 2 {T) . 
Ar_F(x — ih(x)) = —iF(x — ih(x)) + Erreur-(x) 1 



(2.14) 



(2.15) 



\\h(-)Erreur-(-)\\ LHT) <C\\f + \\ LHT) . 
From ( |2.10| i and ( |2.14| i, we see that 

iF'(x + ih{x)) = (1 + i/&'(x)) -1 A/ + (a;) + Goof+{x) 

where 

\\h(-)Goof + (-)\\ LHT) < C||/ + || L2(T) + C||/_|| i2(T) . 

Similarly 

- iF'{x - ih(x)) = (1 - i/i / (a;))- 1 A/_(x) + Goof-(x) 

where 

\\h(-)Goof_(-)\\ L , m < C||/+|U2 ( T) + C||/-IU= (T ). 



(2.16) 



(2.17) 
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2.3 Applying results derived above 

In this section, we are given two functions h, h t £ C 2 (T) (with C 3 norms assumed bounded a 
priori). We again assume that h > and we take the contours T + = {x + ih(x) : x £ T}, 
T_ = {x — ih(x) : x £ T}, and let F(z) be a 2ir~ periodic analytic function on the region 
between T + and T_, smooth on the closure of that region. We take f+(x) = F(x + ih(x)) and 
f-(x) = F(x - ih(x)) for ieT. 

Let 6 be a function on T, with C°— norm assumed bounded a priori. 

We assume that 

\h t (x)\ < Sh(x) (2.18) 

for all x £ supp{9), where S is a given positive number. 

Let u i (x) be a function on T, assumed to have a priori bounded C 2 norm. We want to estimate 

X ee » [ J + (x)9(x) l-*Mf 1 cot ( {X + IH{X)) - {U + lh{u)) 



x(F'(> + i/i(a;)) - F'(u + + ih'(u))du 

+ih t (x)F'(x + ih{x)) | dx, (2.19) 

from below. Applying ( |2.10| i and ( |2.16| >, we have 

g i (x) f 1 / (x + ih(x)) — (u + ih(u)) \ , . , , „, . , , , , , . 

' cot — - — — x(i? ,/ (a;+i/i(ir))-F / (w+i/i(u)))(l+i/i / (M))du 



with 

and 

with 



ttsT 2 \ 2 j 

= ~o x (x){\ + ^/^'(x))- 1 A/ + (a;) + Err+{x), 

\\Err + \\ L 2 (T) < C\\f + \\ L 2 (T): 

ih t (x)F'(x + ih(x)) = h t (x)(l + ih'{x))- 1 Af + (x) + Oops + {x), 

\\OopS + \\ L 2 {supp(g)) < CS (||/ + ||i2( T ) + ||/-||i a (T)) ■ 



Therefore, 

X = n J + (x)6{x){l + iti(x))- 1 {a 1 (x) + h t {x)}Af + {x)dx + Discrep +1 (2.20) 
where 

\Discrep + \ < C(S + 1) (||/+||£ 3(T) + \\f-\\h m ) ■ 

We recall that the derivatives of (1 + ih'(x))^ 1 up to order two are bounded a priori. We make 
the following assumption: 

The C 2 norm of 9{x){a\(x) + h t (x)} is bounded a priori. 

Then, by (2.6 1 we have 

X=^JJ^ {0(x)(l + ih'ix))- 1 ^) + h t (x))} esc 2 |/+0«0 - f+(y)\ 2 dxdy 

+ Errore+ (2.21) 
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with 

\Errore+\ < C(S+l) (||/+||| 2(T) + ||/-||! 3(T) ) . 

If 

9? {9{x)(l + ih'ix))- 1 ^!^) + h t {x))} > 0, (2.22) 

then §221\ yields 

X > -C(S + 1) (||/ + ||£ 2(T) + ||/_||! 2(T) ) • (2.23) 
Thus we have proven the following result. 

Lemma 2.3 Let 9, d\, h, ht be functions on T. Let S be a positive number. We make the following 
assumptions 

1. h>0 on T. 

2. \h t \ < S h on supp (9). 

3. The C 3 norms of h and h t are bounded a priori. 

4. The C 2 norm ofd\ is bounded a priori. 

5. The C 2 norm of 9{x)(l + ih! {x))~ 1 (5\{x) + h t (x)) is bounded a priori. 

6. 5ft {9{x){l + ih'(x))- 1 (a 1 (x) + h t {x))} > on T. 

7. The C°-norm of 9 is bounded a priori. 

Let F(z) be analytic on {z G C : < h(3i(z))}, 2ir— periodic and smooth on the closure 

of the above region. Let 

f+(x) = F(x + ih(x)) 

and 

f-(x) = F(x-ih(x)), 
for i£T. Then the following inequality holds: 

*[ j +mx) L*M[ i cot {( x + Mx))-(u + Mu))\ 

JxeT I 77 Ju&t 2 V 2 / 

x(F'(x + ih(x)) - F'(u + ih(u)))(l + iti {u))du 
+ih t (x)F'(x + ih(x)) j dx 

>-c(s+i) (||/ + ni 2(T) + ||/-||| 2(T) ). 

Similarly, under assumptions analogous to 1,...,7, with 1 — ih! {x) in place ofl + ih! {x) we have 

ftf j(x)9(x)L^[ i a* -ihM)\ 

x{F'(x - ih{x)) - F'{u - ih(u)))(l - ih'{u))du 
—iht{x)F'(x — ih(x)) j dx 

>-C(5 + l)(||/ + ||i 2(T) + ||/_||| 2(T) ). 
Here, the constant C depends only on the a priori bounds assumed in 3, 4, 5 and 7 above. 
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2.4 Garding Inequalities 

First we work on K. Later on, we will switch over to T. 

Here Ar will be the square root of (—A) in one dimension, D = \ on E and 

U a : L 2 {R) -> L 2 (R) 

for a — ±1 is the orthogonal projection onto either positive frequencies (er = 1) or negative 
frequencies (er = —1). Note Ii a = \ (I + iaH), where H is the Hilbert transform (H = 
-isign($,)). 

Commutator estimates: 

Let a be a real-valued function on ML We suppose a' 6 L°°(R), and write a for the operator 
/^ a -/onL 2 (IR). Then 

[a,H]D : L 2 (R) -> L 2 (R) 

is bounded, with norm at most C||a'||£oo( R ), where C is a universal constant. 
That follows from Calderon's theorem on the commutator integral (see 1 16]). 
Since A R = iDH, it follows that 

[a,H]A m : L 2 {R) -> L 2 (R) 
is bounded, with norm at most C||a'||£oo( R ). 
Since ILj- = |(7 + iaH), it follows that 

\\[a,Il a }D\\, ||[a, n CT ]A R j| <C||a'|| i0 o (R) . 

Next, recall that 

25R(7-A R /) >A r (|/| 2 ) 

pointwise (see [5 1). 

If a > everywhere and Aa 6 L°°, then 

»{ / a(x)/(x)A R /(x)dx} > l -f a{x) (A R (|/| 2 )) cfe 
= i / (A R a)(x)|/(x)| 2 dx>-i||ARa||^||/||| 2 . 

z JR * 

Now suppose a, 6 are real, a', b' € A R a, A R 6 £ L°°; a > |6| pointwise. 
Note that A R = £<r=±i II^ArII^ and D = £ CT=±1 n CT £»n CT . 
Therefore, 



aA R + 6£> = ^ {[a,Il a }A R n a + [b,n a }Dn a } 

<T = ±1 

+ ^n CT (aA R + 6D)n ff . 

The terms in curly brackets are operators with norm at most C\\a'\ +C||6'| on L 2 (R), 
with C a universal constant. 
Therefore, 

K((aA R + bD)f, /) =En-or + ^ 5J(II CT (aA R + bD)U a f, /) 

(T 

=En-or + 5R((aA R + 6£>)n ff /, U a f) 

a 

=En-or + ^ 5R((a + c^ArIT,/, !!„/), 



16 



since DIi a = ctArILj; here 

|Error|<C'(||o'|| i ~+||6'|| i ~, lullL2 
Recalling that a + ab > on R, we see that 

R<(aA R + 6D)/,/) > -C{\\a'\\ L oo + \\tf\\ L oo + ||A R a|| L oo + IIAr6IIlo.HI/III,, 

where C is a universal constant. 

Now suppose 0,6 real, a > |6| on R, |o|, |6|, |o'|, |6'|, |o"|, \b"\ < 1. Then 



a:) = -P.V. 

7T 



Now 



/ < / 

^ix-j/i>i {x-yr J\x- y 



{x - y) 

|o(z)| + |o(y)| 



dy< C, 



and 



J|x-l/|<l 



apr) - a(y) 
(a; - y) 2 



P.V. / 

J\x-y\<l 



\>i (x-y)- 

a(x)-a(y) + a'(x)(x-y) 



(x - y) 2 



dy 



< C. 



Similarly for b. Therefore, a', 6', A K a, A R 6 are bounded on R. 
Therefore, we obtain the following result: 

Suppose a, b are real functions on R, such that \a\, \b\, \a'\, \b'\, \a"\, \b"\ < 1; and suppose 
a > \b\ everywhere on R. Then 

5R((oA K + 6£>)/,/)>-q|/||| 2 

for / e Cq (R), where C is a universal constant. 
Rescaling, we obtain the following result: 

Lemma 2.4 Let < 8 < 1, and /"ef a, 6 foe rea/" functions on R, swc/z ffoaf |a| < 5, |6| < 5, 
\a'\ < 1, |6'| < 1, |a"| < S^ 1 , \b"\ < S^ 1 everywhere, and also a > |6| everywhere. Then 

S((oA R + 6D)/,/)>-q|/||| 2 

/or / G Cq (R), where C is a universal constant. 

Next, we transfer our result from R to T. 

Lemma 2.5 Let < 8 < 1, and /ef a, 6 foe real-valued functions on T. Suppose that 



dx 



a(x) 



b(x) 



< S 1 ^ 



for j = 0,l,2 and x E T. 

Suppose also that |6| < a for all x £ T. Then, for any smooth function f on T, we have 

Myj(x)(a(x)Af(x) + b(x)if'(x))dx} > -C Jjf(x)\ 2 dx, (2.24) 



where C is a universal constant. 
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Proof: Suppose first that supp(a), supp(b) C [—1, 1]. We identify T with [— ir,ir], and let 
/ = on M \ [— 7r. 7r]. Thus, we may regard / as a function on BL For x £ [—1, 1], we have 



Af(x) - A R /(x) = / K(x-y)f(y)dy 



for a kernel if, where | |if 1 



< C. 



Hence, in proving ( 2.24 1, we may replace A by Ar. Estimate ( 2.24 1 now follows from lemma 



We have established ( 2. 24} in t he case where supp(a), supp(b) C /, where / is any interval 
of length 2 in T. Conclusion (2.24 1 in the general case now follows by using a partition of unity. 



2.5 Applications of Garding's Inequality 

In this section, we are given functions h, h t £ C 3 (T) (with norms bounded a priori). As usual, 
we assume that h > 0, and we take T + = {x + ih(x) : x £ T}, T_ = {x — ih(x) : x £ T}. We 
let F be analytic on the region of T + iR bounded by T + and T_ and we suppose that F is smooth 
on the closure of this region. As usual, we set f+(x) — F(x + ih(x)), /_ (x) — F(x — ih(x)). 



Let &i(x) and 6{x) be functions on T. As in section (2.3 1, we want a lower bound for 

ai(x) 



X ee 



f + {x)6{x) 



1 fix + ih(x)) - (u + ih(u)) 

o C0t I o 

ueT 2 V 2 



IxGT L 7r 

x(F'(x + ih(x)) - F'(u + ih(u)))(l + iti {u))du 

+ih t {x)F'(x + ih(x)) j dx. 

We assume that the C 2 norm of b\ is bounded a priori and that < 9 < 1 on T. According to 
(fZ7]i and (fTTOjl we have 



(2.25) 



X = ft I f + (x)9(x) {criO)(l + ih'{x))- l kf+{x) + ih t {x)F'{x + ih(x))} dx + Err, 



with 
Also, 

as noted before. Therefore, 



\Err\ < C||/ + ||| 2(T) . 
F'(x + ih(x)) - (1 + ihix))- 1 /'^), 



X = ft I f + {x)e{x){l + ih l {x))- 1 {d l {x)Kf + {x) + ih t { : x)f + {x)}dx + Err, (2.26) 



with 



\Err\ < C\\f. 11 J 



+ Mi 2 (T)- 



We now suppose that: 

The C 2 norms of 9{x)&i{x) and 9(x)h t (x) are bounded a priori. 

Then the same holds for 9(x)(l + ih' (x^a^x) and 9(x)(l + ih'{x))- l h t (x). 
Therefore, formula (|2.6|l gives 



(2.27) 



K / f + (x)9(x)(l + ih / (x))- 1 a 1 {x)Af + dx 
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with 



Also, if iB(x) is real, then 



/ f + {x)3i{9(x){l + ih l (x))- 1 a 1 (x)}Af + (x)dx + Err, 



\Err\ < C\\f + \\l m . 



/ f + (x)B(x)if + (x)dx 
JxeT 



iB{x)\±\f + {x)\ 2 dx 



£CGT 



iB(x)M(f + (x)f + (x)) 
B'(x)\f+(x)\ 2 dx. 



Applying this with B(x) — i^s{9(x)(l + ih'(x)) 1 h t (x)}, we see that 



® / f+(x)0{x){l + ih'(x))- 1 h t (x)if + (x)dx 
= 3?/ f + (x)^{e(x)(l + ih l (x))- 1 h t (x)}if' + (x)dx + Err, 



with 



\Err\ < C\\U\\l 2{ 



T) ' 



Putting ( |2.28| l and ( |2.29| l into ( |2.26| ) we see that 

X = 3? / / + (ar) {A(ie)A/ + (» + B(a;)i/|(x)} + £rr, 



where 
and 

We assume that 

for < j < 2; and that 



\Err\ < C||/ + ||2 2(T) 

A(x) = $l{6(x)(l + ih'{x))- l ax{x)}, 
B(x) = ft{9(x)(l + ih'{x))- 1 h t {x))}. 



< OS 1 '* 



\B(x)\ < A(x) for all x e T. 













ay™ 



From (2.30 1 and Garding's inequality we obtain the estimate 

X>-C||/ + ||1 3(T) . 
Thus, we have proven the following result: 



(2.28) 



(2.29) 



(2.30) 



(2.31) 
(2.32) 

(2.33) 

(2.34) 

(2.35) 



Lemma 2.6 Let 9, h, h t , cf\ be functions on T, and let < 5 < 1 be a number. We make the 
following assumptions: 

1. h>0. 

2. The norm of a\ in C 2 and the norms ofh and h t in C 3 are bounded a priori. 

3. The C 2 norms of 9{x)d\(x) and 9(x)h t (x) are bounded a priori. 
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supp 8 c [—5, 5]. 



5. 



[£Y R{0(a;)(l + ih'{x))- 1 a l {x)} , (±y $l{9(x)(l + ih'(x))- l h t (x)} < 
for < j < 2, x e T. 

6. \&{0(x)(l +ih , (x))- 1 h t (x)}\ < $t{6(x)(l + ih '(x))' 1 ^)} for all x e T. 

7. o < e < i. 

Then 



If 



CTi(x) 



1 



cot 



(.t + ih(x)) — (u + ih(u)) 



x(F'{x + ih(x)) - F'(u + ih(u)))(l + ih'(u))du 
+iht(x)F' (x + ih(x)) j dx 
> ~C"II/+IIl2(t)2! 

where f+{x) = F(x + ih(x)). 

Here, C depends only on the a priori bounds in 2, 3 and on the constant C in 5. 



Remark 2.7 An analogous result holds for f-(x) — F(x — ih(x)). 



3 The unperturbed solutions z and the height function h 

In this section we study the Rayleigh-Taylor function associated to the unperturbed solutions. 
We construct our time-varying domain of analyticity f2(t) = {\^sx\ < h(^tx, t)} adapted to the 
Rayleigh-Taylor function. 



3.1 The Rayleigh Taylor Function 

We consider the analytic solution z(x, t) = (z^x, t),z 2 (x, t)), for t € [— T, T] (T > 0), found 
in (21 ■ We will call z(x, t) the unperturbed solution and we recall that we can choose the origin 
of time so that z(x, t) satisfies the following properties: 

• Z-i( x i t) — x and z 2 (x, t) are 2tt— periodic functions and real for x real. 

• Zi(x, 0) and z 2 (x, 0) are odd functions. 

• 2 1 (0,0) = 0, 

• 9,^(0,0) =0, 
. 9^(0,0) = 0, 

• 9^ (0,0) > ci > 0, 

• d t d x z x (0,0) < -c 2 <0. 

• 9^(0,0) > c > 0. 

Then since \d^z x {x, 0)| < C3 it follows from Taylor's theorem that 

c 4 x 2 < d x z-y{x,Q) < C 5 x 2 , (3.1) 

for |x| < c 6 . 

Similarly, since 

9*9^(0,0) < -02 < 
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and 

19^(^0)1 <C 7 

it follows that 

d t d x z x (x,Q) < -c 8 < 

for \x\ < Cg. 
Since also 

\^d a Zi{x,t)\ < C w 
for \x\ < eg, \t\ < en, it follows in turn that 

d x z x (x,Q) - Cut < d x z x (x,t) < d x z x {x,Q) - ci 3 £ (3.2) 

for |x| < Cg, < t < C14. 

Combining our estimates (3.1 1 and (3.2 1, we find that 

C4X 2 — G\<zt < d x z 1 (x,t) < C$x 2 — ci^t (3.3) 

for |x| < C15, < t < Ci6- Similarly we have that 

c' A x 2 - c' 12 t < d x z x (x,t)< C' 5 x 2 - C' 13 t (3.4) 

for \x\ < C15, — ci6 < £ < 0. 

Also, since 9^^(0,0) = and z x (x, t) has C 3 — norm at most C on {|x| < c, |t| < c}, it 
follows that 

\dlz x {x,t)\<C l7 \t\+Cn\x\ (3.5) 

for \x\ < Ci 8 , |t| < Cig. 

We have also 

\d x z x (x,t)\ < C 20 (3.6) 

for \t\ < c. Finally we notice that 

d x z 1 (x, —t) > C2ot for all x E T and < t < C21. (3.7) 

Now we define the Rayleigh-Taylor function associated to the unperturbed solution as 

n/ f) _ -2ird x Zi{x,t) 
(Jl(X ' t) {d x z x {x,t)y + {d x z^x,t)Y KSX) 

Therefore cr x (x, t) satisfies the following: 

1. <Ji{-, t) is analytic on {x + iy : x E T, \y\ < c;,} with \a x (x + iy, t)\ < C, for all x + iy 
as above and for all t E [— r, r]. 

2. <Ji(x, t) is real for x € T, t € [— r, t]. 

3. crj has a priori bounded C 2 norm as a function of (x, !)eTx [— r, r]. 

4. of (0,0) = 0. 

5. 0^(0,0) =0. 

6. Sgaf (0,0) = -c 2 < 0. 

7. <Vj>(0,0) = ci > 0. 
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For complex x (|Sjx| < c'), we have 

<r°(x,i) = <T?(0,i) +xd x a° 1 (0,t) + id 2 a°((M)x 2 + (\\x\\ 3 ), 

where ||x|| = ||5ix|| + \^sx\ ss distance(x, 2nZ). 
Moreover, 

<r?(0,i) =c!i + 0(t 2 ), 

and we have 

#r?(0,i) = -c 2 + 0(t) 
9^(0, i) = 0(t). 

Hence, for complex x such that \x\ < 10~ 3 and \^sx\ < d we have 

<t?(x, t) = C!i + 0(t 2 ) + 0(si) - y a; 2 + 0(te 2 ) + 0{x 3 ). 
Since also |0(xi)| < j^x 2 + 0(t 2 ), we conclude that 

Kom) - [cit- f x 2 ]| < io- 2 [ Cl t + |x 2 ], 

if |x| <c,te [0,t]. 
In particular, 



ci(l - 10~ 2 )t- y(l + l(T 2 )x 2 < cr?(x,t) < c x (l + l(T 2 )i - y(l - 10~ 2 )x 2 (3.9) 

for real x, if |x| < c, i e [0, r]. Similarly, we have that 

4(1 + 10 _2 )i - y(l + 10~ 2 )x 2 < <rl{x,t) < 4(1 - l(T 2 )i- y(l - l(T 2 )x 2 (3.10) 

for real x, if |x| <c,te [— r, 0]. 

On the other hand, for the Muskat solution in [2 1, (|3.8|> gives 



a-?(x,0)<-c' if ^ < |x| < 7T, x real. (3.11) 
Then, by taking r small enough, we achieve 

4(x,t)<-- if ^ < |x| < tt,x real^G [0,r]. (3.12) 
If h(x, t) is a positive functions defined for x € T, < t < r, and if 

< /i(x, i) < -q, for all such (x, t) (see property 1 above), (3.13) 

then 

|cr?(x ± j/i(x, *),*)- o-?(x, t)| < Ch(x, t). (3.14) 

We will work with a perturbed Rayleigh-Taylor function <7i(x, t), defined and analytic on 
{|3x| < hQRx, t)} for each fixed t € [0, r]; and smooth on the closure of that region. 
We assume that 

|^( f 7 1 (x,t)-a°(x,t))| <r 10 (3.15) 

for < j < 2, and for |3x| < h(9tx, t),0<t<r. 

We have not yet picked the function h(x, t). This we do in the next section. 
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Figure 7: Red: Function h(x,t) at time t = t 2 , with k = 0. Blue: Function h(x,t) at time t = t, 
with k = 0. 



3.2 Two Height Functions 



In this section, we suppose we are given a Rayleigh-Taylor function crj as in the section (3.1 



We pick A » 1 and < r, ft, A << 1 real numbers. We first pick A large enough, then pick 
r small enough (depending on ^4). Having picked A and r, we then take < A < t 100 and 
< k < t 100 . We fix such A, r, A until the end of the paper. The parameter k will be held fixed 
until section 5.7 below, in which will let k tend to zero. From now on, we suppose A, r, A and k 
have been picked as above. 

For i£T and t 2 < t < r, we define 

h(x, t) = A-\t 2 - t 2 ) + (A- 1 - A(t - t)) sin 2 (|) + k. (3.16) 

See Figures (|7]i and (jS}. Thus, 

d t h{x,t) = -2A~H + Asm 2 (J^j . (3.17) 

See Figures ((^i and {lOl. If ||x|| > lOA^tK then sin 2 (x/2) > 10" 2 ||a;|| 2 > A- 2 t, hence 
yields 



and 



\d t h(x,t)\ < 2A~ x t + Asm 2 (^j for all x G T, t G [r 2 ,r], 
< 3Asin 2 (|) if ||a:|| > lOA -1 *^* G [r 2 ,r]. 



(3.18) 



(3.19) 



On the other hand, (3.16 1 yields h(x,t) > \A~ X sin 2 (a;/2) for all x G T, < G [r 2 ,r). Hence, 
(3. 19 1 implies 



|9 t /i(x,i)| <6A 2 /i(x,t) for||x|| > 10A -1 *5,iG [t 2 ,t] 



(3.20) 



Next, let ||x|| < | with c as in inequality (3.9 1. Applying that estimate we learn that 



a° 1 (x,t) + d t h(x,t) - A?h(x,t) 
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Figure 8: Zoom of Figure 1 around the point (0, 0). 




Figure 9: Blue: R-T function at time t = r. Red: Function dth(x, t) at time t = r. Green: Function 
dth(x, t) at time t = t with a larger value of A. 
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Figure 10: Zoom of Figure 3 around the point (0, 0). 
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> (c x f - C 1 sin 2 ^ + (-2A~H + Asm 2 (|))-^ [^(t 2 - t 2 ) + (A' 1 - A(r - t)) sin 2 (|) + 



cit - 2A~ x t - A-2(t 2 _ t 2 ) - A^n 



-Ci + A- (A~i - Ai(r-t))] sin 2 (|) 



The second term in square brackets is positive (since A is large but A 3 / 2 r is small), and the first 
term in square brackets is greater than cr 2 . 
Hence, 



1 



a"(x,t) +d t h(x,t) - A*h{x,t) > ct 2 if\\x\\ <-c,te [t 2 ,t 



Remark 3.1 The reason we restrict to t > r 2 is that, if we allowed t << t 2 , then the first term 
in square brackets above might be negative. 



On the other hand, sup pose that | \x\ \ > \c, t e [t 2 , t]. Then |crj(x,t)| < C, and (|3.17| yields 

d t h(x, t) > cA, while 



3.16 1 implies that < h(x,t) < CA . Hence, in this case 



o-° l {x,t) + d t h{x,t) - A?h(x,t) > ct 2 . 

Thus, in all cases, we have 

<Ti(x,t) + d t h(x,t) - A?h{x,t) > ct 2 for a; G T, t £ [t 2 ,t] 



(3.21) 



Next, suppose ||x|| < 20 A 1 t?, t 6 [t 2 ,t]. Then inequality (3.9i gives <j\{x,t) > ci, while 
< [3.16| ) gives < fc(ar,t) < CA _1 i and ( f3TT7] > gives |5 t /i(ar,t)| < CA -1 *. Hence, 



a%(x,t) - \d t h{x,t)\ - Aih(x,t) >c!t>c!r 2 for\\x\\ < 20A~ 1 t = , t G [r 2 ,r]. (3.22) 



Next, we investigate the x— derivatives of h(x,t), ht(x,t) and crj(x, f). From (3.16 » and (3.17i, 
we see that 

\dlh(x,t)\ < CA- 1 for i e T, ( e [r 2 ,r],0 < j < 3; 



and 
Also, 



|0£ftt(z,t)| < for i e T, ( e [r 2 ,r],0 < j < 3. 

\dia%(x,t)\ <C forxeXte [r 2 ,r],0< j < 2. 



(See section ( |3.1| l). 



Now suppose that \ \x\ \ < 20 A Y t a . According to estimate ( 3.9 1, we have 

ct<o\{x,€) <Ct for\\x\\ <20A- 1 tKte [r 2 ,r]. 



(3.23) 
(3.24) 
(3.25) 

(3.26) 



Since (3.26 1 holds for all x G [— 20A 1 ^2 ; 20A 1 i a ], the mean-val ue the orem implies that 
\d x a%(x,t)\ < CAti somewhere in [-20 A~H^ , 20A~H%}. Hence, by \?>.25\, we have 



(3.27) 



| d x al (x,t) | < CA^, for all ||a:|| < 20A~HKt G [t 2 ,t]. 
From ( f3~25] l, (J3T26J, ( f3~27] i, we see that 

|^cr?(a;,i)| < CA (i^ 3 for all | \x\ \ < 20A~H^ , t G [r 2 , t], < j < 2. (3.28) 
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Next, return to h. From ( |3.16| l, we have 

< h(x,t) < CA~H for all ||x|| < 20A~ 1 tK t G [r 2 ,r] 

and 

\d x h{x,t)\ < CA~Hi for all ||x|| < 20A~ 1 i5, i G [r 2 ,r]. 
From ( [3T23) ,( f3T29| >,( [33Ql >, we see that 



(3.29) 
(3.30) 



i \ 2-J 



|#j/i(x,£)| < C(f5) " for all ||x|| < 20A~HK t G [r 2 ,r], < j < 2. (3.31) 



Similarly, regarding dth(x, t), we learn from (3.17 1 that 

\d t h(x,t)\ < CA~H for all ||x|| < 20A~HK te [r 2 ,r]; 

and 

\d x d t h{x,t)\ < Cti for all ||x|| < 20A~ 1 ti, t G [r 2 ,r]. 
From ( |3.24| >, ( |3.32] > and ( |3.33| l we obtain the estimate 



(3.32) 
(3.33) 



|<9£<9 t /i(x,i)| < CA (t^ 2 3 for all ||x|| < 20A" 1 ii, t G [r 2 ,r], < j < 2. (3.34) 
It will also be useful to estimate d x <7i{z,t) for z = x + iy with ||x|| < 20v4~ 1 £2 and 

\y\<h(x). 

From the assumption 1 in section (3.1 1 and since < h < CA~ x t we have \di +1 a®(x + 
iy, t)\ < C in this range, hence 



\dl<yl{x + iy,t)\ < \dlal{x,i)\+C\y\ < CA (t*) ' + CA -1 * < CA 



i \ 2-J 



Thus, 



|^<T?(x + iy,t)| <CA (t*)' 
for z = x + iy, ||x|| < 20A _1 t3, f G [r 2 ,r], |y| < /i(x), < j < 2. (3.35) 



We look at h(x, r 2 ). From (3.16i, we have 



(3.36) 



h(x, t 2 ) > ^A^t 2 + ^A^ 1 sin 2 (|) for iel 

et 

H(x, t) = - (VV 2 + A" 1 sin 2 (|) ) + ^~ 2 r< + At sin 2 (-) for x G T, i G [-r 2 , r 2 ] 



Note that 



< h(x,T 2 ) < h(x,T 2 ) forallxeT. 

From p.37| i, we have 

dth(x,t) = A- 2 t + Asm 2 (|) forxeT, te [-t 2 ,t 2 ]. 



On the other hand, (3.37 i also gives 



1 



h(x, t)>-( A~ y t 2 + A- 1 sin 2 - for x G T, t G [-r 2 , t 2 ] 



(3.37) 
(3.38) 

(3.39) 
(3.40) 
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From ( |3.39| l, ( |3.40] > we obtain the estimate 

\d t h(x,t)\ < Cr~ x h{x,t) for all x £ T, t G [-t 2 ,t 2 ]. 
Next, suppose | < c; from ( |3.9| and ( |3. 10| > we see that 

cr?(x,t) >ct-Csin 2 (|) for<e[0,r 2 ], 
ai(x,t) >c't-C'sin 2 (|) fori G [-r 2 , 0] 

Therefore 

crj(x,t) + d t h{x,t) - A^Tifoi) > (ci - C sin 2 + (A~ 2 t + A sin 2 

- (a^t 2 + A" 1 sin 2 (|) ) + A" 2 rf + At sin 2 (|) 



(3.41) 

(3.42) 
(3.43) 



-A3 



i 9 1 * 1 9 3 

ct + A 2 t A~2T 2 - A~*Tt 

4 



-C + A- ^A-5 - Alt] sin 2 (|) , 



for t £ [0,t 2 ]. The first expression in curly brackets is greater than \A~ 2 t, and the second 
expression in curly brackets is positive. 
Therefore, 

CT?(x,t) +d t h(x,t) - A?h(x,t) > \a- 2 t if ||x|| <c,t£ [0,r 2 ]. 



On the other hand, suppose ||x|| > c. Then ( 3.39 1 gives dt%[x,t) > cA, whereas (3.37 1 gives 
< h(x,t) < CA- 1 , and we know that \a%(x,t)\ < C. Hence, 

ol{x,t) + dth{x,t)-A*h{x,t) >cA-C-CA~^ >c'A> \a~ 2 t if\\x\\ >c,t£ [0,r 2 ]. 



Thus, in all cases, 



a°(x,t) + d t h(x,t) - Aifi(x,t) > \a~ 2 t forxeT, t£ [0,r 2 ] 



Similar computations yield 



cr?(x, t) + d t h(x, t) - A^h(x, t) > ^A" 2 r for x G T, f G [-t 2 , t 21 



Regarding the smoothness of a®, h, dth, we see from ( 3.37 1 and ( 3.39 1 that 

\dih(x,t)\ <CA- X forieT, i e [-r 2 ,r 2 ], < j < 3; 

and 

|^5tfi(x,<)| <CA forzGT.tG [-t 2 , r 2 ], < j < 3. 
Also, from our assumption on a® in section (3.1 1, we have 



\dia1(z,t) \ < C for^z G T, \%z\ <c,t£ [-t 2 ,t], < j < 2. 



(3.44) 

(3.45) 

(3.46) 
(3.47) 

(3.48) 



Note that (3.48 i holds whenever t £ [-t 2 ,t 2 ] and \$Sz\ < h ($tz,t), since (3.37 1 gives < 



h(x, t) < CA" 1 for all x £ T, t £ [-r 2 , r 2 ]. Similarly, (3.48 i holds whenever t £ [r 2 , r] and 



< h($tz, t), since ( 3. 16 1 then gives < h(Htz, t) < CA 



-i 
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3.3 Rayleigh Taylor in the Complex Domain 



r 2 ,rl 



Let cti(x, t), crj, h(x, t) H(x, t) be as in the sections (3.1 1 and (3.2 1. 
In the time interval [r 2 , r] we define the functions 

a1(x,t) = al{x + ih{x,t),t) forxeT, t€ 

a°_(x,t) = a1(x-ih(x,t),t) for x € T, t € [r 2 , r] 

(t + (.t,<) = <Ti(x + i/i(x, i), i) for x G T, i G [r 2 , r] 

<j_(x,i) = cti(x — ih(x,t),t) for x G T, t G [t 2 , t] 

When < G [— r 2 , r 2 ] we use instead the following definitions: 



a1(x + ih(x,t),t) forxGT, t e [- 



a°_(x,t) = cr?(x - ih(x,t),t) for x G T, t G [-r 2 ,r 2 ] 
o- + (x,t) = ci(x + i?i(x,t),t) for x G T, t G [-t 2 ,t 2 ] 
a-(x,t) = (7i(x - ih(x,t),t) for x G T, t G [-t 2 ,t 2 ] 



(3.49) 
(3.50) 
(3.51) 
(3.52) 



(3.53) 
(3.54) 
(3.55) 
(3.56) 



Note that at t = r 2 , we are giving two conflicting definitions of the Rayleigh Taylor functions. 
That is a defect in th e nota tion, but it will lead to no trouble. 

From (|3.16| and (|3.37|i, we see that < h < CA^ 1 and < h < CA^ 1 whenever h, H are 

(3.57) 



defined. Hence, property 1 in section ([XT}, tells us that 

\<J±(x,t) - <r°(x,t)\ <Ch(x,t) forxGT, te [t 2 ,t]; 



and 



\<j\{x,t) - <rl{x,t)\ < Ch(x,t) forx G T, t G [-r 2 ,r 2 ]; 



(3.58) 



moreover, (3.15 i (modified: use h when t G [r , r], use h when t G [— r , r ]) together with 
( f3~23| ) and ( 346 1, shows that 



(cr± - 4) (x,t)\ < Ct w forx G T, t G [-r 2 ,r],0 < j < 2. 



(3.59) 



Again using ( 3.23[ ) and (3.46 1, and recalling the property 1 in the section ( |3.1| >, we see that 

\dia°±{x,t)\<C forx G T,t G [-r 2 ,r], < j < 2. (3.60) 
From ( |3.59| l and ( |3.60| > we obtain 

|^o-±(x,t)| < C forxeT, te [-r 2 ,r], 0< j < 2. (3.61) 

Next, we study cr±(x, t) for | |x| | < 20 A^ x t^ , t e [r 2 , t\. For such (x, we argue as follows. 
From ( |3.49] >, together with the assumption 1 in section ( |3.1| ), we see that 



and 
Hence 



T° + (x,t) - <Ti(x,t)\ = |cr?(x + ih(x,t),t) - er?(x,i)| < Ch(x,t); 



\d z a°(z,t)\ z=x+ih ( x , t ) -d x a°(x,t)\ < Ch(x,t). 



|ff+(x,0| < \<ri{x,t)\ +Ch(x,t) < Ct, 
thanks to ( |3.26| > and p.29| >; and also 

|0**!W)U=x-HM*,i)| < \d x a°(x,t)\ +Ch(x,t) < CAt*. 



(3.62) 



(3.63) 
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thank to ( ggT) and ( ggg) . 

Differentiating ( |3.49| l, we find that 

0^(7° = [S^crJ («, *) l^ar+i/.fa:,*)] {I + id x tl(x,t)). 

Since ^(a;, <)| < A^t^ < 1 by |3.30| , it follows from ( |3.63| l and \3.64\ that 



(3.64) 



(3.65) 



(3.66) 



(3.67) 



d° x <T±(x,t)\ < CA(tiy " for ||x|| < 20^-^5, i e [r 2 ,r],0 <j < 2. (3.68) 
The main results of this sections are ( |3.61[ ), ( |3.68| > and the following consequences of ( |3.57] >, 

;t) > <r°(x,t) - Ch(x,t) - Ct 10 for x <= T, te [r 2 ,r]; (3.69) 



I Strata;,*) I < CAt^. 
From p.60| >, ( |3.62| >, ( |3.65) >, we obtain the estimate 

\di<rl(x,t)\<CA(ti^ 3 for ||x|| < 20A~H% , t e [r 2 , r], < j < 2. 
Similarly 

|a£cr°(af,t)| < (ts) J for ||ac|| < 20A~HK t G [r 2 ,r],0 < j < 2. 

From ( |3.59) >, ( |3.66| l, ( |3.67| i, we obtain the estimate 

1 \ 2-3 1 

' 1 "< 20A -1 i*, * G [r 2 ,r], < j < 2. 



( [338] ), ( [339l >; 



$ta±(x, 



and 



(3.70) 



3ftcr ± (a,i) > (Xi(x,t) - Ch(x,t) - Ct 10 forxeT, te [-r 2 ,r 2 ]. 

We also need to estimate < 3a±(x, t). From (3.57 i,(3.58 1, (3.59 1 and the fact that er° (2;, t) is real 
for x £ T, we have 

\%<r±(x,t)\ < Ch{x,i) + Cr w for x E T, t e [r 2 , r]; (3.71) 



and 



|9fo-±(»,t)| < C7»(af, t) + Ct 10 forxe T,te [-r 2 ,r 2 l 



(3.72) 



3.4 Sobolev Embedding in Q(t) 

Let h(x) = h(x,t) for some fixed t € [t 2 , t], or else let h(x) = H(x,t) for some fixed < e 



3.2 



[0, r 2 ]; here, i), h(x, t) are as in section 

The goal of this section is to check that a simple Sobolev-type embedding theorem holds on 
the domain 

= {(eC : |3fl < h(UC)} 

with constants independent of the parameters A, r, k used to define h(x, t) and h(x, t). 
We write 

r± = {C € C : - ±fc(«C)}. 

Note that |/i'(a;)| < C for all ieE, with C independent of A, r, k (See ( |3TT6> , ( f337] ». Let F(Q 
be continuous on O closure and holomorphic on 0. Assume that F(C) is 27r— periodic and real for 
real £. We will prove the following assertions: 
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X\. Suppose F\y + belongs to Lip(a), with < a < 1. Then F belongs to Lip(a) as a function 
on Sl closure , and the Lip(a)-norm of F is at most C times that of F\ T+ , where C depends 
only on a. 

X 2 . Suppose F belongs to the Sobolev space H i (T + ). Then F belongs to C 3 5 (fi) and 



\F\\ c s, in) <C\\F\\ HHr+) 



for a universal constant C. 



Proof of X\\ It is enough to show that F € Lip(a) on 0, and to bound its Lip(a), and to bound 
its Lip(a) norm there, assuming that F\y + has Lip(a) norm less or equal to one. 
Let z, z' G ft. By the Cauchy integral formulas we have 



F(z)-F(z') 



1 

2TTi 



no 



i 



i 



C - z C-z' 



1 

2ttI 



f(0 



1 



1 



C-z c 



We show that both the integrals on the right-hand side are dominated by 

C{a)\z-z'\ a . 

This will complete the proof of X\ . We examine only the first integral here; the second is com- 
pletely analogous, and we have ||-F| r _ ||n P (a) = I |-F|r + 1 |u P (a) = 1 since F(() is real for real 

c 

Fix points Co, Co € T+, with 3?Co = ^z and 3?C = ^z'. 

Note that \( - ('\ < Cij^C — ^C'l forC C € T+. In particular, |Co - Col < C x \z-z'\. 
We subdivide the contour T + into 



-pnear 
1 + 



and 
Since 



{C e r+ : |3f?c - < 100(d + l)\z- z'\} 

r'f = r+ \ r^ ear . 

Jr. 



by routine contour integration, we have 
1 1 



L 



no 



c-z c 



dC 



(F(C)-F(( )) 



1 



1 



:/ (F(O-F(Co)K 
/ (^(Co)-^))-^- 

,/pnear C, — Z 



C-z c-2 
rfC 



pfar 



dC 

(F(C)-F(C )) 

(no-nco)) 



dC 



C-z' 

1 



C-z C - z' 



dC 



(3.73) 



Since 



l-P(C) - ^(Co)| < C|C - Co|° 
<C'\3tC-*KCo\ 
=C*'|3*C-^| Q 

<cic-«r, 
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it follows that 



Similarly, 



(m-nco))-^- 



™ C, — z 



< C"\z-z'\ c 



(3.74) 



(3.75) 



Also, 



<C, 



since we may deform T 1 ^ 31 to a contour of length less or equal than C\z — z'\, and distance larger 
or equal than \z — z'\ from z. 
Therefore, 



(F(Co)-F(C )) 



<C\F{C )-F{Q\ 



< c"|Co - CoT < C'l^Co - R&P = - < c"\z - 

For C e T^ 1 ', we have 

|F(C)-F(Co)|- ^--^ 

c - ^ c - % 

Since also |C - Col < C|SRC - 3%Co I = C|5RC - 5Rz| < C|C - z|, it follows that 

1 1 



(3.76) 



< C|C - Co| c 



\F{Q-F{^)\ 



|C-z| 2 ' 

s tl 

<c|z-z'nc-zr- 2 . 



Consequently, 



(nO-^Co)) 



1 1 



(3.77) 



(3.78) 



Putting p774) ,...,( |3"778) into ( |3~73| l, we learn that 

F(0 1 1 



c 



c 



< C|z-z'| c 



This completes the proof of X%. 

Proof of X 2 . Suppose F|r + € -ff 4 . Then F and its derivatives up to order 3 are bounded 
on T + , hence also on T_ (since F(£) — F(£)). Therefore, by the maximum principle, F and 
its derivatives up to 3— order are bounded on f2. Also, since F\r + E H 4 , we know that F"'(£) 
belongs to Lip(|) on T + . Applying X\, we conclude that F'" belongs to Lip(i) on Q. Thus 
F e C 3 ' 5 (£l), completing the proof of X 2 - 

When we apply X 2 in later sections, we will typically take F(() = -z M (C t) — z„(C t), where 
i))/i=i, 2 an d (2u(Ci *))/*=!, 2 are two solutions of the Muskat equation. 



4 Perturbing the Muskat solution 

Here we expose the necessary a priori estimates of the difference of two solutions of the Muskat 
equation in order to obtain the theorem [43] of subsection|4.7| 
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4.1 Differentiating The Muskat Equation 



In this section, we suppose we are given a pair of functions z^(x, t) (p, = 1, 2), z\(x, t) — x and 
Z2(x, t) periodic of period 2ir in x, analytic in the complex variable x in the region 

fi(t) = {|9far| < h(Vtx,t)}, 

and smooth in (x, t) on the closure of this region. We write 

r±(i) = {x±ih(x,t) : x G T} 

for the two contours bounding the above region. 

Later, we will take h(x, t) (of this section) to be either our previous h(x, t) or our previous 
h(x,t), depending on whether we work with t G [t 2 ,t] or with t G [— t 2 ,t 2 ]. For now, we 
simply suppose that h(x, t) is a positive, smooth function of x for each fixed t. 

We suppose that our z M satisfy the Muskat equation 



sin {zi(x,t) - zi(x + y)) (d x z,j,(x,t) - dyZ^x + y,t)) 
j cosh {z%(x, t) — z 2 {x + y)) — cos {z±{x, t) — z%{x + y)) 



dy. 



(4.1) 



Note the sign in (4. 1 1 and note that there is no factor ~ in front of the integral. 

In (4.1 1, x and y G T are real. Later we will use the analyticity of z^(x, t) in x to move into 
the complex plane. 

We will assume the following COMPLEX CHORD-ARC CONDITION: 

Let x, y G Then 



(4.2) 



\cosh(z 2 (x,t) - z 2 (y,t)) - cosOiO,*) - *i(y,t)) | > c Cj4 [||3?(x - y)|| + |$(x - y) 

Let fc be a large enough positive integer constant (her e, fc = 4 is enough). 

We apply the fc— th derivative d x to both sides of (4. 1 1. The goal of this section is to see what 
results. 

We obtain 

sin (z±(x, t) — Z\[x + y)) 
y€ T I cos h (z2(x,t) - z 2 (x + y)) - cos (zi(x, t) - zi(x + y)) 



= E < k ^) / d " 



k+k'=k 



xd x {d x z ll (x,t) ~ d x z t ,(x + y,t)}dy, 



(4.3) 



where c(fc, fc) are harmless coefficients and c(0, fc) = 1. 

We fix t, y and write lu(x) = Zu(x, t) — Zn(x + y, t). The integrand in (4.3 1 is then (with 
fc + fc' = fc) : 

sin(zi(.T)) 



x \ cosh(z 2 (a;)) — cos{z\{x)) 

a* 



= £ c^J^sHM*))} 
where c(fc' J , k b ) are harmless coefficients 



1 



cosh(z 2 (x)) — cos(fi(a;)) 



Now d!f {sin(z! (x))} is a linear combination (with harmless coefficients) of terms 



sin^ (*i 0*0)1 II l d *^( x ) 
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where fcj > 1 for each j, k\ + ■ ■ ■ + fe" « = k" and sin*-™") denotes the th derivative of the 

1 — J ' 1 mil 



sine function. 
Also, 



c> fc 



1 



_ cosh(z 2 (:r)) — cos(zi(ir)) 
is a linear combination, with harmless coefficients, of terms 

m 

[cosh(z 2 (x)) - cos(zi(a;))]~ (m+1) J[ d h J {cosh(z 2 (x)) - cos(£i(a:))}, 

3=1 

with kj > 1 for each j, and with ki + ■ ■ ■ + k m = k b . 

Therefore, the integrand in ( |4.3| l is a linear combination, with harmless coefficients, of terms 

m 

x [cosh(z 2 (x)) — cos(5i(a;))] _ ^ m+1 ^ 9^' {cosh(z 2 (x)) — cos(5i(a;))}, 

with fcj, kj > 1 for each j, and with &•[ + ■■• + fc^, +&! + ••• + fc TO + k' = k. 
Furthermore, for kj > 1, 

9^{cosh(z 2 (a;))} 
is a sum with harmless coefficients of terms 

rrij 



■1=1 



with in j > 1, with kji > 1, and with Yli^i kji — kj- 
Similarly, 

is a sum with harmless coefficients of terms 

rrij 

with TOj > 1, with each kji > 1, and with Y^=i kji = kj- 

Consequently, the integrand in ( |4.3| l is a linear combination with harmless coefficients, of 
terms 



zf^{x] 



Z x (x) 



i=l 



m' m j 

x [coM^))-cos(Mx)T im+1) j[ I [coshW (z 2 (x)j\ JT [ 

m" ( m 'j 

xjT [co S K') (il(a;)) ]rj[^) (x ) 



Z (X, 



) = 1 



(4.4) 
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where all the 
and 



k' 



TO 8 

E 



k], k'jv k 'h > 1. a11 the K< ™" ^ 1 ' 



*$+EE*$i + EE*ft = *> rn' + m" = m. 

3=1 i=l j=l i=l 



(4.5) 



(4.6) 



Let us call each expression in square brackets in (4.4 1 a BASIC FACTOR. Thus, the integrand in 



( |4.3| l is a linear combination, with harmless coefficients, of products of basic factors. We assign 
to each basic factor a weight, as follows. 



r*(*0 



i (x)] has weight 1. 
(x)] has weight 1. 



[sin^™ \zi(x))] has weight 1 if m) is even, if is odd. 
[cos 1 -"^ ' (f i (x))} has weight 1 if m'j is odd, if m" is even. 

• [cosh*™^ (z2{x))\ has weight 1 if m'j is odd, if m'j is even. 

• [cosh(z2(a;) — cos(5i(x))]~^ m+1 - ) has weight — 2(m + 1). 

The significance of this notion is that if 4>(x) is a basic factor with weight w, then \<f>(x)\ = 
0(\\y\\ w ) when \ \y\ \ is small, uniformly in x. That is because z^ fe ^ = d^z fi (x, t) — d^.z fl (x + y, t) 
and because of the COMPLEX CHORD- ARC CONDITION. 

We define the TOTAL WEIGHT of a product of basic factors to be the sum of the weights of 
those basic factors. 



We now check that (4.5 i and (|4.6|l implies that the sum of the weight of the basic factors in 



(4.4 1 is greater than or equal to zero. 

To see this, we note that the total weight of 



3 , 

^\z 2 {x))\{(z k ^{x)) 



1=1 



is 1 



to', odd + m 'j > 2, since m'j > 1 by (4.5 1. 
Similarly, the total weight of 



i=l 



is l m " odd + TO j > 2, since m" > 1 by (4.6 1. 
Also, the total weight of 



Bin< ro '>($i(a;))J [] 



(I) 



is l m » even + m" > 1, since m* > 0. 



Consequently, the total weight of the terms (4.4 1 is greater than or equal to 

1 + 1 - 2(m + 1) + 2m + 2m" = 0, 



since m + m = m by (4.6 1. 

This completes the verification that the total weight of ( |4.4| > is bigger than or equal to zero. 
We classify the terms (|4.4[) as follows. 
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• Suppose [£p (x)] or [z^ (x)] appears in the product (4.4 1. Then we say that the term 
\AA\ is dangerous. 

• If [zp' (x)] or [i® (x)] appears in a term (4.4 1, then we say that the term (4.4 1 is safe. 

• If the only [z[ (x)], appearing in a term (4.4 1 have k < k — 1, then we call term 
(|4.4|i easy. 



Thus, any term (4.4 1 (for which ( |4.5| l and ( |4.6| l holds) is dangerous, safe, or easy. These conditions 
are mutually exclusive thanks to (|4.5|l and (|4.6[) and because we will take fe > 4. 



Let us identify the dangerous and the safe terms of the form (4.4 1. Since 



3=1 3=1 «=1 3=1 «=1 



we see that the only possible dangerous term ( |4.4| i is the term arising from fc' = fc, all fc", fc^, 

= 0, m' = 0, 



( , - 0. Since all the k' Tl , k"^ > 1 by (4.5 i and (j4.6|l, we must have k' = k, 
m" = 0. That is, the only dangerous term of the form (4.4 1 is 



^ +1 )(x)l [sinOMx))] 



1 



cosh(z 2 (x)) — cos(£i(x)) 



This term appears in the integrand of ( |4.3[ l with a coefficient 1. 
Next, we look for all the safe terms ( |4.4"| ). 
We distinguish several cases. 



Case 1: The basic factor [z[ k (x)] may have k' = k — 1. In this case (4.5 i and ( |4.6| > imply 



that there is exactly one nonzero number among the fc', k'^, k'-{, and that number is 1. Since all 
the kj, fcj-j, k'-^i are at least 1 by (4.5 1 and ( |4.6) , there remain only the following possibilities: 

• = 1, k\ = 1, m' = 0, m" = 

• m tt = 0, m' = 1, m' x = 1, k' n = 1, m" = 

• m tt = 0, m! = 0, m" = 1, m!{ = 1, fc^ = 1. 

These possibilities give rise to the following term of the form (|4.4|>: 



£<*>(*)] [cos(^(x))] 
sin(£i(x)) 



1 



cosh(£ 2 (^)) ^ cos(zi(x)) 

1 



cosh(z 2 (x)) — cos(zi(x)) 



sinh(z 2 (x))][£ 2 (x)] 



z®(x)l [ S in(zr(x))] 



cosh(z 2 (x)) — cos(Si(x)) 



[sin(zi(x))][£i(x)] 



Case 2: A basic factor [z\ j (x)] in (4.4 1 has fej = fc. Then (4.5 i and ( |4~6*1 ) show that fc' = 0, 
= 1, m' = 0, m" = 0, m = 0. This gives rise to a term (4.4 1 of the form 



[z'^x)] [cos(£i(x))] z[-\x) 



1 



cosh(z 2 (x)) — cos(ii(x)) 
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(k' ) 

Case 3: A basic factor [z 2 3 (x)] in (4.4 1 has k[ }i = k. Then (4.5 1 and |46]) show that k 1 = 



rrv — 0, m! = 1, m[ = 1, k' n = k, m" = 0, and m = 1. This gives rise to a term (4.4 1 of the 
form 



[z'^x)] [sin(zi(a;))] 



1 



cosh(z 2 (a;)) — cos(£i(x)) 



[smh(z 2 (x))][z^(x)}. 



(k" ) 

Case 4: A basic factor [z\ 3 (x)} in (4.4 1 has fc" = k. Then (4.5 1 and d46l) show that k' = 



m* — 0, m' = 0, m" = 1, m![ — 1, fc^i = k, m = 1. This give rise to a term (4.4 1 of the form 



%{x)] [wa{z x {x))] 



1 



cosh(i 2 (iz:)) — cos(zi(a:)) 



[Mh{x))][^\x)]. 



These are all the safe terms. 

We have now identified all the dangerous and the safe terms. All the remaining terms are 
easy. 



We now substitute into ( |4.3[ ) the results of the above discussion of the integrand of (4.3 1. 
Recalling our definition of z\{x), z\(x), and making the change of variable u = x + y in each of 
our integrals, we obtain the formula: 



d t [dtz^x,t)] 



sm(zi(x,t) — zi{u,t)) 



i ueT cosh(z 2 (x,t) - z 2 (u,t)) - cos(zi(x,t) - zi(u,t)) 
-(coeff) 



[c£ +1 z M (x,i)-^ +1 z>,i)] du 



[d xZl (x,t) d uZl (u,t)] C os{ Zl {x,t) - Zl (u,t)) [d k_ z ^ t) _ fig^fct)] du 



-{-(coeff) 



uST cosh(z 2 (a;,t) - z%{u,i)) - cos(zi{x,t) - zi(u,i)) 

[d x z 2 {x, t) — d u z 2 (u, t)] sinh(z 2 (:r, t) — z 2 {u, t)) sin(zi(x, t) — zi(u, t)) 



ttST 



(cosh(z 2 (x,t) - z 2 {u,t)) - cos{z x {x,t) - Zl (u,t))) 2 

x [d£ Zll (x,t) -8%Zp(u,tj\ du 



+(coeff) 
+(coeff) 
+(coeff) 



[d x zi(x,t) - d u zi(u,t)][sin(zi(x,t) - zi(u,t))} 2 r . 



x (cosh(z 2 (x,i) - z 2 {u,t)) - cos(zi{x,t) - zi(u,t))) 2 



[d^z^x, t) - d£z M (u, t)] du 



[d x z,(x,t) - d u z,(u,t)} co S ( Zl (x,t) - Zl {u.t)) [qHm _ d K ZlM] du 



ue j cosh(z 2 (x,t) — z 2 (u,t)) — cos(zi(x,t) — zi(u,t)) 

[dxZ^x^t) - d u Zfj,(u,t)}smh.(z2(x,t) - z 2 (u,t))sm(zi(x,t) - zi{u,t)) 
„ eT (cosh(z 2 (.T,<) - z 2 (u,t)) - cos(zi(x,t) - z 1 (u,t))) 2 

x [d^z 2 (x,t) - d^z 2 (u,t)] du 



-(coeff) 



[d xZfl (x,t) - d u z,j,(u,t)][sm(zi(x,t) - zi(u,t))} 2 
ueT (cosh(z 2 (a;,i) - z 2 (u,t)) - cos{zi{x,t) - zi(u,t))) 2 

I/max 

+ ^c£EasyTerr<(z,t), 



[d x L z 1 (x,t)-d^z 1 (u,t)] du 



(4.7) 



i/=i 



where 



Easy Term^(x, f) 
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TT [efrzfat) - d^ Zl (u,t)] TT \d$z 2 {x,t) - d$z 2 {u,t) 
'»&i=i <=i L J 

x[sm( Zl (x,t) - z 1 {u,t))} m3 [cos{z 1 {x,t) - zi(u,i))] m4 

x[sinri(z 2 (x,t) - z 2 (u, t))] m * [cosh(z 2 (x, t) - z 2 {u,t))] me 

x (cosh(z 2 (x, t) — Z 2 (u, t)) — cos(zi(x, t) — z\{u, t)))"™ 7 du, (4.8) 

where the mi, fcj, fc^ may depend on //, i/, although we do not indicate that in the notation, and 
where 

all the hi and k\ are > 1, and 7717 > 1, (4.9) 

and 

mi + m 2 + ma + 7715 — 2m7 > 0; (4.10) 

and 

all ki and fc^ are < k - 1. (4.11) 

Now we recall that we assumed z fJi (x,t) are analytic in f2(i), and that (zi(x, t), z 2 (x, £)) are 
assumed to satisfy the COMPLEX CHORD-ARC CONDITION. 

The integrands (4.7 1, (4.8 1 have no singularities; they are holomorphic in (x,u) £ f2(t) xf2(t). 
Hence, we may deform the contour in (4.7 1, ( |4.8| l from T to T± (i), and the formulas holds, not 
merely for x real, but for all x in the closure of f2(i). 

We record this result as a lemma. 

Lemma 4.1 Let k > 4. Le? /i(x, i) fee a positive smooth function of x £ T /or each fixed t. Set 
fi(t) = {i£ C/2ttZ : |3?x| < /i(3fcc,t)}, T± = {x±i/i(x,t) :ieT}. Let zxfat), z 2 {x,t) 
be smooth functions of (x, t) on {x G Q,{t) closl "' e , t 6 Time interval}. 

For each fixed t, suppose z±{x, t), z 2 (x, t) are analytic in f2(t). 

Assume f/iflf z\{x, t), z 2 {x, t) satisfy the Muskat equation 



{ME) : dtzp{x,t) 



ye 



sm(zi{x,t) - Zijx + y)) (d x z fl (x,t) - d x z^{x + y,t)) 
1 cosri(z 2 (x,i) - 22(^ + 2/)) - cos(zi(af,t) - Zi(x + y)) 



dy, 



for (jl= 1,2. 

AssHme flfao f/ia? 21 (x, i), z 2 {x, t) satisfy the COMPLEX CHORD-ARC CONDITION 
{CCA) : |cosh(z 2 (x,t) - z 2 {y,t)) - cos (zi (x, i) - Zi{y,t)) \ 
>c CA m{x-y)\\ + \Q{x-y)\} 2 . 

for x, y € fi(t). 

TTien, /or £ G T + (t), we have that 



{I) d t [dfz^C, t)] = Dangerous f) + £ Sa/fe^C, t) + £ cJjawyJKC, *)> 

J=l v=\ 



where 

{II) Dangerous „(Cji) 



sin(zi(C,i) - z 2 (w,t)) 



iu£r+(t) 



cosh(z 2 (C,t) - z 2 (w,i)) - cos(zi(C,t) - zi(w,t)) 



x (af +1 z AI (C,<)-^+ 1 z AI ( W ,i))d W 
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(777) Safest) 

(d{Zi((,t) - d w z\ (w,t))cos(z 1 ((,t) - z 1 (w,tj) 



(coeff) 



cosh(z 2 (C,i) - z 2 (w,t)) ~ cos(^i(C,t) - z\ (w,t)) 
wer + (t) 

x (O^iCt) - dt Zfl (w,t))dw 

(IV) Safest) 

(d{Z 2 (C,t) - d w z 2 (w,t))smh(z 2 (C 1 t) - z 2 (w,t)) sm(zi(C,t) - zi(w,t)) 

\2 



= (coeff) J 



(cosh(z 2 (C,t) - z 2 (w,t)) - cos(zi(C,t) - zi(w,t))Y 

wer + (t) 

x (d^(C,t) - d^(w,t))dw 

(V) Safest) 

icocj j j i (d(Zi((,t) - d w Zi(w,t))s\v?(z 1 (C,,t) - zi(w,t)) 



(cosh(z 2 (C,i) - z 2 (w,tj) - cos(zi(C,t) - z 1 (w,t)) f 
wer + (t) 

x (d^((,t) - dt Zft (w,t))dw 

(VI) Safest) 

= (coeff) J 



_ (d(Zn((,t) - d w z^(w,t)) cos(zi(C,t) - zi(w,t)) 
( 

wer + (t) 



(cosh(z 2 ((,t) - z 2 (w,t)) - cos(zi(C,i) - z 1 (w,t))) 

x (d^ Zl (c,t) - d^ Zl (w,t))dw 

(VII) Safe* 

_ f (d c z^(C,t) - d w z fJ ,(w,t))sin(z 1 ((,t) - z^w, t)) sinh(z 2 (C, t) - z 2 (w,t)) 

J (cosh(z 2 ((,t) - z 2 (w,t)) - cos(zi(C,t) - z 1 (w,tj)) 

wer + (t) 

x (d}z 2 ((,t) - d^z 2 (w,t))dw 

(VIII) Safe%((,t) 

= (coeff) [ (dc z ^(CQ - d w z l _ l (w,t))sm 1 (z 1 (C,t) - zi(w,f)) 

J (cosh(z 2 (C,t) - z 2 (w,t)) - cos(z 1 (C,t) - z 1 (w,t))) 2 
wer + (t) 

x (d^ Zl (C,t) - d^ Zl (w,t))dw 

and each term Easy^((, t) has the form 
(IX) Easy^Ct) 

= J Uid^z^-d^z^w^^dfz^-d^z^t)) 

wer + (t) i=1 i=1 

x [sin(^(C,i) - 0i(w,i))] m3 [cos(0i(C,t) - ^iKt))]™ 4 

x [sinh(z 2 (C,t) - z 2 (w,t))} m "[cosh(z 2 (C,t) - z 2 (w,t))] m ° 

x [(cosh(z 2 (C,t) - z 2 (w,t)) - cos(zi(C,i) - Zi(w, t)))]~ m7 dw, 
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where 



and 



also 



(X) mr > 1, 1 < h < k - 1, 1 < A;- < k - 1, 
(XI) mi + m2 + ?7i3 + 7715 — Im^ > 0; 

i i 

Here, the mi, hi, k[ may depend on fi and v, though we have suppressed this dependence in the 
notation. 

4.2 Perturbing Easy and Safe terms 

Let z — (z fJi (x, t))^=i,2 and z ~ (zSx, i)) AI= i l 2 be two solutions of the Muskat equation (ME) 
as in the lemma ( |4.1) , 

We assume that z^(x, t) and z„(x, t) are defined on T x T, where the time interval T contains 
some given to 6 [r 2 , r]. 

We make the following assumptions on and z^(x, t): 

1. z fl (x,t ) and z^Xjto) extend to functions analytic on fl = < h($lx,t )} and 

smooth on the closure of Q. 

2. 

\dl%(c,t )\ 2 d$t( < c, 

Cer± 



where r± = {3£ = ±ft.(5i£, t )}. Here, h(x, t) is as in section 3.2 

3. 

| cosh (z 2 ((,t ) ~ z 2 (w,t )) - cos (2 x (C,to) -£i(w,*o)) | 
>c C a [||5R(C-^)|| + |3(C-^I] 2 , 

for (, w E il. 

4. 

IK",*) -^(-,*)llH*(n(t)) < A, 

with A as in section [3~2l . 
Here we recall that 



E/ I/(C)I 2 ^C+E/ I5c 4 /(C)I 2 ^C- 



2 

ff 4 (fi(t)) . , , 

± ^Cer±(t) ± Jcer±(t) 



Note that our unperturbed solution z is assumed to have 4+10 controlled derivatives in 
hypothesis 2, whereas the perturbation z — z is assumed to have only 4 controlled derivatives in 
the hypothesis 4. In this section, we write c, C, C", etc. to denote constants determined by the 
constants C and cqa m hypotheses 2 and 3. These symbols may denote different constants in 
different occurrences. 

From the previous hypotheses, we conclude that 

/ \d 4 ( z^C,t )\ 2 dnC<C, (4.12) 
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< C for < j < 3, C G ^ 



closure 



| cosh (z 2 (C,*o) - z 2 (w,*o)) - cos (^(Cto) -Zi(w,*o)) | 



>c C a [||»(C-«;)|| + |3(C-w) 



for £,w € £1 



closure 



(4.13) 



(4.14) 



Let Safe^iC,, to) (j = 1, 6)) and Easy£(C, to) be as in the lemma (4.1 1, arising from the 
Muskat solution (zi(C, t), Z2(C, t)), taking t = to in that lemma. Similarly, let Safe^((, to) and 
Easy^((,to) arise from the Muskat solution (z^^t), z 2 ((,t)). Our goal in this section is to 

estimate Easy^{Q, t a ) - Easy^(C, t ) and Safe^((, to) ~ Sa/e M (C, t ) . We begin with the 
Easy terms. 

From (1)...(5), we obtain the estimates 



<c*a o<i<3, cer ±; 



(4.15) 



and 



K 3 (z M - ^)(C,*b) - - ^)K*o)| < C||R(C - «;)||M [a c 4 (z M - z M )] (C) (4.16) 

for r+, where M is a Hardy-Littlewood maximal function. 

We prepare to estimate 

|[cosh(z 2 (C,< ) ~ z 2 (w,t )) - cos(zi(£fo) - zx(w,t ))] 

- [cosh(z 2 (C,t ) - z 2 (w,t Q )) -cos(z 1 (C,<o) -^(iMo))]! 

for C,w e r+. 

To do so, let £, 77, £, 7? be complex numbers, with |£|, \r)\, |£|, \rj\ < C. We set £ s = 
s£ + (1 — s)£ and ?7 S = srj + (1 — s)?7 for < s < 1. We then have 



ds 



[cosh(£ s ) - cos(?7 s )] 



= |sinh(£ s )(£ - f) + sin(?7 s )(77 - r))\ 



< C\t.\\t -CI + C'l^llr? - 77I < C"(|C| + ICDIC -CI + C(M + \r,\)\ V r,\. 



We integrate in s. 
We take 

£ = Z 2 ((,t ) ~ Z 2 (w,t ), 

1 = z 2 (C,t ) - z 2 (w,t Q ), 
V = zi(C>*o) - «i(Wj*o)j 
Q = £i(C,*o) -2i(«>,*o) 

for |C — ttjj < c, where c is a small constant and thus obtain: 

|[cosh(z 2 (C,io) - z 2 (w,t )) - cos(z 1 (C,t) - z 1 (w))} 

- [cosh(z 2 (C,t ) ~ z 2 {w,t )) -cos^Ct) -2iM)]| 
< C"[|z 2 (C^o)-z 2 (u',to)| + |z 2 (C,io)-^2( w ' < o)|]x|[2;2(C^o) - Z2(iMo)] - fe 2 (C,*o) - l 2 (w,to)]| 
+C"[|2i(C,to)-«i(«',*o)|+Ui(C)*o)-li(«',*o)|]x|[«i(C,*o) ~ zi(w,t a )] - [li(C,*o) - *o)]| 

for (,we r + , \( - w\ < c. 



(4.17) 
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Since 



|zi(C*o) -Zl{w,t )\, |^i(C,*o) -2l(w,i )|, |*2(C>*0) -^2(w?,to)|, U 2 (C,*o) ~ Z 2 (w,t ) \ 

<C\\Vt((-w)\\, 

for (, w E r + ; and since 

I t-^i (C> *o) — ^1 (w, t )] — fei(C> *o) — ^1 *o)] 1 = |[2i(C>*o)-2i(C>*o)]-[2i(w,*o)-Zi(w,to)]| 
< max|9 fi (^i -2i)(^,*o)| • |R(C~ «0| < CA|«(C-w)| 



for (,w e r+, |C — w| < c, thanks to (4.15 1; and since, similarly, 

I MC) *o)-^2(w, to)]" feCC, i o)-Z 2 (w, t )] | = |[z 2 (C, to)-I 2 (C) *o)] - i z 2(w, t )-Z 2 (w, t ) 

< max |a £ (z4 - z 2 )(£,t Q )\ • |K(C - w)| < CA|K(C - w)| 



for £, w G r + , |C — w| < c, we conclude from (4.17 i that 



|[cosh(z 2 (C,io) ~ z 2 (w,t )) - cos(zi(C,t ) - zi(w,t ))] 
- [cosh(z 2 (C,t ) - z 2 (w,t )) -cos(^ 1 (C,to) -li(^,*o)X 

2 



< CA 



sin 



( — w 



(4.18) 



for £, id € r+. We have proven this for |£ — w\ < c. For |£ — w\ > c, £, w 6 T + the 
estimate ( |4.18| l is trivial from ( |4.15| l. Together with the chord-arc condition ( |4.14| i, this yields the 
following estimate 

[cosh (z 2 (C, h) - z 2 (w,t )) - cos (zi (C, t ) ~ z x (w, to))] -1 



[cosh(z 2 (C,t ) ~ z 2 (w,t )) - cos (^(Cio) -&(w,i ))] 1 

-2 



< CA 



£ — it; 



forC, to e r+. 

Next, from ( 4.15| > we obtain easily the estimates for (, «i e L 



|cos(>i(C,*o) - zi(w,t )) - cos (^(C, to) -2i(w,i ))| < CA 
|cosh(z 2 (C,t ) - z 2 (w,to)) - cosh(z 2 (C,t ) - z 2 (iu,to))l < CA. 

Next, for T + , we have 

|sinh(z 2 (C,to) - z 2 (w,t )) - sinh (z 2 (C*o) - z 2 (w,t ))\ 

<C\[Z 2 - Z 2 ] (C, to) - [Z2 ~ Z 2 ] (W, t ) I 

<C'A||K(C-«;)|| 

thanks to ( |4.15| l. 

Because of the periodicity we conclude that 

|sinh(z 2 (C,t ) - z 2 (w,t )) - sinh (z 2 ((,t ) - z 2 (w,t ))\ 



(4.19) 



(4.20) 
(4.21) 
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< cx 



£ — w 



(4.22) 



for (,w e r+. 
Similarly, 



I sin (zi (C> *o) - (w, t )) - sin (z x ((, t ) - z^w, t )) | 



< CA 



Hill 



(4.23) 



for £, u; G T + . 

From ( 4. 15] >, we see that (for (, to £ T + ) we have 



< CA 



(4.24) 



when 1 < k <2. 

From ( 4. 16[ l, we have (for (,to € T + ) that 



j to) - fljU M («;,to)] - [^(C.to) - ^(w,*o)] 



< C 



C — id 



M[0 C 4 (^-^)] (C,t ). 



(4.25) 



We now recall the form (IX) of the term Easy%((, t) in lemma (4.1 (with k = 4. Changing the 



factors in (IX) one at a time to their analogues for z, and applying estimates (4.19 (...(4.25 I, we 
learn that 

\EasyZ({,t ) - Easy>*((,t ) 



< 



< CA + ^ CAf - ^0, to))] (C), 



(4.26) 



for £ g r + ; here we use also (XI) in lemma (4.1 1 to see that we do not pick up negative powers 
of|sin(£=^l 



It follows from 4, ( |4.26| l and the Hardy-Littlewood maximal theorem that 



Cer+ 



Easy2(t,t ) - Easy»((,t ) 



(4.27) 



Next, we estimate 



Saf^(C,t )-Safe^C,t ) 



for j = 1, .... 6; see lemma (4.1 1 for the definitions of these quantities. 

For instance let us examine Safety ((, to)- Recalling (IV) from lemma (4.1 1, we have 



Safe%((, to) = c / K{(, w) [<9 4 ^(C, to) - d 4 z^w, to)] dw, (4.28) 
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for £ e T + , where 

K((,w) = 

[9 c z 2 (C,io) - d w z 2 (w,t )} ■ [sinh(z 2 (C,io) - 2 2 (w,to))| • [sin (;?i(£, i ) - zi(w,t Q ))} 
[cosh (z 2 (C,to) - z 2 (w,t )) - cos(zi(C,t ) - 2i(w,t ))] _2 

= [Factorl] ■ [Factor2] ■ [FactorS] ■ [Factor A]. (4.29) 

Similarly, starting from the Muskat solution z in place of z, we obtain 

SafetiUo) = c I K((,w)[d 4 c z fl ((,t )-diz^w,t )]dw, (4.30) 
Jwer + 

for £ E r+, where 

|T(C,iu) = [Zactorl] ■ [£actor2] • [£actor3] • [ Factor A }. (4.31) 

We write 

= c K(C,w)[d^(z li -z li )(C,t )-d^(z li -z tl )(w,t )]dw 
Jwer + 

+c [ [K((,w)-K((,w)} [dlz^,t )-diz^wM]dw 
Jwer + 

= Terml((,t ) + Term2((,t ), (4.32) 



forC G r+. 



It is convenient that the 4— th derivatives of z„ enter into Term2, rather than the 4— th deriva- 
tives of Zu. 



From our estimates (4.13 i and (4.14i, we see that, for |£ — w\ < c, 

Factorl = A x (C, to) • (C - w) 
Factorl = A 2 ((, w) ■ (£ - to) 
Factor?, = A 3 ((, w) ■ {( - w) 
FactorA = A 4 (C, w) ■ (C - to) -4 

(4.33) 

with the derivatives of Ax,...,A^ up to order 1 less or equal than C in absolute value. 

For ||3?(C — to) || > c, Factorl, Factor2, Factor3 and FactorA we have derivatives up to 
order 1 bounded by C. 

Therefore, 

K({,w)=A((,w)-({-w)- 1 

for |C-to| < c, with the C 1 - norm of A(C,to) on T + x T + bounded by C. On ||3?(C-ty)|| > c 
the C 1 — norm of the kernel K(£, w) is at most C. Recalling how K{C,, w) enters into Terml in 
(4.32 1, we see that Terml((, to) is a singular integral operator applied to 

- z M )(-,i ). 

Consequently, 

N \Terml((,t )\ 2 d^A <CX. (4.34) 
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To estimate Term2 in ( |4.32| i, we first estimate 

\K(£,w) - K((,w)\ 
= | [Factor!] ■ ■ ■ [Factor A] - [ Factorl ] ■ ■ ■ [ FactorA ] 



J^J [Factor i] ■ [F actor j — F actor j] ■ Y\[ Factori ] 

J=l i<3 i>3 



< J^J \Factori\ ■ [Factor j — Factor j\ ■ [ F actor i [. 



(4.35) 



Our estimates (4.13 i, (4. 15 i, ( |4.24[ ) and the definition of [Factorl] in (4.29 1 (and its analogue for 
[Factorl]) show that 



\Factorl\, \ Factorl 
[Factorl — Factorl 



Similarly, ( |4.13| l, ( |4.22| i yield the estimates 

[Factor2[, \Factor2 



\Factor2 - Factor2 
Also, ( |4~T3) l and $~23\ tell us that 

|Factor3|, I Factor 3 



| Factor 3 — FactorS 
From 3, ( |4.14| ) and ( |4.19| l, we have 



< CA|C-w|- 



< C\C-w\, 

< C\\(-w\. 



< C\C-w\, 

< CA|C-w|- 



[Factor^ [ FactorA ] < C\( - w[~ i , 
[Factorl - FactorA [ < CA|C-w|~ 4 , 



when |C- w\ < 1CT 2 . 

When||R(C-tu)|| > c, 3, gig, <[4~19> yield 

|Factor4|, | Factor 4 1 < C, 
[FactorA- Factor 4 | < CA, 

Putting ( |4.36] >,...,( [4~40] > into ( |4.35| l, we see that 

|js-(c,t«) - jc(c,w)| ^cak-^i- 1 , 

if |C- to| < 1CT 2 . Also, 

\K(C,w)-K(£,w)\<C\, 
if \( — w\ > c. On the other hand, 2 gives 

|% M (C,t )-feK*o)| <c|c-H- 

Consequently, 



(4.36) 



(4.37) 



(4.38) 



(4.39) 



(4.40) 



|tf(C,«0 -=£(C>«0I • fe(C,*o) - diz^(w,t )\ < cx 
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for (,we r + . 

Hence, by definition of Term2(£, to) in (4.32 1, we have 



\Term2{(,t )\ < C\ 



for C G T+. 



Together with (4.32 1 and (4.34 1, this yields the estimate 



Cer + 



Safe%((,t ) - Safe^Ch) 



The other Safe^((, to) — Safe^((, to) may be handled similarly. Thus, we have the estimate 



go/<(C,to)-ga/g?(C,to) 

'cer+ ■ 

4.3 Auxiliary Functions 

We assume we are in the setting of section ( |4.2| i. We define 

a((,w,t) = 

sin (^(C,*) ~ zi(w,t)) 



d3?C < CX 2 



(4.41) 



delist) 



cosh (z 2 (C, t) - z 2 {w, t)) - cos (z^Ct) - z^t)) [^(COP + I^Ci)] 220 *v 2 

(4.42) 

for C, we n closme (t), where 



T+(t) = {(eC : C = s + *), x G T} 



and 



0(f) ={(eC : |9fC| < ft(RC,t), a; G T}. 

Here i G [r 2 , r]. If i £ [— t 2 , t 2 ] we replace /i(a;, £) by 7i.(a:, i). 
Note that we refer to the unperturbed solution z in defining a. 
We recall the following properties of the unpertubed solutions z (see section [3~T| : 

Zx(C, t), z 2 (C, t) are C 3 -function on n(t) closme , t G [-t 2 , t] and analytic in C G for fixed t. 

(4.43) 

kittM \z 2 (C,t)\ <C for |9fC| < c , te [-r 2 ,r]. (4.44) 
j?! (C, 0), 2 2 (C, 0) are odd functions of £ and real for real (. (4.45) 
We also recall the COMPLEX CHORD-ARC CONDITION 

| cosh (z 2 ((,t) - z 2 (w,t)) - cos (z x (C,£) - ZiO,i)) I 



> CCM [H^C-^H + I^C-io) 



forC, w G Q(i), £ G [-t 2 ,t]. 
We define 



o(C,f)= / a(Cj w,t)dw. 



(4.46) 



(4.47) 
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We notice that since a(£, w, t) is analytic in w £ 0,{t) we can deform the contour of integration 
in ( |4.47| i to obtain 



a(c,t) 



a(C, u>, t)<iu>. 



Let us assume that 



a(C,t) € c 2 (n(ty 



with ||a||c- 2 co(t)-'— ) < C, 



(4.48) 



(4.49) 



where C is a controlled constant. We will prove this fact in the next section (see ( |4.64[ )). 

For £ = 0, t = 0, we have 2„(0, 0) = 0; hence a(0, u>, 0) is an odd function of w, thanks to 
(4.45 1 again. Integrating over w £ [— 7r, 7r] and comparing with (4.48 1, we see that 

2(0,0) = 0. 



From ( |4.50| >, ( |4.49| l we obtain the estimate 

|5(C,t)|<C|C| + C|t| 

for C £ rj(i) closure , t £ [-r 2 ,r]. 

We take ( = x + ih(x, t) with t £ [r 2 , t]. 
Then ( |4.51| ) gives 



(4.50) 



(4.51) 



|5(cc + z/i(xc,t),t)| < C|a;| + C?^(x,t) + Ct 



(4.52) 



From section (13. 2b follows 



for a; e [-7T, 7r], t G [0, r]. Hence, (4.52 1 yields 

\d x h(x,t)\ ■ \a(x + ih(x,t),t)\ < CA^x 2 + CA- l \x\t + CA^bJ^x, t) 



Now 



< CA- X x 2 + CA~H 2 + CA- l h{x,t). 
A~ 1 x 2 < Ch(x,t) 



(4.53) 



for X £ [— 7T, 7T ], t £ [r 2 , t\. Therefore, (|4.53) yields 



|9 x /i,(a;,t)2(a; + i/i(a;,t),t)| < C/i(x, t) + C4" 



(4.54) 



for x £ T, i £ [t 2 ,t]. (Compare with ( |3.2l| ), ( |3.22| ». 

Also, since we assume z u (C, £) are real for real C, we see from (4.42 1 that a(£, w, t) is real 
for real £ and real w. Hence, (4.48 1 shows that 5(C, t) is real for real Q. logether with (4.49 1, this 
yields 



\Qa(x + ih(x,t),t)\ < Ch(x,t) 



(4.55) 



forx £ T, t £ [t 2 ,t]. 

Let us derive an analogues estimate for the function H(x,t) in section (3.2 1. The proof of 
< |4.55[ > yields also 

\<Sa(x + ih(x,t),t)\<Ch{x,t) (4.56) 

for a; £ T, t £ [-t 2 ,t 2 ]. 
Also ( |4.51[ ) gives 



|a(x + ih(x, t),t)\< C\x\ + Ch(x, t) + C\t\ 
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for x G [-7T, 7r], t e [-T 2 , r 2 ]. 
From ( |3.37) >, we see that 

for the same range of x, t. 
Hence, 



\d x h(x,t)\ <CA^\x\ 



\d x ti(x, t)a(x + ih(x, t), t)\ < CA^x 2 + CA- 1 \x\\t\ + Ch{x, t) 
< CA^x 2 + CA-H 2 + Ch{x, t) 

for x G [-7T, 7r], t G [-r 2 , r 2 ]. 
Since 

CA- X x 2 + CA-H 2 < CA- 1 ! 2 + CA^t 4 < Ch(x,t) 
for x G [-7T, 7r], * G [-t 2 , t 2 ] (thanks to (|3.37}), it follows that 



\d x h(x, t)a(x + ih(x, i),t)\ < CK(x, t) 



(4.57) 



for a; G T, t G [-T 2 ,r 2 ] 



Our basic results on the size of a and its imaginary part are (4.54 1 (4.57 1 



4.4 Perturbing Auxiliary Functions 



Let z, z be two solutions of Muskat, that satisfy the assumptions in the section (4.2 1. We also 
assume that z satisfies the additional assumptions made in the section \A3\. 



We define a((, w, t) and 5(C, t) as in \AA2) and ( |4.47 1. 

For all t for which the perturbed Muskat solution z(-, t) is defined, we let a(£, w, t) and a(£, *) 
be the analogous expressions for z, i.e., 



sin(zi(C,t) - zx(w,t)) 



cot 



£ — w 
2 

(4.58) 



cosh(z 2 (C,t) - z 2 {w,t)) - cos(«i(C,t) - «i («;,«)) [c^i(C,*)] 2 + [9 c z 2 (C,t)] 2 
for 

|9fCl < l$H < M^*M) if* G [t 2 ,t] 

|9fCl < fi(RC,*), 1^1 < fi(Mu>,t) if* G [-t 2 ,t 2 ]. 
We will show that, since z(C, i) G C 3 (n(t) closure ) the function 

o(C, w, t) is C* 1 in w for (,io£ f7(t) c/osure and has C^-norm at most C. (4.59) 
To prove ( |4.59| l, we may restrict attention to 

= {C, w G O(i) dosure , * G [0, r], |C - w\ < c}, 

for small c. Outside this region there is no problem, thanks to the complex chord-arc condition 

Then 

. / , > x / xx sin (zi(C, t) — zi(w, t)) , . „ . , 

sin ^c* - z 1 (w,t)) = v ; — ^ c,t - *i K* 
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= G 1 ((,w,t)-(C-w), 

with the C 2 — norm of Gi(C w, at most C for fixed t. Examining the asymptotics as w — > (, 
yields 

G 1 (t,(,t) = d c z 1 (t,t), 

hence 

Bin(*i(C,t)-ziKt)) = [a c z 1 (C,t)+Ga(C,v> > t)(C-w)]{C-n>), (4-60) 
with the C 1 — norm of ^(C, w, t) at most C for fixed t. Similarly, 

i / / > i\ / ,« -, cosh (z 2 (CO - z 2 (w,t)) - 1 2 

cosh (2:2 (CO -^KO) - 1 = ——— ; ^2 0*2 (CO - z 2 {w,t)) 

0*2 (CO - 22(w,i)J 

= G 3 (C,tM)(C-«0 a , 

with the G 2 — norm of G 3 (C to, t) at most C for fixed t. Examining the asymptotics as w —> (, 
gives 



g 3 (cco = ^(^ 2 (co) 2 - 



hence 



cosh (z 2 (C - z 2 (w, t)) - 1 



^(d c z 2 ((,t)f + G 4 ((,w,t)((-w) 



with the C 1 — norm of C?4(C w, t) at most C for fixed i. Similarly, 

"1 



1 — cos (21 (CO — z \{w,t)) = 



(9 cZl (C0) +G 5 (C W ,0(C-^) 



with the G 1 — norm of G 5 (C w, t) at most G for fixed t. Adding, and recalling the COMPLEX 
CHORD-ARC CONDITION ggg) , we see that 



cosh (z 2 (C0 _ ^2(10,0) _ cos(zi(C0 ~ z i( w ;0) 
(9 c z 2 (C 0) 2 + (^i(C. 0) 2 1 ■ [1 + G 6 (C w, t) ■ (C - w)} (C - wf 



(4.61) 



with the G 1 — norm of Ge{C, w, t) at most G for fixed t. Dividing ( 4.60 1 by (4.61 1, and noting that 

1 . G 6 



1 



we obtain 



1 + G 6 • (C - w) l + G 6 -(<;-w) 

sin (2:1 (CO - zi(w,t)) 

cosh (z 2 (C0 ~ z 2(w,t)) - cos (zi(C _ ^1(^)0) 



2fl fi(C,*) + Gy(c t) . (c _ w) 

(d c z 2 (C0) 2 + (^i(C0) 2 



1 



29 c zi(C0 



1 



(d c z 2 (C0) 2 + (^i(C0) 2 (C-«0 

with the G 1 — norm of G?(£, w, t) at most G for fixed 



(C-u;) 

+ G 7 (C,w,0, 
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Since also 



2d c z 1 ((,t) cot (t~ w 



(d ( z 2 (c,t)) 2 + (d <Zl (ct)) 2 



2 



+ G 8 {C,w,t), 



(d c z 2 (c,t)) 2 + (d (Zl (Ct)) 2 C - w 

with the C 1 norm of G$(£,w,t) on £2 bounded by C, we conclude that a((,w,t) belongs to 
C 1 (£2), with norm at most C. 

Again recalling that there is no problem outside £2 (i.e., when |£ — w\ > c even after (, w are 
translated by multiples of 2ir), we conclude ( |4.59[ ). 

Next, we notice that 



S(C,t) = / a((,w,t)dw = / a(£,w,t)dw (4.62) 



ift€ [r 2 ,r]; 

a(£, t) = / a((,w,t)dw = / a((,w,t)dw (4.63) 

J$iw=h(3tw,t) JwET 

if te[-r 2 ,T 2 }. 

We shall study the regularity of the function a(£) for £ G r + and 

T + = {CeC: ( = x + ih(x), x G T}. 

We recall that /i(a;) is a smooth function and we will prove that a has C 2 -norm bounded by a 
controlled constant. (We replace h by h if we are in the interval [— t 2 , t 2 ].) 

Remark 4.2 JVof/ce fto sznce z((, t) G C 3 (f2(t) c/o ™ re ), /ram ( |4.62| i ana! ( |4.63| l we see that a is 
a C 2 -function. However, this argument is not enough to control the norm 

l|fi(- + *M-.t))llc»cr) 

by a controlled constant. We would obtain a constant which depends on — ^7 — tt and we can 

not permit it. 

In fact, we will prove that a has bounded norm in C 2+s (Q(t)) for any S < 1/2. To see this, 
it is enough to show that 

a|r+(t) € C 2+s with norm at most C(5), for any 6 < -. (4.64) 



(The reduction from Q(t) to T + (t) is possible, thanks to Xi in section 3.4 1. 
To prove ( |4.64[ ), we first note that 

ntr +\ pi/ / sm( Zl ((,t)- Zl (w,t)) 

a{Q,t)=P.V. — — — — — -T7 — — — — —dQ (4.65) 

Jv+{t) cosh(z 2 (C,t) - z 2 (w,t)) - cos(zi(C,i) - Zi(w,t)) 



for £ G T + (t), since 

P.V. 

'r+(t) 

as we noted in section |2] 



P.V. I cot ( ^— - |dffl = 0, 
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It is useful to make the change of variable 



and to set 



for jit = 1, 2. 
Then 



£ =x + ih(x, t), 
w =x + y + ih(x + y,t), 

z^(x) = z^(x + ih(x,t),t) 



z\(x) — x and ,§2(2;) belong to C 3 ' 5 (T), with norm less or equal than C. 



(4.66) 



(4.67) 



Hence also 



z^(x) - z^(x + y) 



z'{x + sy)ds 



(4.68) 



belongs to C 2 5 (T x [— x, ir]) with norm less or equal than C. 
Our change of variable ( |4.66| l gives 

~, ,. h , . W s pi/ r sinjz^x)- z^x + y^jl + id^jx + y^)) 

a(x + in(x,t),t) = P.V. I — - ,„ , , —. rr . . -. , —. rrdy (4.69) 

J_„ cosh(z 2 (a;) - z 2 (x + y)) - cos(zi(a;) -z 1 {x + y)) 

for x ET. 

It is convenient to introduce the analytic functions 



ftm(C) = 



sin(C) 

c : 



Jcsh(C) = C ° Sh c ( P ~ , 9cs(0 = 



1 - cos(C) 

C^ 



for C € C and 



?L(cc_) = 

<?l(C,C) 



ffsn(C) -flsn(C) 

c-c : 



c-c 

8 (C)-g cs (0 

c-c ! 



for c, C e c. 

We can then write 



y 1 sin(li(x) - £1(2; + y)) = g sn (zi(x) - z x {x + y)) 

and 

y~ 2 [cos\i(z 2 (x) - z 2 (x + y)) - cos(zi(x) - z±(x + y))] 



z~i(x) — z\{x + y) 



(4.70) 



(4.71) 



(4.72) 



g C sh(z 2 (x) - z 2 (x + y)) 



z 2 {x) - z 2 (x + y) 

y 



+ g cs (zi(x) - zi(x + y)) 



zi(x) - h(x + y) 



y 



(4.73) 
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Thanks to p37) and ( |4~68) , the right-hand sides of p772) and gj3j belong to C 2 5 (T x [— k, %]), 
with norm at most C. Also, the left-hand side of ( |4.73| l has absolute value larger or equal than c, 
thanks to the chord-arc condition. 



Therefore, we may express the integrand in ( |4.69[ ) in the form 



Q(*,y) 



where 



n( v _ [LHS of@J2$](l + id x h(x + y,t)) 
W> v) - [LHSof(|473)l 



belongs to C 2 ' 5 (T x [— tt, n]), with norm at most C. That is, 



a(x + ih{x,t),t) = P.V. 



* Q(x,y) 



dy = 



* Q(x,y) - Q(x,0) 



dy. 



(4.74) 



with 



< C. 



IIQII 

It follows easily that the function x M- a(x + ih(x, t),t) belongs to C 2+S (T) with norm less or 
equal than C(S) for each 6 < \. This completes the proof of ( |4.64[ ). 
To conclude this section, we prove that 



sin(^i(C,i) - Z\{w,t)) 



sinfei(C)*) -&y{w,t)) 



cosh(z 2 (C t) — z 2 (w, t)) — cos(zi(C, t) — Zi(w, t)) cosh(z 2 (C, t) — z 2 (w, t)) — cos(z 1 (<^, t) — 2 1 (w, t)) 



A 



e(c,*)cot 



Q — w 



(4.75) 



for (,w e r+(t), where 6 G C°(T + (t)) and 9 B G C°(T + (t) x !+(*)) for fixed £, with norms 
less or equal than C. 

To see this, we again use the change of variable 

C =x + ih(x, t), 
w =x + y + ih(x +y,t); 

and we introduce z^{x), ZAx) for fi = 1, 2, where z^(x) is given by ( |4.66| l and $lAx) is the 
analogous expression with z^ in place of z^. 

From ( |4.71| i and ( |4.72| i and the analogue of ( |4.72[ ) for the % , we see that 



y 1 [sin(zi(x) - zi(x + y)) — sin{z t {x) - z^x + y))] 



= [g sn (zi(x) - zi(x + y)) - g sn (h( x ) ~ h( x + v))] 
+ [9sn(h( x ) - h( x + y))] 



(zi(x) - z 1 (x + y)) 



(zi(aQ - zi{x + y)) _ (z^x) - Ij^aH- y)) ' 

y y 

9L(zi( x ) + y), ii(ar) -z 1 (x + y))] pi - i x ){x) - {zi - I x )(a; + y)] 

zi{x) - h{x + y)~ 



+ [g sn (zi(x) - zi(x + y))} 



(h - h)( x ) - {h - h)( x + y) 



Since z„ 



belongs to C 3 5 with norm less or equal than CA, we know that 



(4.76) 



belongs to C 25 with norm less or equal than CA, as in the proof of ( |4.68| l. Therefore, ( |4.76| l 
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shows that 

y~ 1 [sm(z 1 (x) - zi(x + y)) - sin(i 1 (a;) - z ± (x + y))] 

belongs to C 2 5 (T x [—ir, n]), with norm less or equal CX. (4.77) 

A similar argument involving the function g\ sh and g\ s in ( |4.71| i shows that 

y~ 2 [{cosh(z 2 (a;) - z 2 (x + y)) - cos(zi(x) -zi(x + y))} - 
{cosh(z 2 (a;) - z 2 (x + y)) - cos^O) - z x (x + y))}] 

belongs to C 2 ' 5 (T x [—n, tt]), with norm less or equal than CX. (4.78) 
We have seen that 

y^ 1 sin(zi(x) — z±(x + y)) and y~ 2 {cosh(z 2 (x) — z 2 (x + y)) — cos(zx{x) — Zi(x + y))] 
belong to C 2 ' 5 (T x [—it, tt]) with norms less or equal than C, and that 

|y _2 {cosh(z 2 (a;) — z 2 (x + y)) — cos(zi(x) — Z\{x + y))}\ > c on T x [—tt, tt]. (4.79) 

From < |4.77) , ( [4.78) 1, ( |4.79| > it now follows easily that 

sin(£i(x) - zx{x + y)) 

cosh(z 2 (x) - z 2 (x + y)) - cos(ii(i) - z x (x + y)) 
sin(i 1 (a;) -z^x + y)) 



cosh(z 2 (x) — z 2 (x + y)) — cos(z 1 (x) — 'z x {x + y)) 



(4.80) 



for a function Q» e C 2 5 (T x [— 7r, tt]) with norm less or equal than C. 

Our result ( |4.80[ ) is much sharper than our assertion ( |4.75| ) In particular, we have proven 
44775] ). 

4.5 Perturbing the Dangerous Term 



Let z, z be as in the section ( |4.2| >. Suppose that z satisfies the additional assumptions made in the 
section (14. 3b. 



Let <z(£, w, t), a(C, w, t), a(£, i), a(£, i) be as in section (4.3 I. 
We recall from section ( |4. 1} that 

Dangerous ^(^t) = 



cosh {z 2 ((,t) - z 2 (w,t)) - cos (zi (C, t) - zi (w, t)) C 10 
t«er + (t) 

(4.81) 

forCer+(t). 

For the Muskat solution z in place of z, the analogous expression is 

Dangerous (Cji) — 

w ,^/ in( 7 l(c ;^- l( rfA, — t^t flTVc *) - *)]*« 

cosh (z 2 (C, t) - z 2 (w,t))- cos (MCt) - z^w^t)) t 
wer+(t) 

(4.82) 
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forCer+(t). 

Our goal in this section is to compute Dangerous t) — Dangerous (£, t), modulo small 

M 

errors. Let £ 6 be given. From (4.81 1 and (4.82 1, we have 

Dangerous ^{C,t) — Dangerous (£, f) = 



«er+(t) 



cosh (z 2 (C,i) - Z2(w>*)) - cos(zi(C,t) - zi(tM)) 

f sin (zi(C,*) - zi(w,t)) 

\cosh(z 2 (C,*) ~ z 2 (iM)) - cos(zi(C,i) - Zl («>,*)) 
™er + (t) 

sin (^(Ct) ~ z^w.t)) 

cosh (z 2 ((,t) -z 2 (w,i)) - cos (zi(C,t) - z^w.t)) 

xid^z^Ct) - di +1 z^w,t)]dw. 

Substituting the definition of a(£, w, t) and applying ( 4.75| l, we conclude that 

Dangerous ^{C,t) — Dangerous (C;i) = 



dw 



(0c*i(C,*)) 2 + (^(C,*)) a 



cot 



5 c 4+i (z„ - zjit, t) - d^ +1 (z, - z^{w, t) 



dw 



wer + (t) 

4+1, 



+ J a((,»,i)[f(^-^)(C.')-rfe-^)(»,t) 
C — w 



dw 



A|8(C,t) cot 
tuer+(t) 
Recall that we define 



+ e^C^i) [9 c 4+1 ^(C,i) - &}+%(w,t)]dvi. (4.83) 



A r+(t) F(C) 



2tt 



cot 



C — w 



[F'(C) - F'(w)]dw 



(4.84) 



for holomorphic function F. Also, in the terms involving a(£, to, i) in (4.83 I, we recall (4.62 1, 
and we integrate by parts in w. Thus, ( |4.83| l becomes 

Dangerous fj,(^,t) — Dangerous = 



-2tt 



8<*i«,t) 



(9 c z 1 (C,i)) 2 + (9 c z 2 (C,t)) 2 



[A r+(t) (^(-,t)-9^M))] 



+a(c,t)^ 1 (* M (c,t)-z M (c,t)) 

[9«,o(Cj w ^ - ^){w, t)dw 



wer + (t) 

+ X I \e((,t)cot( C —^) +0\( 1 w,t)\[d 4 c +1 z((,t)-d 4 w +1 z(w,t)}dw. (4.85) 



«er+(t) 
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Recall from assumption 2 in section ( |4.2| i that 

\d?%(C, t) - 8*%{w, t)\ < min{C|C - w\,C}. 

(Here, it is important that we refer to z, not to z.) 
Since also 

|6(C,*)|, |e»(c,«;,t)| <c 

by ( |4.75| >, we have 



e(c,*)oot (j-^) + e*(c, «/,<)} • [a c 44 %(c,*)-d c 4+1 z>,t)]^ 



for any Ce 
Also, since 



|9 w a(C, «;,t)| <C 



by (4.59 1 and assumption 4 in section [472] (recall, we assume A < 1), we have 



[d w a((, w, t)]d* (z M - t)dw 



wer + (t) 



< CA 



for all £ £ r _)_(<), thanks to assumption 4 in section (4.2 1. 
Substituting ( |4.86| l and ( |4.87| > into ( |4.85[ ), we see that 



Dangerous t) — Dangerous (£, t) 



-2tt 



(d (Zl (C,t)) 2 + (%* 2 (C,*)) 2 



[Ar + ( t )(a 4 ^(^)-a 4 z M (-,t))] 



+a(C,t)a 4+1 (2 M (c,t)-^(c,t)) 

+ Error PDT ((,t) 

forC G r+(i), with 

|Error P7 3 T (C,t)| < CA. 
The basic results of this section are ( |4.88| l, ( |4.89| l. 



< CA, 
(4.86) 



(4.87) 



(4.88) 
(4.89) 



4.6 The Main Energy Estimate 

In this subsection, we suppose that z is an unperturbed solution of Muskat problem, while z is a 



perturbed solution, as in the last several sections. Let a\, o~\, a±, h and h as in the sections 3.1 
3T2landl3~3l 

Our goal is to estimate (from below) the quantity 



dt 



Cer + (t) 



|^(^-^)(c,*)r^c. 



To do so, we first recall (I) from lemma (4.1 1. Applying that result to both our Muskat 
solutions, we see that 

*„)(£*)] = 
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[Dangerous t) — Dangerous (£, t)] + jtJSafeffc, t) — Safe^iC, t)] 



j'=l 



l<v<v max 

forCer + (t). 

According to ( |4.27| >, ( |4.41| i and ( |4.88| >, ( |4.89| l, this implies that 



-2tt 



9 c *i(C,*) 



(9 c2l (C,t)) 2 + (a c z 2 (c^)) 2 



ft$(**.-**)(c>*)] = 

[A r+( t) *) - feO, *))] +o(C, *)9 4+1 MC, *)-^(C, *)) 



forC G r+(i), where 



Recall our assumption 



and note that 



+ Error(£, t) 
[ \Error(C„i)\ 2 dK(;<C\ 2 



Cer+(t) 



|^-^)(C,*)| «C< A 2 



(4.90) 
(4.91) 

(4.92) 



1 d 

2 <ft 



Cer + (t) 



|^(«/*-^)(C,*)| <«C 



= 

Cer+(*) 
where 



if i G [r 2 , r] and 



r+(t) = {x + ift(a;) : x G T} 



r + (i) = {x + ift(a;) 



if £ G [— t 2 , t 2 ]. (In this range we use ft in place of ft, in the previous expression.) 
We conclude that 

'Cer+(t) 



1 d 

2 d£ 



= 3? 



-27T2i(C,t) 



Cer+(t) 



K(C,0) 2 + (4(C,t)) 2 

+ Error 



(4.93) 



where 
and 



F(o = a c 4 (z Al -z Al )(c,t), 



I Error I < CA 2 



(4.94) 



(4.95) 



Recall | |4.54) ,...,( |4.57| ), ( |4.49[ > and ( |4.64[ ). Our assumptions on A, A, r, k in section [pimply 



< A < r 



(4.96) 
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Recall from section|43j equations ( j4.49| l, j4.55) , j4.56\ and ((437), that 
5(C, i) is a C 2 — function of £ € T + (t), with norm at most C; 

|9fi(C,t)| < Ch{$tC,t) forCer + (t) ; 

|<9 x /i(x,t)a(C,*)| <C/!(i,t) + Ci _1 T 2 for C = x + ih(x, t) eT+(t). 



(4.97) 
(4.98) 
(4.99) 



Estimates (4.98 1, (4.99 1 hold for t e [r 2 , r]. When i e [-t 2 , t 2 ], we use fi. in place of h. 
Note that 

Q{a(x + ih(x, t),t) ■ (1 + id^/iO, t)) -1 } 

= [3a(x + ifc(x,t),i) - d x h{x,t) ■ Ra(x + ih(x,t),t)] ■ (l + (^(x, i)) 2 )" 1 . 
Hence ( |4T97| >,...,( |499l show that 

d(x + ih(x, t), t) ■ (1 + id x h(x, t))^ 1 = cjr(x) + iaj{x), 

where 



and 



or(x), aj(x) are real, C 2 — functions on T, with norm < C 

Cfe(^,t) + CA- 1 r 2 
(a?) | < — 77T-T-, — rrrr; — forallxeT. 



\ + {d x h{x,t)Y 

If t G [ — t 2 , r 2 ], then we use h in place of h. 
Let 
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/(x) = F(x + ^(x, t)) = ^ - z M )(C, t) 



£= x-\-ih(x,t) 



Thus, (4.92 1 gives 



Note that 



|/(x)| 2 dx < CA 2 



/'(x) = (1 + id x h(x, t))F'(x + ih(x, t)). 



(4.100) 
(4.101) 
(4.102) 

(4.103) 

(4.104) 
(4.105) 



We have 



Cer+(*) 



^(05(0^(0^ 



/(x) • {[a fl (x) + m/(x)] ■ (1 + icU(x, t))} {(1 + iS a ft(s, i))- 1 /'^)} dx 



Cer+(t) 

/(x) • [a fl (x) + ia I {x)]j'{x)dx 
1 d 



or(x) 



(|/(x)| 2 )dx + 5R / f(x) ai (x)if(x)dx 



1 



a'fl(x) • - (|/(x)| 2 ) dx + / _ /(x)a 7 (x)i/'(x)dx. 



(4.106) 



The first term on the right is dominated by CA in absolute value, thanks to (4.101 1, (4. 104 1. By 
(4.101 1, (4.102 1 and the sharp Garding inequality, the second term on the right is bounded below 
by 



-» / f{x) [Ch(x,t)+CA- 1 r 2 ] (l + (d x h(x,t)) 2 )- 1 Af{x)dx-C J |/(x) 

xGT i£T 
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Putting the above remarks and ( 4. 104 1 into ( |4. 106| l, we see that 



5ft / F(0~a((,t)F'(()d( 
J Cer+(t) 

>-C\ 2 -Cft[ J{x)[h(x,t) + A- 1 T 2 ](l + {d x h(x,t)) 2 y 1 Af(x)dx. (4.107) 
JxeT 

It is convenient to use (|2.6& to conclude that 



U, 7w ^ t)+A "'^ i + wkV A/( ' ) * 



'\2 



< C / |/(x)| 2 cfa < C'A 



Together with (4.107 1, this yields the estimate 



st / F(c)a(c,t)F'(C)dRC 

JCer + (t) 



7(x)[>,(x, t ) + A-V ^ 1 + -^; ) '»» A/ W ^. 



> -CA 2 - C5i 
Recalling that 

A/O) = (1 + id x h(x, t))A r+ ( t) F(x + ih(x, t)) + Error(x), 

with 

\\Error(x)\\ L , m < C\\f\\ L 2 (T) 

(see ( |2.10| i), we conclude that 



31 / F(C)a(C,t)F'(C)«C 
JCer+(t) 

> -CA -CK / /(x) [h{x,t)+A t \ rTTxui — 7^2 A r+{t) F{x+ih(x,t))dx, 

JxeT 1 + {o x n(x, tj) 

where again we use ( 4. 104| to control the term involving Error(x). Recalling the definition 



(4.103 1 of /, we can conclude that 



31 / F(C)a(C,^'(C« 
JCer+(t) 

> -CX 2 - CM ( Fj(j [h(x, t) + A- x t 2 ] A r+(t) F(C)d5RC- 
Substituting this inequality into ( 4.93| >, we see that 



1 d 

2dt 



Cer + (t) 



|^(^-^)(C,i)| 2 ^C 



>3? / F(C) ■ {[a+($lt,t) - Ch(x,t) - CA-W 2 ] A r+(t) F(C) 
Cer+(t) 

+ {id t h{®.C, *)) F\C)\d^C - CA 2 , (4.108) 

where we recall that 

° + m\t)= J^^7—r 2 force r + (t). 
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4.6.1 The main energy estimate I 



Now suppose t G [r 2 , r] . We prove a lower bound for the right-hand side of ( 4. 108 i 
We set up the partition of the unity 1 = 6^ 



8 out on T with 6 in , 9 out > 0, 

supp6 in c{\\x\\ < 
supp6 out C{||x|| > 10A-M}, 



and 



Pi 



< c 



(a- 1 **)' 



for < j < 3. 

Of course and # out depend on t, but we have fixed t above. 

Our plan is to insert a factor in or 6 out in the integrand in ( 4.108| l and study the two resulting 
integrals. 

First we study the term involving 9 out . We prepare to apply lemma 2.3 Here, we take 

h(x) = h{x,t), h t (x) = d t h(x,t), = e out , u\ = a+ - Ch - CA~ 1 T 2 md S = 6A 2 . We 



check the hypotheses of lemma 2.3 with constants of the form C(A). We write C(A) to denote 
constants depending only on A. 



Hypothesis (1) merely says that h(x, t) > for t E [r 2 , r] 
Hypothesis (2) asserts that 

\d t h(x,t)\ < 6A 2 h(x,t) 



which is immediate from ( |3.16[ ). 



on suppQ out . This assertion follows at once from ( |3.20) >. 

Hypothesis (3) asserts that the C 3 — norm of h(x, t) and d t h(x, t) as functions of x for fixed 
t, are dominated by C(A). This assertion follows at once from ( |3.23| > and ( |3.24] >. 

Hypothesis (4) asserts that the C 2 -norm of the function fii is bounded by C(A). We know 
that the C 2 -norm of Ch(x, t) + CA^t 2 (for fixed t) is bounded by CM" 1 (see ( |3~23| ». Also 
the C 2 -norm of ct+ is at most C by ( |3.61| l. Therefore, hypothesis 4 holds. 

Hypothesis (5) asserts that the C 2 — norm of the function 



x 1 — y 



0out{x){l + ih!{x 1 t)y 1 {d l {x) +d t h(x,t)) 



is at most C(A). 



Since the C — norm of h(x, t) (for fixed i) is at most CA , it is enough to show that 



x 1 — y 



%ut(x)(ai(x) + d t h(x,t)) 



(4.109) 



has C 2 — norm at most C(A). 

We estimate the derivatives of this function up to 2— order. 

For < 10A~ 1 t2 , we have Q ou t(x) = 0. Hence, the function in question vanishes in that 
interval. 

For 1 1 a; 1 1 > 20A~ 1 t 2 , we have 6 out (x)=l. 

Since we have already checked that the C 2 — norm of &i(x) is at most C, and since the 
C 2 — norm of d t h (for fixed t) is at most C(A), we see that the derivatives up to order 2 of 
the function in (4.109 1 are at most C( A) for | \x \ \ > 20 A^ 1 ^ . 

It remains to examine the region {lO^l^ 1 ^ 2 < \ \x\ \ < 20A~ 1 t 2 }. 

First we recall that 



\dla + (x)\ <CA{t^) 2 -i 



(4.110) 
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forO < j < 2, \\x\ \ < 20A~HK by ( |3~68T >, and 

\di{d t h(x,t)}\ < CA{t^ f- j (4.111) 
forO < j < 2, < 20^-^2, by (|3T34j», and also 

\di6 0Ut (x)\ <C(A-Hi)- j (4.112) 

forO < j < 2, H xji < 20^-^2. 
From ( |3.31| >, we have 

\dih(x,t)\ < C(t^) 2 - j (4.113) 
for <j <2,\\x\\< 20A- 1 t^. And also, obviously 

I^^V 2 ! < C(t5) 2 ^' (4.114) 

forO < j < 2, ||z|| < 20A"H2, sincet E [t 2 ,t]. 
From ( |4.110| i, ( |4.113| l, ( |4.114| ), we see that 

{04. (a) - Cft(a?,t) - CA^t 2 }] < CAiti) 2 ^ 
forO < j < 2, < 20A-^3, i.e., 

I finite) I < CA^) 2 "-' 

forO < j < 2, < 20A~Hi. 
Together with (4. 1 1 1 1, this yields 

|fl£{ffi(aO+fit/i(;r,t)}| < CL4(^) 2 - j 

forO < j < 2, < 2(L4" 1 ii 

In turn, this estimate and ( |4.1 12[ ) yield 

\di [e out (x){a 1 (x) + d t h(x,t)}}\ < cA 3 (t^) 2 -i 

for < j < 2, < 20A~ 1 ti , completing our verification of hypothesis 5. 
Hypothesis 6 asserts that 

M{d out (x)(l + iti(x,t))- 1 (a 1 (x,t)+d t h{x,t))} > 

on T, i.e., 

1 + 9 {h'(x\)) 2 i^( x ' t )+ d t h ^,t) + dxhfatySa^t)} > 
on T, which will follow if we can prove that 

9?CTi(x,i) +d t h(x,t) + d x h(x,t)$sai(x,t) > (4.115) 



on T. By ([3T69j>, ([3771], we have 

3?&x(i,t) > a°(x,t) - C/i(x,t) - CVTV 2 (4.116) 

for x E T, and 

|Q#i(x,*)| < C7i(M) + CVTV 2 (4.117) 
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for x eT. We recall that 



Also, 



~27rzUx) 
\d x h{x,t)\ < CA~ X 



(4.118) 



by (|X23l. 

The above remarks tell us that 



3f?CTi(a;,t) + d t h(x,t) + d x h(x,t)^sai(x,t) 
> a^(x, t) + d t h(x, t) - Ch(x, t) - CvTV 



(4.119) 



for all ieT. 

The right-hand side of (4.1 19 1 is positive for all ieT; see (3.21 1. Therefore, (4.1 15 1 holds, 
completing the verification of hypothesis 6. 

Note also that < 9 ou t < 1, hence hypothesis 7 is satisfied. 

Thus, hypotheses 1,...,7 hold here, with constants C(A). Applying lemma (2.3 1, we obtain 
the following inequality. 



J? 



F{Q9 out {^Q 



-2irz[((,t) 



Cer+(t) 



(4(C,i)) 3 + (4(C,*)) s 



+ 



Ar +(t) ^(C) 



> 



(id t h(M(,t))F'(C)}d$lt 
-i- J c 



Cer±(t) 



l^(C)l 2 «C- 



(4.120) 



Next, we prepare to apply lemma (2.6 1, taking = h(x) = h(x,t), h t — dth( x, t), 
tj\ = (T_|_ -CTi-CA" 1 ?- 2 and (5 = 20^-^3. We now check that the hypotheses of lemma \2.6\ 
hold here, with constants C(A). 

Hypothesis (1) simply asserts that h > 0. Since t E [t 2 ,t], this follows at once from ( |3.16[ ). 

Hypothesis (2) asserts that (for fixed t) the C 2 — norm of &\ and the C 3 — norms h(x, t) and 
d t h(x, t) are bounded by C(A). This is contained in our previous results ( |3.61) >, ( |3.23| ) and ( |3.24] >. 

Hypothesis 3 asserts that the C 2 — norms of 9i n (x)a±(x) and 6i n dth(x, t) are at most C(A). 

Recall that supp9 in C {||x|| < 20A~ 1 t^}, where we have the estimates 



dx r 
d^ 3 



\ 6 ln {x) <c(AtH^ 3 for0<j<2, 



dx r 
d^ 3 



dx t 
d^ 3 



dx r 

d} 3 
dx 



0in{x) 

d t h(x,t) 
<r+(x,t) 
h(x, t) 



< C{A) (A-H* 



< C(A) lA-H 



2-j 



2-3 



for < j < 2 (see fl333fr ), 



for < j < 2 (see fiMfo , 



A-\ 2 



<c(t*) 

<C(**) 



2-i 



2-J 



for < j < 2 (see COT) ), 



forO < j < 2, i g [r 2 ,r). 
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Hence, 



dx r 
dx 



6 in (x)d t h(x,t) 



<C(A)(t^ J for0<j<2, (4.121) 
<C(A)(t^Y 3 for0<j<2. (4.122) 



In particular, the C 2 — norm of 9i n ^i(x) and 8i n (x)dth(x, t) are at most C{A), completing the 
verification of hypothesis 3. 

Hypothesis 4 is satisfied since supp C < 20 A -1 }. 

Hypothesis 5 asserts that 



dx 



^{eMii + idxhix^))- 1 ^)} 



< C{A)5 



1-j 



and 



dx 



3? {e in {x){l + id x h(x, t^dthix, *)} 



< C(A)6 



1-3 



(4.123) 



(4.124) 



for < j < 2, x G T. 

From( (3~23| we see that 



dx 



{(l + idxh&t))- 1 } 



< c < cs- J 



for < j < 2, x G T. Together with ( |4.121| > and ( |4.122| >, this implies the estimates 



dx 



{e in (x)(i + idxhix^t))- 1 ^)} 



< C(A)5 2 - J 



and 



dx 



5i {9 <n (a:)(l + f))- 1 ^^, t)} 



which are stronger than ( |4.123| ) and ( |4.124| i. This completes the verification of hypothesis 5. 
Hypothesis 6 asserts that, for all x G T, we have 

\to {9 m (x)(l + id w h{x, t))- x d t h{x, t)}\<U {9 m (x)(l + id x h(x, t))- 1 ^)} 



i.e., 



i + (d x h{x,t)) 

i.e., 



^ - 2 |»{(1 - iS^s, t))$M*, t)}| < txt^TT K 1 ~ *))*!(*)} > 



i + (a a fc(s,t)) s 

|5 t ft(a:,i)| < 3f?ffi(a;) +a x /i(x,t)3?5-i(a;) 



(4.125) 



in supp 6 in . 

In view of ( |4.116[ ), ( |4.117| i, ( |4.118| l, estimate ( [4.125) will follow, if we can prove that 



\d t h(x, t)\ < al{x, t) - C'h{x, t) - C'A 



' 4-l T 2 



(4.126) 
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in suppQin. 

However, ( |4.126| l follows at once from ( |3.22| i. This completes the verification of hypothesis 

6. 

Note also that < 9 in < 1, hence hypothesis 7 is satisfied. 

Thus, all the hypotheses (1),...,(7) of the lemma (2.6 1 hold here, with constants C(A). Apply- 
ing that lemma, we obtain the following inequality. 



J? 



Cer+(t) 



(^(C,*)) a + (4(C,*)) 2 



Ch(^C,t)-CA- 1 T 2 



Ar +(t) nC) 



+ 



(i$/i(KC,t))F'(0}<SiC 
^(C)I 2 ^C- 



Adding this to ( |4.120| i, we obtain our basic lower bound 

-27r4(C,t) 



Ch(^RC,t)~CA- 1 T 2 



Cer+(t) 



Lte(C,*)) a + (4(C,*)) a 



Ar +(t )^(C) 



+ 



(id t h(^c,t))F'(C)}d^C 

>-C(A)J2 I \F(Q\ 2 d®C 
± */Cer±(t) 

This estimate holds for is [r 2 , r] . 

Together with ( 4. 108| > and ( 4.94 1, this gives the estimate 



1 d 

2dt 



Cer + (t) 



|^-^)(C,*)| *RC 



>-C\ 2 ~C(A)J2 ! |fl£(z„-z„)(C)| a dRC 
± JCer±(f) 

for t £ [t 2 ,t], under the assumptions made on z, z. 

Recalling from assumption 4 in section 4.2 that the integral on the right is at most A 2 , we 
conclude that 



1 d 

2dt 



\d 4 c (z^-z^((,t)\UViC>-C(A)\ 2 , 



(4.127) 



'Cer+(i) 
ifie [t 2 ,t]. 

This estimate holds provided z„ is a Muskat solutions satisfying 



E / K*M-**)(C)| 2 dHC + £ / |a c 4 (^-^)(C)| 2 «C<A 2 (4.128) 
4. JCer+m . Jcer+m 



and 



Here, 



Cer±(t) 
A<r 2 . 

r±(i) = {ir ± i/i(a:, £) : x £ T}. 



(4.129) 
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4.6.2 The main energy estimate II 



Next, we fix t E [— t,t 2 ]. We prepare to apply lemma (2.3 1 for 6{x) — 1, h(x) = h(x,t), 
ht(x) = d t h(x,t), S = r _1 , and 

-27r4(C,t) 



(^(C,*)) 2 + (4(C,*)) 2 



£=a;+i?l(:c,i) 



C7i(x,i) - CVTV 



(4.130) 



We will check that the hypotheses of the lemma (2.3 1 hold, with constants of the form C{A). 

Hypothesis 1 asserts that h(x, t) > for all x E T, which is immediate from ( |3.37| >. 

Hypothesis 2 asserts that \d t h(x : t)\ < r^ 1 h(x, t) for all x € T, which is precisely ( |3.41| >. 

Hypothesis 3 asserts that, for fixed t, the C 3 — norms of h(x, t) and dth(x, t) are bounded by 
C(A). This is immediate from ( |3~46] > and ( |3T47] >. 

Hypothesis 4 asserts that the C 2 — norm of the function d\ is at most C(A). The C 2 — norm of 
the function cr+(x, t) is at most C; see (3.61 1. The C 2 — norm of the function 

x^Ch(x,t)+CA- 1 r 2 



is at most CA~ 1 ; see ( 3.46 1 



These remarks and (4.130l show that the C 2 — norm of d\ is at most C, proving hypothesis 4. 
Hypothesis 5 asserts that the C 2 —norm of the function 



(1 + id x h(x,t)y 1 (a 1 (x) +d t h(x,t)) 



(4.131) 



is at most C(A). We have just seen th at the C 2 — norm of <j\ is at most C. Moreover, the 
C 2 — norm of dth(x, t) is at m ost CA, by (3.47 1. Since also the C 2 — norm of (1 + id x ti(x, t))^ 1 
is at most C, thanks to (3.46 1, we conclude that the function (4.131 1 has C 2 — norm at most C(A), 
thus proving hypothesis 5. 
Hypothesis 6 asserts that 



on T, i.e., 
on T, i.e., 



8£{(1 + id x h(x,t)y 1 (a 1 (x) + d t H(x,t))} > 

SRcti (x) + d x h(x, t)SSai (x) + d t h(x, t)>0 
5Jcr + (x,t)+ d t h(x, t) - Ch(x, t) - CA- x t 2 



+ d x h(x,t)Qa + (x,t) > 



(4.132) 



onT. 

From p.70| > we have 



n<r+(x,t) > o%(x,t) - Ch(x,t) - Cr 10 
for all x ET; and from ( |3.72| i, we have 

\^a+(x,t)\ < Ch{x,t)+CT 10 



for all x ET. 

Also, from ( |3.46| >, we have 

for all x E T. 



\d x h(x,t)\ < CA- 1 
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The above remarks imply that, for all x £ T, we have 

$ta+ (x, t) + d t h(x, t) - Ch(x, t) - CA^t 2 

+d x H(x, t)$sa+(x, t) 

> al(x,t) + d t h(x,t) - Ch(x,t) - CA^t 2 . (4.133) 

The right-hand side of (4.133 1 is positive; see ( |3.45[ ). Hence, ( |4.133 1 implies (4.132 1, completing 
the proof of hypothesis 6. 

Thus, hypotheses 1,...,6 hold, with constants C(A). Hypothesis 7 holds trivially, since 9=1. 

We may now apply the lemma in section (|2.3|>. We obtain the estimate 



no 



-27Mi(C,t) 



Cer+(t) 



(*i(C,t)) 2 + (4(C,*)) s 



+ (id t h(U(,t))F'(()}d$tt 

>-c{A)t-^ I l^(C)l 2 «C 

+ Jcer ± (t) 



(4.134) 



This estimate holds for t G [— r 2 , t 2 ]. Moreover, we have estimate ( 4.108 I, in which /i should 
be replaced by h since we are working with t G [— r 2 , r 2 ]. Hence, (4.108 1 and (4. 134 1 yield the 
estimate 



1 d 

2dl 



Cer + (t) 



|9 4 (^-^)(C,f)| «C 



> -CX 2 - C'iA)^ 1 V / |^(^ - 2 M )(C)| 2 ^C (4-135) 
± ^Cer±(t) 



fort G [-r 2 ,r 2 ], with 



T±(i) = {a: ± i) : x G T}. 



Recalling from (4.94 1, and from assumption 4 in section d472b that 



we conclude from ( |4.135| ) that 



|9 4 (^-^)(C)| 2 <ffi£<A 2 



1 d 

2dt 



1x2 



Cer+(t) 



W/. -**)(£*)! d$t(>-C(A)T-'\ 



(4.136) 



ifte [-r 2 ,r 2 ]. 

This estimate holds provided is a Muskat solution satisfying 



E / K*M-**)(C)| 2 dRC + £ / |^(^-^)(C)| 2 dKC<A 2 (4.137) 



and 
Here, 



± -/Cer±(t) 
A, r, A, k are as assumed in section [372] 

T ± (<) = {x±ih(x,t) : x £ T}. 



(4.138) 



Our basic energy estimates are (4.127 1,..., (4.129 1 for £ G [t 2 ,t] and (4.136 (,...,(4.138 1 for t G 

h-V 2 ]- 
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4.7 Conclusion 

The previous estimates allows us to achieve the following result. 

Theorem 4.3 Let z(x,t) be a solution of the Muskat equation in the interval t £ [tieasti T ] C 
[— t 2 , t] and let z(x, t) be the unperturbed solution. Assume that z(x, t) satisfies 

• z\ (x, t) — x and z-i (x, t) are periodic with period 2ir. 

• t) is real for Q real. 

• z(C, t) is analytic in £ G Q(t). 

. z(C,t) - (C,o) e H 4 (a(t)). 

• Complex Arc-Chord condition. 

| cosh(z 2 (C, t) - z 2 (w, t)) - cob(*i(C, t) - *i(u>, t))\ > cc A [\mC -w)\\ + mC - w)\} 2 , 
for C w G 

i/ere in f/ze definition of£l(t) we use h(x, t) ift G [t 2 , t] ant/ i) i/i G [— t 2 , t 2 ]. Tnen 

|||z(,t)-z(,t)|| 2 ff4(0(t)) >-C(A)A 2 , 

«/ 1 G [r 2 , r] fl [ti eastl t], and if also 

\\z(-,t) - z(-,t)||H4 (nW ) < A. 

/n addition, 

Jj.U*V' b > ~ -V' ^ll// 4 (o(t)) 
(fi £ [— r 2 ,r 2 ] H [tfeajf, t], and if also 

\\z(-,t) - z(-,t)\\ H 4 {n{t)) < A 

ana! A < r 2 . 

Proof: We have to estimate the quantity: 

d 



i\\z{-.t)~z{-,t)\\ 2 HH nm>~ C i A )^ 1 ^ 



j t \\z(;t)-z(;t)\\ 2 HHn{t)) , (4.139) 



where 



± Jcer±(t) ± Aer±(t) 

(4.140) 

The term in ( |4.139[ ) coming from the second term on the right hand side on ( |4.140| i is bounded 
using ( |4.127[ ) and ( |4.136| >. The rest of terms are lower order terms which are easy to control. 

5 Galerkin approximations and main theorem 

We prove existence of solution of the Muskat equation in H 4 (fl(t)) and the main theorem |TTT| 
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5.1 Conformal Maps Depending on a Parameter 

Let m > 100 and h > be given, and let h(x, t) be a real-valued function on T x [0, 1] with 

Mh(x,t)\<i 

forZ,j < 100m + 100. 

We write c, C, etc. to denote constants determined by m and h alone. 
Let Vo(x + iy) be the solution of the Dirichlet problem 

AV (z) = for -2/i < 3?z < 
V r (a;) = forseK 

V (x-2ih ) = -^3 {cot j ^^ 2 ^ j forxeM . 

Then Vo is 27r— periodic and real-analytic; hence we may regard Vo as harmonic, real analytic 
and 27T— periodic on a strip 

Ho = {-2h -c < 3?z < c }. 
We will solve the Dirichlet problem on the region 

Cl(t; S) = {x + iy : -2h Q - 5h(x, t) < y < 8h(x, t)} 

for < 5 < ci (small enough ci). 

As an approximate Poisson kernel, we try 



P(x + iy,x + ) = ^-3? 



cot 



(1 + (Sh\x + ,t)) 2 ) - Vo(C) ■ (1 + {5ti{x + ,t)f), 



for x + iy e £l(t; 5), x + e R, where 

( = (x + iy) — (x + + i5h(x + , tj) — i5h'(x + ,t) sin[(x + iy) — (x + + iSh(x + ,t))}. 

Since <5 < ci for small enough c\, we know that ( e H n , so P(x + iy, x + ) is well-defined. 
Suppose / e C m (T). Then 



u t (x + iy) 



P(x + iy,x + )f(x + )dx + 



is harmonic in £!(£, S), and one checks that 

lim u t (x + i5h(x,t) — ie) = f(x) + / Ki(x, x + ,t)dx + 
lim u t (x + i5h(x,t) — 2h i — ie) = / K 2 (x, x + , t)dx + 

£^0+ Jj 



where 



and 



\d l t d J x d^ + K i (x,x + ,t)\ < CS for?, j, k < m and i = 1,2, 
\did* + Ki(x,x+,t)\ < CS for j, k < m and i = 1, 2. 



(5.1) 
(5.2) 

(5.3) 
(5.4) 
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We indicate the main steps in verifying the above. 

We introduce a new variable p, and we study how £ behaves when we take y — 6h(x,t) — p. 
We work on the domain 

{(x,x + ,p,S,t) \x-x+\ <7T,8e [0,c], t € [0,1]} 

for a small constant c. 

Let us write G±, G2, etc. to denote functions on this domain having at least 10m derivatives 
of absolute value at most C. 
Then 

( = Z-iSti(x,t) sin(Z) 



with 



First of all, note that 

( = Z- (l-iSti(x + ,t) 



Z = [x + iSh(x, t) — ip] — [x + + iSh(x + , t)]. 
sin{Zy 



z 



= [0 - x+) - ip] 
and therefore 

uniformly in x, x+, p, 5, t. 
Next, by writing 



we find that 



moreover, 



. h(x,t) - h(x+,t) 

1 + 10 -. r ; 

(X - X + ) + ip 



l-iSh'(x+,t) 



sin(Z) 



\C\ = \(x-x+)-ip\-(l + 0{S)) 
"sin(Z) - Z 



sin(Z) = Z + Z 6 



Z 3 



( = Z - iSh'(x+,t)Z + SZ 3 G! 



Z = (x — x + ) ■ [1 + i8h'(x+, t)} — ip + S(x — x + ) Gi- 

Hence, 

C = [1 + (6h'{x + ,t)) 2 } ■ (x - x+) - Spti (x+,t) - ip + S(x - x+) 2 G 3 + Sp 2 G 4 
Next, we refine ( |5.7| i in the case \x — x + \ < pi . In that case, we have 

\S(x - x+) 2 G 3 + S P 2 G 4 \ < CSpi 

and 



\5ph'(x+,t)\ < CSp, 



while 



\[l + (6h\x + ,t)) 2 ](x-x + )-ip\ >cp. 
Therefore, ( |5.8| l implies that 

ICI = |[1 + {5h'(x+,t)) 2 ]{x - x + ) - 5pti{x +1 t) - ip\ ■ (1 + 0(5p^)), 

uniformly in x, x + , p, 5, t, under the assumptions \x — x + \ < pi . 



(5.5) 
(5.6) 



(5.7) 



(5.8) 



(5.9) 
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Now we are ready to analyze 9cot(£/2). For |5ft£| < |tt, we have 

with H analytic and real-valued for real (. 

On the other hand, £ is a C 90m+90 — smooth function of x, x+, p, 8, t, equal to (x — x+) when 
p = 8 = 0. Therefore, 

^ cot u) =-|^+pG 5 +5G 6 . (5.10) 
From ( |5.7| i,..., f5.10| > and the properties of Vo discussed above, it is now a routine task to check 

Thus, u t (x + iy) is approximately equal to / on the upper boundary of fl(t; 8), approximately 
zero on the lower boundary of S) and (exactly) harmonic in f2(t, 8). 
Moreover, Ki(x, x+, t) (i = 1, 2) are given by explicit formulas. 

There are analogous explicit formulas for a harmonic function on f2(t; 8) approximately equal 
to a given function g on the lower boundary of fi(i; S) and approximately equal to zero on the 
upper boundary. 

Our explicit formulas, together with a Neumann series, show that the Dirichlet problem 



V(z, t) harmonic in f2(t; 8) for each fixed t 

V(x + i5h(x,t),t) = 1 

V(x - i8h(x, t) - 2ih ,t) = -1 

admits a solution V that is C m on the region 

E(S) = {(z,t) : zeJl(^) d ~, te [0,1]}. 

In particular, the above Neumann series converges, since we take 8 < c\ for a small enough 

Cl- 

Taking a harmonic conjugate U of V, we may suppose U also belongs to C m on the region 
S(£). 

The analytic function $ t = U + iV is a conformal map of f2(f ; 8) to the strip {|Sz| < 1}. 

Moreover, $ t depends C m -smoothly on t, and + 2tt) = $t(z) + A(t) (all z e 8) 
for a smooth, real-valued function X(t). 

A trivial rescaling now shows that f2(t ; 8) maps conformally to a strip < by a 

conformal map <f>J, that depends C m — smoothly on < and satisfies + 27r) = 3>f(z) + 27r; 
moreover, is C m — smooth. 

We can now easily remove our small 8 assumptions. 

Let 

= {x + iy : h-(x,t) < y < h+(x, t)} 

for t G [0, 1]; where h + and h_ are real, smooth, and 27r-periodic. 

We assume that h- (x, t) < h + (x, t) for all x £ M, t G [0, 1]; and we suppose that (x, t) 
and h + (x, t) are real-analytic in x for each fixed t. (The real-analyticity assumptions ought to be 



removed, but our h+ and h- happen to be real-analytic. See section 3.2 ) 

Given to € [0, 1], a conformal mapping defined on a neighborhood of f2(io) closure carries 
fi(io) to a strip. This same conformal transformation carries fl(t) to a domain close to that strip, 
whenever t is close enough to to. Rescaling in the t— variable we reduce matters to the case 
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considered above. Hence, O(t) is mapped conformally to {|3z| < h)(t)} by a conformal map 
z i-» $(z. i), that depends C" l -smoothly on (z, t) for z e r2(t) closure , i near t . We have also 
$(z + 27r,i) = $(z,t) + 2?r. 

Thus, we have our smooth family of conformal maps, defined in a small neighborhood of 
any given to- These maps carry the upper boundary (respectively, the lower boundary) of fi(t) 
to the upper boundary (respectively, the lower boundary) of H{t). Since such <i> f are uniquely 
determined up to translations, it is trivial to patch together our results on small i-intervals. 

Thus, we obtain the following result. 

Lemma 5.1 Let m > 100 be given. For t G [0, 1], let 

Sl(i) = {x + iy : h-(x,t) < y < h + (x, t)} 

where h + , h- are smooth, real-valued functions on M x [0,1]. We suppose that h-(x,t) < 
h + (x,t) for all (x,t) 6lx[0, 1], and that h + (x,t), h^(x,t) are 2tt— periodic and real-analytic 
in xfor each fixed t. 

Then there exist a C m —smooth, positive function h*(t) defined on [0, 1], and a C m — smooth 
map 

$ : {(z, t) e c x [o, i] : z e n c,osure (t)} -> c, 

such that, for each fixed t £ [0, 1], z t— > maps f2(i) conformally to {|3CI < ^"(^)}- 

Moreover, + 2n, t) = $(z, t) + 2it for each (z, t) in the domain of '(f), 

5.2 Changing Coordinates 

Let h(x, t) be a real-analytic, positive function on T x [t — 5,t + 5], for h(x + 2ir, t) = h(x, t) 
for all x, t. Let 

r±(t) = (CeC: SC = ±W,i)} 
n(t) = {C e C : |9fC| < h($C,t)}. (5.11) 



Then, by lemma 5 . 1 there exist a smooth positive function h (t ) on [to — 8, to + 8] and a smooth 
map 

<& : {(C,t) : C e ^ dosm ' e W, I* — *o I = I^Cl < M*), I* " *d| < <0 (5-12) 
with the following properties: 

$(C,*) = (C,*) = (*(£*),*) (5.13) 

for each (£, i). The map 

C->0(C,t) (5.14) 
is a biholomorphic map from Q(t) to the strip {|3Cl < ^(*)}> f° r eacn ^ 

0(C + 2tt, t) = <ft(C, t) + 2tt (5.15) 

: r+(t) ^ {3C = ft(t)} (5.16) 



for each fixed t. 

The inverse of $ is the map 



<p(£,t) is real for real (. (5.17) 

*(C,t) = (V-(C,t),*) = (C,0. (5.18) 
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We would like to solve the Muskat equation 

sin(zi(C,i) - zi(w,i))[(? C 2 Al (C,<) - d w z^(w,t)] 



u>er+(t) 

for C € t E [to — 6', t ], with 



cosh^ (C, t) - z 2 {w,t)) - cos(«x(Ci*) _ 



with 



t) holomorphic in £ G fi(t) for fixed t, 



2„(C. *) - C*i, /< e i? 4 (r± (t)) for fixed t, 



where Sij is the Kronecker delta, with 

Zi(Cj i) — C an d ^2(Ci periodic with period 27r, for fixed t. 

z p (C, i) is real for real £, any i G [to — S',to], 
and satisfying the initial condition 

z„(£*o) = *£(C)> 



(5.19) 

(5.20) 
(5.21) 

(5.22) 
(5.23) 

(5.24) 



where z°(-) satisfy (5.20 1,.. .,(5.23 I, as well as the chord-arc condition 

| cosh(z 2 ((,t ) - z 2 (w,t )) - cos(zi(C,to) - 2i(w,i ))| 



(5.25) 



>c CA [||K(C-«;)|| + |9f(C-ti;) 

forC, w G fl closure (t a ). 

We assume that z?(£) satisfy a Rayleigh-Taylor condition, to be explained below. Our goal 
is to transform the Muskat problem (5.19 1,... ,(5.25 1 to an equivalent problem in (C,t) coordi- 
nates, and to show that the Rayleigh-Taylor condition for (5.19 (,...,(5.25 i implies an analogous 
Rayleigh-Taylor condition in ((, t) coordinates. 

For ( 5. 19 !,...,( 5.25 1, we introduce the 'Generalized Rayleigh-Taylor function' 

RT((,t) = 
-2^ c Z!(C,i) 



If 



/ 



P.V. 



-[i + tdxh^Ct)]- 1 

(d c zi{£,t)) 2 + [d c z 2 {(,t)y 

sin(zi(C,f) - gifovQ) 

cosh(z 2 (C, t) — z 2 (w, t)) — cos(zi(£ t) — z 1 (w, t)) 



)dw 



wer + (t) 



+ 



^(RC.tjJfl + i^W.*)]" 1 



forCGr+(t). 

Our Rayleigh-Taylor assumption regarding ( |5.19| l,...,( |5.25[ ) is as follows 



(5.26) 



(5.27) 



(strict inequality) for ( G r+(t). Let us check that assumption ( 5.27[ ) holds for z M (C, t) as in the 
previous sections, using either of our previous functions h(x, t) or H(x, t) for the function h of 
this section. (If we use h(x, t), then we assume that t G [t 2 , t]; if we use H(x, t), then we assume 



70 



that t £ [— t 2 , t 2 ]). We wi ll a lso as sume that A < A 1 t 2 . After we check this, we return to the 
task of transforming (5.19 (,...,(5.27 i into (£, t)— coordinates. 



Let's start with the case of the function h(x, t) (rather than h{x, t)). 
Then, for Q — x + ih(x, t) e T + (t), we have 



-27r9 c ^i(C,t) 



-2nd x zi(x, t) 



and 



(a cZl (c,t)) 2 + (d c z 2 ((,t)) 2 {d xZl {x,t)y + {d x z 2 {x,t)y 

—2Trd x zi(x,t) —2izd x Zi{x,t) 



< Ch(x,t), 



< CA-'t 2 



(d x z 1 (x,t)) 2 + (d x z 2 (x,t)) 2 {d x z x (x,t)) 2 + (d x z 2 (x,t)) 2 
/ —Om-Fi^Jr A 

[l + id^xA)}- 1 



Therefore, 

> -Ch(x,t)~CA- 1 r 2 



{d^^t)) 2 + {d x z 2 {x,t)) 2[ 

-2-Kd x z x (x, t) 



-Ch(x,t) - CA-\ 2 + ZTT : \ rtW ft + (d x h{x,t)) 2 y\ 

{d x z 1 [x,t)y + {o x z 2 {x,t)Y 



(5.28) 



Recall that 



5(C,t) 



sin(zi(C,<) - (to,*)) 



i»er + (() 



cosh(z 2 (C ; £) ~ z 2(w, i)) — cos(zi(C, i) — zi(w, t)) 



(da.iu)) 2 + (d ( z 2 (ct)) 2 

Thanks to our observation in section (12.21) that 



cot 



C — w 



)} 



dm. 



(5.29) 



P.V. 



, C - w 

cot — 



eku = 0, 



we have that 

a(c,t) = p.v. 



sin(2i(C,t) - zi(t«,i)) 



ujer + (t) 



cosh(z 2 (C,0 ~ z 2 {w,t)) - cos(«i(C,*) - zi(w,t)) 



(5.30) 



Comparing (5.29 1 with ( 4.58| > and (4.62 1, we see that a is as in section (4.4 1. Hence, (4.75 i 
together with ( |4.54| >, ( |435| , tell us that 



1,2 



withx = c g r+(t). 



Thanks to (4.56 1, (4.57 1 and (4.75 1, analogous estimates hold also for H(x,t) in place of 
h(x, t). 

Consequently, 



{5(C, t)[l + id x h(x, t)}- 1 } > -Ch(x, t) - CA 



-\ 2 . 
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Recalling ( |5.30| l, we conclude that 

sin(zi(C,i) - z 1 (w,t j) 



3 < P.V. 



™er + (() 



cosh(z 2 (C,0 ~ z 2 (w,t)) - cos(zi(C,£) ~ zi{w,t)) 



)dw[l+id x h(x,t)}~ 1 



1^2 



> -Ch{x,t) -CA- l t 



(5.31) 



for a; = M(, ( G 

Analogous estimates hold for 7i(a:, t) in place of /i(x, t). 
We note that 



1 + (d x h{x,t)y 



(5.32) 



Putting ( |5.28| ), ( |5.3Q[ ), < |5 .3 1 1 >, ( |5.32| i into definition ( |5.26| l, we see that, for C G r + (t), a; = 3?C 
we have: 



JzrtCtj^ci + ^i.t)) 2 )- 1 

x -Ch(x, t) - CA^t 2 + 



-2ird x z x {x,t) 



{dxz^x.t)) 2 + (d x z 2 {x,t))'- 



+ d t h(x,t) 



(5.33) 



The term i n squ are brackets in (5.33 1 is strictly positive, thanks to ( |3.21| i. 

Thus, ( 5.27 1 hold s for t G [t 2 , r] , under the assumptions of the previous sections. 

Similarly, (5.27 1 holds also for t G [— t 2 ,t 2 ], with h(x,t) in place of h(x,t), under the 
assumptions of the previous sections. 



We have succeeded in verifying the Rayleigh-Taylor condition (5.27 1 for our Muskat so- 
lutions fr om se ction |4] We now return to the business of transforming our Muskat problem 
( 5.19 i,. ..,( 5.25 1, into (C, *)- coordinates. 

We make the change of variable 



so that also 

Recall, t = t. We define 



( = <p((, t), w = <j)(w, t), for fixed t, 
C = V(C> *)> w = VK^j *)> f* or fixed t. 

2„(C,t) = far = 



Note 

Similarly, 
Thus 
and 
Also, 



0^ M (£t) = 9^(C,i)a c z M (C,<) for (C,t) - *(C,**). 
d^z^(w,t) = d li tp(w,t)d w z fl (w,t) for (iu,t) = \I>(u>,t). 

9 c z M (C,i) = 9 c ^(C,«)^(C,t) for (C,t) = $(C,t); 
d w z^(w,t) = d w 4>(w,t)d^z M (w,i) for (u>,£) = $(C,i)- 
du> = d^ip(w,t)dw 



(5.34) 
(5.35) 

(5.36) 

(5.37) 
(5.38) 
(5.39) 
(5.40) 
(5.41) 
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for fixed t. 

The chain rule shows that for ((, t) = t), we have 

d t z^t) = d t [z M o $] (c,t) = ft{i M (0(C, *),*)} = dMtt) + [d t 4>{^t)]df Zll {lt) 

= dMl i) + W) o *(C, t)]S^(C, *)■ (5-42) 

In view of the above remarks, our Muskat equation ( |5.19| l is equivalent to the following equation 
in (£, t)— coordinates 

d t z^c,t) + p t 0) o *(c,t)]a e ^(c,t) - 

sin(zi(C,<) - zi(w,i)) 



cosh(z 2 (C,i) - z 2 {w,t)) - cos(li(C,t) - li («),*)) 
s™=/i(t) 

x {[(0 C $ o *(C,t)]^^(C,t) - [(cW) o *(tB,t)]a(D« M («»,«)} d^{w,t)dw (5.43) 
Since f) and t) are inverse functions, we see that 

[(d w ct>)oy(w,i)]-[d^(w,i)] = l. 

Hence, defining 

B((,w,i) = [(d ( <t>) o*(C,<)][9^(*,t)] - 1, (5.44) 

we see that 

B(C,w,t) is smooth on < f e [i Q - 5,t + <5]} (5.45) 

73(C, w), t) is holomorphic in (£, t2>) £ < ft(t)} for fixed t (5.46) 

and 

B(C,C,f) = 0. (5.47) 



Remark 5.2 Note that B does not depend on our Muskat solution; it is entirely determined by 
the function h(x, t) and the coordinate change (£, t) = t). 



Thanks to ( |5.44) ,...,( |5.47 1, we can write ( 5.43[ ) in the equivalent form 



[(M) o *(££)] + / M^-^t))B(lw,t) ;) ^ _ 

7 cosh(z 2 ((,t) - z 2 (w,t)) - cos(li(C,t) - 2i («>,£)) " 
cosh(z 2 (C,t) - Z 2 (u>,t.)J - cos(2i(C,t) - zi(w,i)) 

9u>=/i(t) 

(5.48) 

This holds for 1 3?C | < ft(t). 
We define 

i(C,t) = -[(M°*(C,*)] (5-49) 
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for |3CI < h(i), i e [t -S,t + 5}. 
Thus 

t) is smooth in (£, i), and holomorphic in ( for fixed t. (5.50) 



Remark 5.3 Note that A(£, t) does not depend on our Muskat solution. It is determined entirely 
by h(x, t) and $(C, t). 



We may rewrite ( |5.48| l in the equivalent form 



--h(i) 



sin(zi(C,t) - Zi(w,i))B(C,w,i) 

cosh(z 2 (C,i) - z 2 {w,i)) " cos(li(C,t) - Zi(w,t)) 



sin(£i(C,<) - £i(tw,t)) 



Qui=/i(t) 



cosh(z 2 ((,t) - z 2 (w,t)) - cos(fi(C,i) - fi («;,£)) 



i) - dujZ^w, i)]dw 
(5.51) 



We call < g3j) the TRANSFORMED MUSKAT EQUATION or 'TME'. 

Our additional conditions ( |5.20[ ),...,( |5.25| l are easily seen to be equivalent to the following: 



t) is holomorphic in {|S*C| < h(t)} for fixed £. 
t) ~ fa, ,Gi? 4 = ±M*)}) for fixed t. 

^i(C) *) — C ar, d -^2 (C; i) 316 27T— periodic for fixed t. 
z"n(C, t) is real for real (. 



(5.52) 
(5.53) 

(5.54) 
(5.55) 



The initial data £/i(C,io) = ^2(0 sre given and satisfy (5.52 (,...,(5.55 1, as well as the chord-arc 



condition 

| cosh(z 2 (C,to) - Za(w,to)) - cos(5i(C, < ) - Zi(w,t ))\ 
>c[m(-w)\\ + \Z((~w)\} 2 

for \<S(\, \<3w\ < h{t ). 
Here, of course, 

^(C)=^°*(C,*o). 

To deduce the equivalence of ( 5.23[ ) and ( 5.55 ), we use the fact that 



A((, t) is real for real £; and 
B((, w, t) is real for real io 



(5.56) 



(5.57) 



(5.58) 



and we change the contour in (5.51 1 to = 0}, i.e., our contour may be taken to be the circle 



T. 



We use also the fact that 



A((, t) is 2ir— periodic in ( for fixed t; and that 
.B(C) w, t) is 27T— periodic in both £ and u) for fixed t. 



(5.59) 
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Observe that ( |5.5 1 [ >,... ,( |5756] > are equivalent to ( |5.19) ,...,( [5!25) , 

It remains to formulate a 'Generalized Rayleigh-Taylor Condition' for ( |5.51) , ( |5.56) , and to 
prove the equivalence of that condition to ( |5.27| i. 

We define (for 3C = Hi)) 

RT(U) = 



-2-Kd t h{Q,i) 



I 



(0 f fi(C,t)) 2 + (d^HCW 

sin(zi(C,i) - zi(w,t)) 



Qw=h(t) 



cosh(z 2 (C,i) - z 2 (w,t)) - cos(5i(C,<) - z^w.t)) 



{B(£,w,t) + l)dw 



■Z{A((,t)} + h'(t), 



(5.60) 



We will see that ( 5.60| l is the natural Rayleigh-Taylor function for the problem ( 5.51 !,...,( 5.56 1. 
Our 'Generalized Rayleigh-Taylor Condition' for ( |5.51| l,...,( |5.56| l is 



RT((,t Q ) > 



(5.61) 



(strict positivity) for all ( such that 3£ = h(t ) 



We hope to produce a solution to (5.51 (,...,(5.56 1, provided our initial data zj; satisfy the 



chord-arc condition (5.56 1 and the 'Generalized Rayleigh-Taylor condition' (5.61 1. 

This will be done in a later section; for the moment, we simply note that our solution to 



(5.51 (,...,(5.56 1 will later be shown to exist in a short time interval [to — 6' , to] where 5' will be a 



5.61 



5.18 



small enough positive number determined by the constants in our assumptions on the initial data. 
For the rest of this section, our task is to prove that condition ( |5.27[ ) implies condition 
To do so, we first make a few elementary remarks on the maps <1>, ^, (j>, ip in (5.12 1,... 

Then we return to ( |5.27| > and ( |5.61| l. 

First of all, since £ H> t) maps P + (t) to {3£ = h(t)} for each t, we have 

d c 4>(t,t) = p(Ct)[i + id x h(x,t)]- 1 

for x = 3£ = h(x, t); here, 



p(C,t) >0 



for x = Ji^, 3£ = h(x, t). The operator 

d 



d t + [d t h(x,t)} 



corresponds under $ to 8$ + [h'(t)] 



8 



for x = $t(, Q?C = Kx, t), (C, t) = $(C, t); 



here 



t) is real whenever ( 5.65 i holds 



(5.62) 
(5.63) 

(5.64) 
(5.65) 
(5.66) 



Indeed, (5.62 1, (5.63 1 hold because r+(i) maps to {Q( = h(t)} with positive orientation (for 



each fixed t); and ( 5.64[ i,...,( |5.66| l hold because a moving particle that stays on the upper boundary 
r+(£) as t varies will stay on the upper boundary S£ = h(t) when viewed in the (£, t) coordinate 
system. 



Let t) be a holomorphic function of ( for fixed i; let F = F o vp. Then (5.64 1 gives 

8 t F((, t) + [d t h(x, t)]id c F((, t) = d;F((, t) + h'd^FiC, t) + b(C> i)d^F(C, t), (5.67) 
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whenever ( |5.65| ) holds. 

On the other hand, since F((, t) = F(<f>(£, t),t), the chain rule gives 

d t F(c,t) = dfF(c,i) + [Mic, t)]d £ F((,t) 



whenever (£, t) = $(£, t), in particular, whenever (5.65 1 holds. 
Subtracting ( |5.68[ ) from ( |5.67| i, we find that 

[id t h(x,t)]d c F(<;,t) = [ih'(t) + h(c,i) - d t <t>((,t)]d £ F(c,i), 



(5.68) 



(5.69) 



whenever ( |5.65[ ) holds. 
The chain rule gives 

d c F((,t) = [d c <t>(t,t)]d i F(t,t) 

whenever (5.65 i holds. Putting this into ( 5.69| l and recalling that F((, t) may be any function 
holomorphic in £ for fixed i, we see that 

idth(x,t)[d c <f>({,t)] = idfhd) - b(c,t) - 9 t <£(c,t) 

whenever (5.65 1 holds. Recalling ( |5.62[ i, we now see that 

id t h(x,t)p(C,t) 



Next, we note that since z M (C, t) — ^(^(C, i),t), the chain rule gives 

5c^(C,*) = [^(C,*)]^(C.*) 



(5.70) 



whenever (£, t) = <!>(£, t), in particular, whenever (5.65 1 holds. 
Consequently 



-27ra c «i(C,t)[a c ^(C,t)] 



2TTd £ h{C,t) 



(d c zi(t,t)y + (d c z 2 (c,t)y (^^(ct)) 2 + (^£ 2 (c,i)) 2 

whenever ( |5.65[ ) holds. Recalling ( |5.62| i, we conclude that 

-27rS t «i(C,t) p((,t) _ -27rfl^i(C,*) 



(a C 2i(C. <)) 2 + (^(C, *)) 2 1 + i5 *M*> «) (C, 0) 2 + (0^2(C, 0) 2 

whenever ( |5.65[ ) holds. 



(5.71) 



Next, substituting the definition ( 5 .44 1 of B(£,w,i), and changing variable by £ = (f>(£,t), 
w = 4>(w, t), t — i, we find that 



P.V. 



sin(ii(C,i) - zi{w,i)) 



Qw=h(i) 



cosh(i 2 (Ci*) - z 2 (w,t)) - cos(£i(C,t) - zi («>,*)) 



(P(C,to,i) + l)cM> 



P.V. 



sin(£i(C,<) - zi(iu,t)) 



ai«=/i(t) 
= P.V. 



cosh(52(C,t) - h{w,t)) - cos(ii((,£) - zi(io,t)) 
sin(zi(C,t) - zi(iu,t)) 



«jer+(i) 



cosh(z 2 (C, f) — 22(10, t)) - cos(zi(C, t) — zi(w, t)) 



(d c <j>)oy((,t)][d a i>{w,t)]dw 

)dwd c <f>((,t). 



76 



Recalling ( 5.62 1, we conclude that 

sin(li(C,t) - zi(w,t)) 



P.V. 



Qw=h(i) 

1 + id x h{x,t) 



cosh(.z 2 (£,t) _ h{w,t)) - cos(fi(C,t) - zi(w,t)) 



(B((,w,i) + l)dw 



P.V. 



wer + (t) 



sm(zi(C,t) - z 1 (w,t)) 

cosh(z 2 (C,0 - z 2 {w,t)) - cos(zi(C,i) - zi(w,t)) 



)dw 



(5.72) 



whenever ( |5.65| ) holds. 

Next, we apply ( |5.49| > and ( |5.70| i. Thus, whenever ( |5.65[ ) holds, we have 

A(c,t) = -Wet) = ^Vlf^ lh '^ ^ 5 - 73 > 

1 + io x n(x, t) 



Now, substituting \5.1\\ , \5.12\ and ( |5.73| into the definition \5.60\ of RT((,t), we find that, 
whenever ( |5.65[ ) holds, we have 

-2nd cZl (C,t) p(C,t) 

(d (Zl (c,t)r + (a c z 2 (c,t)) 2 i + id x h{ x ,t) 



If 



+3 



P(C,*) 



1 + id x h(x, t) 



P.V. 



sin(zi(C,t) - zi(io,t)) 



-)dw 



cosh(z 2 (C,0 - z 2 {w,t)) - cos(zi(C,t) - 

— — 6(C,*)\ + (5.74) 



™er+(t) 
id t h(x,t)p(C,,t) 



1 + id x h(x, t) 

Note that the /i'(t) terms cancel, and that b((,t) does not contribute to the right-hand side of 
( [5.741 ), thanks to ( |5.66] >. 

Recalling also ( 5.63[ ), we see that ( |5.74[ ) may be rewritten in the equivalent form 

RT((,t) = 
-27ra c «i(C,t) 1 



p(C,t) 



1 + id x h(x, t) 



P.V. 



(d cZl ((,t))2 + (9 c z 2 (C,t)) 2 l + tfl a /i(i,t) 

sin(zi(C,f) - 

cosh(z 2 (C,i) - z 2 (w,t)) - cos(zi(C,t) - zi(w,t)) 



)dw 



5R 



d t h(x,t) 
1 + id x h(x, t) 



whenever ( |5.65[ ) holds. 

Comparing this equation with the definition (5.26 1 of i?T(C, t), we conclude that 

RT((,i) = p(t,t)RT(t,t) (5.75) 

whenever ( |5.65[ ) holds. 

Given (, t with = W), we set (£,t) = *(£,*) and x = 5?C- Then Q?C = f), 
a; = 3?C and (£, t) = $(C, *), i.e., (5.65 1 holds. Therefore, (5_75 1 applies. We recall from ( 5.63 i 
that p((, t) is s trictly positive. If (5.27 1 holds for ^(C, t), then RT((, t) is also strictly positive. 
Therefore, by (5.75 I, we have RT((, t) > (strict inequality). This is the desired 'Generalized 
Rayleigh-Taylor Condition' ( |5.61| l. Thus, as claimed, ( |5.61| > follows from ( |5.27[ >. 
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5.3 Generalized Muskat 



In this section we prove local existence for the 'Transformed Muskat equation' ( |5.51| l,... 
Let's start with some definitions. 



Define 



Let h(t) be a positive, smooth function defined for t 6 Iu m e = [to — to). 

r±(t) = {(ec : 9fC = ±M*)}- 

tl{t) = {C G C : |9fC| < h(t)} 

n = {(C,*) : C G m ClOSUre , t G Iume] 



closure 



e n(t) closure , t e i time } 



Assume: 



and 



Assume that 



Assume that 

Assume that 
such that 



Let A(£, t) be a smooth, complex-valued function on Q. 

A((, t) is real when ( is real; 

A(C + 2n,t)=A(C,t); 

£ i — y t) is holomorphic on fi(t) for fixed £. 
Let £?(C, be a smooth, complex-valued function on + . 

-B(C, w, i) is real when £, w are real; 
B(C + 2vr, iu, *) = B((, w + 2vr, t) = B((, w, t); 
(£, w) H> -B(C, w, t) is holomorphic on x for fixed t; 
B((,(,t) = 0. 

For /x = 1, 2, let z°(C) be a holomorphic function on Sl(i )- 

z?(C + 2^) = zO(C) + 27r and z 2 °(C + 2tt) = z 2 °(C); 
z u(C) i s rea l whenever £ is real. 



z°(C) extends to a continuous function on f2(io 



. closure 



(*?(0-0 r , ^2°(C) 

1 ±(*o) 



r±(t ) 



belong to i? . 



We assume also that z° satisfy the chord-arc condition: 



((535). 

(5.76) 

(5.77) 
(5.78) 
(5.79) 
(5.80) 
(5.81) 

(5.82) 
(5.83) 

(5.84) 
(5.85) 

(5.86) 
(5.87) 
(5.88) 
(5.89) 
(5.90) 

(5.91) 
(5.92) 

(5.93) 

(5.94) 



| cosh(z 2 °(C) - 4{w)) - cos(z°(C) - 4(w))\ > c C A [\mC -w)\\ + |9f(C - w)\} 2 (5.95) 

for n(t ) closure . 

Also, we define the Rayleigh-Taylor function for the initial data as follows. 
For £ 6 r_|_(to), we define 

i?T°(C) = 
-27r(z?)'(C) 



((^)'(C)) 2 + (fe )'(C)) 2 
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J cosh(z^(C) - z%(w)) - cos^C) - z$(w)) 
wer+(t) 

+ ^A(C 7 t ) + h , (t). (5.96) 
We assume that the following Rayleigh-Taylor condition holds for the initial data: 

RT°(C) > for all C & T+(t ) (5.97) 

Note that we assume strict inequality here; we do not allow RT° (£) to vanish at any point ( e 

r+(to). 

Next, we introduce the notation of a 'controlled constant'. 

A controlled constant is a positive real number determined entirely by the following data. 
An upper bound for the C 4+100 — norm of h on I t ime- 
A positive lower bound for h(t a ). 

An upper bound for the C 4+1 °— norms of ^4(£, t) and B((, w, t). 

An upper bound for the H 4 —norm of — ( and (C) on r±(£o). 

The constant cca in the Chord- Arc condition. 

A positive lower bound for the Rayleigh-Taylor function RT°(() on r + (io). 

We write c, C, C, etc. to denote controlled constants. 

These symbols may denote different constants in different occurrences. 

We assume that the length S of the time interval I t ime above, is less than a small enough 

controlled constant. (5.98) 

Under the above assumptions, we would like to prove the existence of a solution to the following 
GENERALIZED MUSKAT PROBLEM. 

Find C 1 functions z M on f2 satisfying the following conditions: (5.99) 

(C, t) is real whenever ( is real. (5. 100) 

zi(C, t) — C an d z 2 (C, t) are 2ir— periodic function in £ for fixed t. (5.101) 

z M (C, t) is holomorphic on fl(t) for fixed t. (5.102) 

Zi(C, t)-( and z 2 ((, t) belong to H 4 (T ± (t)) for fixed t, (5.103) 
with norm bounded by a controlled constant. 

*m(C>*o) = *£(C) (5.104) 

| cosh(z 2 (C,t) - z 2 (w,t)) - coB(*i(£t) - *i («;,t))| > C' CA [||R(C -tu)|| + |9(C - w)|f 

(5.105) 

for (C, to, t) e o+. 

For t e /time, C € r+(i), define 

JZT(C) = 

f -2 7 r(a c z 1 )(c,t) 1 
l((^ 1 )(C,t)) 2 + ((^ 2 )(C,i)) 2 i 

J cosh(z2w{(,t) - z 2 {w)) - C0B(*i(£t) - ziH) L V ; J 
t«er + (t) 



+9< 
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+ <3A{(,t) + ti(t). 



Then 



RT((, t) > for t G lame, C e r+(t). 
For (£, t) G fi, the following equation holds: 



(5.106) 
(5.107) 



A(C,*) 



/ 



mer + (i) 



sin(zi (C, t) - zi (w, t))-B(C, w, t) 
cosh(z2(C,*) - z 2 {w,t)) - cos(zi(C,t) - Zi(ti;,i))' 



-dw 



►0c**.(C,*) 



sin(zi(C,t) - 0i 



[^C^(C,t) - flu, 2^(10, i)]<iifl. 



cosh(z 2 (C,i) - 2 2 (u>,t)) - cos(zi((,t) - Zl(l0,t)) 

(5.108) 

In the above integrals, the integrands are holomorphic in w G for fixed ( and i; in particular, 
the integrands have no poles. Hence, we may change the contours in these integrals from T + (t) 
to T, without changing the values of the integrals. 

In the next section, we prove existence of solutions to the above GENERALIZED MUSKAT 
PROBLEM. 

Thanks to the results in section ( |5.2| i this yields the desired local existence of solutions of our 
original Muskat equation, with desired holomorphic and H 4 properties, etc. 



5.4 Galerkin Approximation 



Let h(t), r± (t), t , 6, Q(t), fl, A((, t), B((, w, t), z"(Q be as in section (5.3 1. (5. 109) 



For a large enough positive integer N, we define z^ (£) from z®(Q, by using the projection 



n 



-N 



We define z\^\Q by stipulating that 

zi JV1 (C)-C = n Ar [z?(0-C] 

and 

zf ] (0 = n N [z° 2 (0}- 

For N l arge en o ugh, th e functions zjf^ satisfy the chord-arc and Rayleigh-Taylor condition 



(5.110) 



(5.111) 



(5.112) 



(I5T05J, (|5.106[), (|5.107[) with controlled constants. Also, since z[ N] (() and z[ N] (£) - ( are 



trigonometric polynomials, those function are holomorphic on fl(to), they extend to smooth func- 
tions on the closure of fi(to), and their H 4 — norms on r to are at most the H 4 — norms of z 2 and 
z \ (C) — C> respectively, and are thus dominated by a c ontro lled constant. 

Consequently, the conditions imposed in section (5.3 1 for z° hold also for z\f\ with con- 



trolled constants, provided N is large enough. From now on, we suppose that N is large enough 
that this holds. 



In place of the Generalized Muskat equation (5.108 I, we will solve its Galerkin approxima- 
tion. That is, we will solve the equation 



d t z ll {C ) t) = u N [j(-,t)}(c,t) > 



(5.113) 
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where 



wer + (t) 



J((,t) = 

sin(zi(C,t) - zi(w,t))^(C, w,t) 

cosh(z 2 (C,t) - z 2 (w,t)) - cos(zi(£,t) - z\{w,t)) 



dw 



sin(zi(C,t) - zt(w,t)) 



cosh(z 2 (C,t) - z 2 {w,t)) - cos(z!(C,t) - Zi_{w,t)) 



wer + (t) 

We impose the initial condition 



^(C,io) = 4 2V] (C) 



[fyzniCt) - d w z^(w,t)]dw 
(5.114) 
(5.115) 



Here, (£ , t) depend on N, though our notation does not indicate that. 
We hope to solve ( |5.113) , ( |5.114) , ( |5.115| l for t G /« me . 

Note that ( |5.113[ ), ( [5.114) , ( |5.115| l is simply a system of ODE's. We solve these ODE's, 
starting at time to, and proceeding backwards in time, until either we reach time to — 6 (the left 
endpoint of the interval Itime), or else the conditions (5.99 1,.. .,(5.107 1 (with suitable controlled 
constants) are not satisfied. 

We will see that we get all the way back to time to — 5. Thus, we obtain a solution of 
( |5.113[ l, ( |5.114[ i, ( |5.115[ ), satisfying conditions ( |5.99]>,..., ( f5. 1Q7| > with controlled constants; these 
controlled constants are independent of N, so by ( 5. 103[>, we m ay pick a sequence N v — > oo for 
which the corresponding solutions of (5.113 I, (5.1 14 1, (5.1 15 i tend weakly to a limit z^(£, t) in 
H* (r± (t)) for fixed t (rational). 

The resulting z M (£ t) is then easily seen to solve the GENERALIZED MUSKAT PROBLEM 
in section \53\ . 

It all comes down to an a priori estimate 



d 
dt 



|^(C,*)| dC>-C, 



(5.116) 



Cer±(t) 



to be proven as long as our solution to ( |5.113| >, ( |5.114[ ), ( |5.115| l satisfies ( |5.99| l,...,( f5.107| >. Note 
that Zfj,(^,t) = Zfj, (£, t) so the integrals over T + (t) and T_ (t) in ( |5.1 16[ ) are equal. 
We sketch a proof of ( |5.1 16| l. We have 



1 d 
2dt 



|9 c %(C,t)| dC 



Cer±(t) 



= 3? 



(C,t) {d t d A (Ztl (c,t) + iti(t)dl +1 z^t)}dc 



Cer + (t) 



Si 



Cer + (t) 



Cer + (t) 

= ft 

Cer+( f ) 



9 C 4 ^(C, t) {5 c 4 n w [J(-, t)}((, t) + iti(t)dl +1 z^ t)}dC 

9 c %(C,t) {n N [dp(;t)]({,t) + ih'{t)dl +1 z^,t)}dc 
d* Zll ((, t) [d^ j(c, t) + ^'(i)c» 4+1 ^(c, t)}d( 



(5.117) 
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since <9 4 z M (£, t), is a trigonometric polynomial in the range of 11^. 
Thus, Hn has disappeared from the right-hand side of ( |5.1 17| i. 

The above manipulations are made possible by the fact that the contours T± (t) are horizontal 



straight lines. That is why we took the trouble to make the conformal mapping in sections [57T 
and 



Using ideas from section 4. 1 we can compute d^J(C, t) modulo errors with controlled L 2 — norm 



on r + (t). To do so, we first recall our assumptions that the i? 4 — norms of Zi(C, t) — C and z 2 (£, t) 
on r + (t) are bounded by a controlled constant C, and that the chord-arc condition holds with a 
controlled constant c. Consequently, for £ e r+(t) and y e [— tt, 7r], the following estimates 
hold: 

|sin(«i(C,t)-2i(C + y,*))| <C|y|, |co8(zi(C,t)-«i(C + y,t)| <C (5.118) 
|sinh( 2l (C 2 ,t)-2 2 (C + 2/,t))| <C|»|, |cosh(z 2 (C,*)-z 2 (C + 2/,t))| <C (5.119) 
|5^(C,t)-^(C + y,i)| <C|y| forO<fc<2, M = l,2 (5.120) 
|^(C,t)-^(C + l/.*)l<C|y|M„(0 for /*= 1,2 (5.121) 

where M M (£) is the Hardy-Littlewood maximal function of 9 4 z M |r + (i). an d consequently 



|M M (C)| 2 dC < C. (5.122) 

r+(t) 



Also, our assumptions on the function 5(0 w, t) tell us that 

B(C,y,t) = B(C,C + y,*) fQrCer + (t),ye[-7r,4 (5.123) 

satisfies 

\d l c B(t,y,t)\<C\y\ forO<Z<4. (5.124) 

We now proceed to the computation of dfj(£, t), for ( e T + (t). 
By the proof of lemma [4~T] we have 



Jwer + (t) cosh(z 2 {(,t) - z 2 {w,t)) -cos(zi(C,i) - *i («>,*)) 

6 ^,„ax 

= Dangerous^, i) + ]TSafe^C,t) + £ c£Easy£ (C, t) , (5.125) 

3=1 v=\ 

with Dangerous^, Safe^ 1 and EasyJJ defined as in lemma|4.1l with k = 4. 
Using our estimates ( |5.118| l,...,( |5.122| ), one shows easily that 

|Easy^(C,*)| < C foreach//, z/. (5.126) 

To see how to estimate the "Safe" terms, we take for example 



Safest) 



(coeff) 



[9 c z 2 (C, t) - d c z 2 (C + y, t)} sinh(z 2 (C, *) - z 2 (C + y, t)) sin^ (0 t) - z 1 (( + y, t)) 



(cosh(z 2 ((,t) - z 2 (( + y,i),t) - cos(zi(C,t) - zi{( + y,t))Y 



x (a c %(C,t)-a c %(C + y,i))dy 



(5.127) 
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Since z\ (£, t) — C an d z 2 (d t) have H 4 — norms at most C, and since the chord-arc condition holds 
with a controlled constant c, one shows easily that the expression in curly brackets in ( |5.127| i has 
the form 6> ^' y ' t - ) , with estimates 

\6(t,y,t)\<C 

and 

\d y e(t, y ,t)\<c. 

Consequently, ( |5.127| i and the L 2 — norm boundedness of the Hilbert transform tell us that 



/ |Safe£(C,t)| 2 dC<C. 



T+(t) 

We can treat the other "Safe" terms in ( |5.125| l in the same way. 

Together with ( |5.126| l and the definition of the "Dangerous" term in ( |5.125| l, this now implies 

that 

tuer+W cosh{z 2 ((, t) - z 2 {w,t)) - cos( Zl ((, t) - Zl (w,t)) ^ 
+ Error M (C,t), (5.128) 

with 

/ |Erro r;l (C,t)| 2 rfC<C. (5.129) 

Thus, modulo errors with bounded I?— norms, we have computed the fourth £— derivative of the 
last integral in ( |5.114| ). 

We now turn our attention to the integral involving B((, w, t) in ( |5.1 14| i. Let < k < 4. As 
in the proof of lemma[4~T] we find that 



«fc f sin(^i(C,t) - zi(w,t)) , <m . 

Of / —, jr r 7 rr 7 r- r -. rr±S(C,W,t)dlV (5.130) 

Jwer + (t) cosh(z 2 (C,t) - z 2 {w,t)) - cos(zi(C,£) - «i(tw,t)) 
may be written as a linear combination (with harmless coefficients) of terms of the form 



[sii^zx (C, t) - z x (C + y, t))f [cos( Zl (C, t) ~ Zl (C + y, t))\ m ' 
[sinh(^(C, t) - zi (C + y, t))f [cosh(^ (C, t) - zi(C + y, t))] m " 

v' V 

d l c B(C,y,t)f[ [df Zl {^t)-df Zl {C + y,t)\ fj [dfz 2 ((,t)-dfz 2 (( + y,t) 

v=\ i/=l 

[cosh(z 2 (C, t) - z 2 (C + y, *)) " cos(zi(C, t) - zi(C + 1/, i)F m rfy (5.131) 
where m > 1, each k' v , k" > 1, and also 

k' + k" + 1 + v' max + v'; nax - 2m > 0, 

and 
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These terms ( |5.131| l can be estimated using ( |5.118| l,...,( |5.124| i, as we have just done for the 
terms Easy£(C, t) and Safe^ in fl5.125[ ), 

If each k' v and fc" is at most 3, then we can prove that the integral ( |5.131| l has absolute value 
at most C (this is the analogue of the "Easy" terms). If one of the k' u or k" is equal to 4, then we 
can prove that ( |5.131[ ) has L 2 — norm at most C as a function of £ € T + (this is the analogue of 
the "Safe" terms). Note that the analogue of the Danguerous term from lemma 4.1 is an Easy 
term in the present context, because four derivatives fall on the known, harmless function B. 

Consequently, for ( |5.130| >, we learn that 



< C 




ofc f sin(zi(C,f) -zi(w,t)) 

Of I —. 7- — r ~. rr -. 7- — : -. —B{(,W.t)dw 

Jwer+(t) cosh(z 2 (C,i) - z 2 iw,t)) - cos(zi(C,i) - Zi{w,t)) 
forC € T+(t), < k < 3; and 

o4 f sm(zi(C,t)-zi{w,t)) 

O c I — — r -, rr -. y- r 7 -r B{Q,W ,t)dw 

Jwer + (t) cosh(z 2 (C,i) - z 2 (w,t)) - cos(zi(Q,t) - zi{w,t)) 

has L 2 — norm at most C on r+(t). 
Therefore, 

smjz^Ct) - Zl (w,t))B(<;,w,t) \ 

j- r -. 77 -. 77 r -, rrdW0cZ u ((,,t) > 

Toer+(t) <nsh{z 2 {C,t)-Z2(w,t))-(X)s(z 1 ((;,t)-zx(w,t)) ^ J 

sin(zi(C, t) - z x {w, t))B((, w, t)dw g5z ^ +Error ^ t \ 

«er + (t) cosh(z 2 (C,t) - z 2 {w,t)) ~ cos(si(C,*) - zi(w,t)) c M M 

(5.132) 

with 

/ |Error^(C,t)| 2 dC < C. (5.133) 
Jr + (t) 

We have succeeded in computing the fourth £— derivative of another of the terms in ( |5.1 14| >. 

Finally, since the i? 4 -norms of zi((, t) — C ar, d Z 2(C t) on r+(t) are less than a controlled 
constant C, and since A((, t) has 4 derivatives bounded by C, we have 

dl{A((,t)d cZfl (ct)} = 

A(t,t)d 5 ( z^,t) +Error{ t (C,i) ) (5.134) 

for £ € r + (i), with 

/ |Error»(C,t)| 2 rfC < C (5.135) 
Jr+(t) 

Now, comparing ( |5TT28| ), gT29), ||5T32), < |5Tl33) , and (|5T34|, ( |5TT35) with the definition 
Q5 . 1 14] > of J, we learn that 

9 c 4 J(C,t)- 

\ Ler + (t) cosh(^ 2 (C,i)- z 2 ( w ,i))-cos(z 1 ((,t) - Zl (w, t)) J < > 

„. , r+(t) cosh(z 2 (C,<) - z 2 (w,t)) - cos(«x(C,t) - WMJj ^ 

Err M (C,i) forCer+(t), (5.136) 
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where 



/ |Err M (C,t)| 2 dC < C. (5.137) 
Jr+(t) 

34 



We have succeeded in computing d^J((,t) modulo an error with L — norm at most C on 

r+(f). 

We may rewrite dfj((, t) in the form 



[ Jwer + (t) cos^Z2[C,t) - z 2 {w,t)) -cos(zi((,t) - zi{w,t)) J 



2 , , n „ ^ ^2 — V 2 



u. G r + (t) [(d C 2i(C,t)) +(9 c z 2 (C,0) 

sin(zi(C,i) - Zi(iM)) 



cosh(z 2 (C,£) — Z2(w,*)) - cos(zi(C,i) — z\{w,i)) 



d^z fi (w,t)dw 



+ Err M (C,i) forCer+(i), (5.138) 

with Err^C, t) as in \5.U1) . 

Let us estimate the second-to-last integral in ( |5.138| >. 

The expression in curly brackets in that integral has C 1 — norm at most C, as a function of 
w 6 this follows by arguments in sections 4.3 and 4.4 



Therefore, integrating by parts in w, we obtain an integral of the form 

/ [Quantity of absolute value at most C] ■ d^z^(w, t)dw; 

Jwer + (t) 

such an integral has absolute value at most C. 

Consequently, ( |5.138| l may be rewritten in the form 

d c 4 j(C,t) = 

L(C,t) + P.V. [ S in(^(C, t) z l{ w t)) (g((>, t) + l)d W 1 

\ Ler+it) cosh(z 2 ((,t)-z 2 (w,t)-cos(z 1 ((,t)-z 1 (w,t))) f i^' J 



(d cZl (C,t)) 2 + (d c z 2 (C,t)) 2 ' 'J w er +W 1 ' ' ' i? 



+ Enf(C,i) forCer+(t), (5.139) 



with 



Now let 
Then 



|Enf| dC < C. (5.140) 



*) 

4. 



3.(0 = (5-141) 

/ |i^(C)| 2 rfC<C. (5.142) 

Jeer, ft) 
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We have 



P.V. I cot ( ^— - I d^{w,t)dw 



wer + (t) 



2 



-p.v. f cot (^) K(o - *»] ^ 

Jwer+(t) \ 1 / " 



/tuer + (t) 

= 2ttAF a1 (C) forCer + (t), (5.143) 

thanks to ( |2.1| i. Also, applying ( |2.16| l (with our present in place of h(x) in ( |2.16| >) we see 
that 

S c %(C,i) = i^(C,f) = -iAi^(C) + Go / M (C) (5.144) 

forC G r+(i), where 

/ \Goof^C)\ 7 dC<C-r^^2 [ mOfdC (5.145) 



We recall that > c; see our definition of controlled constants in section 5.3 Also, since 
i^(C)=^UO. we have 



/ |F M (c)| 2 dC= / \F^(C)\ 2 dC<c 

Jr,(t) Jr_(t) 



T+(t) j r _(t) 

by ( |5TT42> . 

Hence, ( |5.145| l yields 

/ \Goof^(0\ 2 dC<C. (5.146) 
Jr+(t) 

Now, putting ( |5.143| >, ( |5.144| > and ( |5.146| ) into ( |5.139| >, and recalling ( |5.140| >, we see that 

- i { A(( t) + P.V. [ sm(z 1 ((;,t)-z 1 (w,t))(B((,w,t) + l)dw \ Ap 

\ ' ' 'L e r + (t) cosh(z 2 (C,t) - z 2 (w,t) - cos( 2l (C,t) - z 1 {w,t))) j M 

27r9 f l(C ' f) ,AF,( C ) 

+ Erreur M (C, t) for C € r+(t), (5.147) 



with 



/ |Erreur M (C,t)| 2 dC<C. (5.148) 
Jr+(t) 

(Here, we have also used the fact that the quantity in curly brackets in ( |5.139| l is at most C, thanks 
to our assumptions on A(£, t), B((, t) and on z M (C, t).) 
Also, from ( |5.144| i and ( |5.146| ), we have 

ih'(t)dl Zfl ({,t) = ih'{t)Fl(Q = h'{t)KF^Q) + Goof e * tra {Q,t) (5.149) 

forC e T+(t), with 

f \ G oof^ tra ((:,t)\ 2 d<:<C. (5.150) 
Jr+(t) 
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From ( |5.147| ) ( |5.150| l, we see that 

sin(z 1 (C,t) - Zx(w,t)) (B(C,w,t) + 1) dw 



-iA(£,t) - iP.V. 

27Td C Z 1 (C,t) 



+ (t) cosh(z 2 ((,t) - z 2 (w,t) - cosOi(C,i) - zi(w,t))) 



m r,+tfM(» (^ + h ' {t) ^(O+Errr^C) (5.151) 
(d cZl (C,t)) +{d c z 2 (C,t)) J 

forC G r+(t), with 

/ |Errr^(C,*)| 2 dC < C. (5.152) 
Jr+(t) 

From $141) , $5AA2) and ( [5151) , < |57T52> , we find that 



/ #«„(c,t) [a c 4 j(c,t) + ^'(i)5 c 5 z M (c,i)] dc 

^Cer + (t) 



> 31 /" F^()R(()AF^Od( - C, (5.153) 



'Cer + (t) 
where 

sin(zi(£, i) — i)) (-B(C, w, t) + 1) du> 
'«>er+(t) cosh(z 2 ((,t) - z 2 (w,t) - cos(zi(C,i) - z 1 (w,t))) 



R(0=-iA(C,t)-iP.V. f 
2itd c z x {£,t) 



—*+ti(t). (5.154) 
(d <Zl ((,t)) 2 + (d c z 2 (Ct)) 2 

We note that R(Q has C 2 -norm at most C on T + (t). (Compare with < |4.58) , ( |4.62| > and ( |4.64| i. 
From ( |2.6[ ) and our present estimate ( |5.142| >, we see that 



J? 



/ F M (C)i?(C)AF A1 (C)dC > / F„(C)B(i2(C)) Af p (C)dC-C. (5.155) 
•/Cer+(t) ^Cer+(t) 



Moreover, 3t(R(()) > on thanks to the Rayleigh-Taylor hypothesis $537) and ( |5T98> . 

(Recall: We are solving backwards in time, and we stop as soon as RT fails.) Since also the 
C 2 — norm of 3fJ(i?(£)) on r+(i) is at most C, the inequality \2.6\ now tell us that 



F M (C)5R(J?(C))A^(C)dC>-C / |^(C)| 2 dC>-C, (5.156) 
7r+(t) 7r+(t) 

by ( |5~T42l >. 



From < [5.153) , < [5.155) , ( |5.156| ), we learn that 



3? / [d 4 J(C,t) + ih'(t)d 5 c z^,t)] d( > -C. (5.157) 

Jr + (t) 



Recalling ( |5.117| i, we conclude that 

/ \d 4 c z^(,t)\ 2 dC>-C. (5.158) 
Jcer.it) 



d 
dt 



Since we are solving backwards in time, this is our desired energy inequality. 
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We conclude that our Galerkin problem can be solved in a suitable time interval [to — S, to], 
maintaining control of the iJ 4 — norm, the chord-arc constant, and a positive lower bound for the 
Rayleigh-Taylor function. This control is uniform in the parameter N used to define our Galerkin 
approximation. Here, 8 may be taken to be a small enough controlled constant. 

It is now a routine exercise to deduce the following local existence theorem for the General- 
ized Muskat Problem, 

Lemma 5.4 Let h(t), I Ume , Q(t), fl, fi+ A(C,t), B((,w,t), z?(C), 4(C), and RT°(() be as 
in ( |5.76| l,...,( |5.98[ ). Then there exist functions zi(Q,t), z 2 (C,t) on fl, such that ( |5.99| l,...,( |5.108| l 
hold. 

This lemma provides only local existence, since hypothesis ( |5.98| l demands that the length of our 
time interval is less than a small enough controlled constant. 

5.5 Local Existence for Muskat Solutions in the Complex Domain 

In this section, we give the following local-in-time result. 

Theorem 5.5 1. Let h(x, t) be a positive real-analytic function on T x [t — 5, to + 8]; set 

r±(t)={(ec : 9fC = ±M«C.*)} 

and 

n(t)={(eC : |3CI <h(diC,t)}. 

Let 2°(C) = (^i(C)j z 2(0) be continuous on fl(to) closure and holomorphic on f2(to)- Sup- 
pose that Z®(C) — C and Z 2 (C) are 2ir— periodic and take real values for real C- Also, 
suppose that z®(() — C and £ 2 (C) belong to H i (T±(to)). 

2. Assume the chord-arc condition 

|cosh(z 2 °(C) - z° 2 (w)) - cos(z°(C) - 4H)\ >c[\mC + mC - w)\} 2 , 

for c,we n(t ) clos " re . 

3. Finally, assume the Rayleigh-Taylor condition i?T(C) > for all C € r_|_(t), where RT is 
defined by ( |5.26| l in terms of{zi{C > ), z 2 (C))- 

Then there exist positive numbers 8' , C , cf and, functions (zi(C, t), ^(C t)), continuous on 

{(C«)eCx [t -5', t ] : C G m c,omre } , 

for which the following hold. 

A. For each fixed t € [to — 8', to], the functions Zi(C, i) — C an d 22(C) are lit— periodic, 
holomorphic on Q(t), real-valued for real C. and have H A —norm at most C on T±(t). 
Moreover, 

|cosh(z 2 (C,t) - z 2 (w,t)) - cos(«i(C,t) - z 1 {w,t))\ > cf [||K(C- w)]\ + |3(C - w)\] 2 , 

for c,we n clos " re , te[t -8', t ]. 

Also, i?T(C, t) > (strict inequality) for all C G r+(i), i G [(o - 5', io]. 
(Here, RT((, t) is defined from zi(C, t), 22 (d i) fry ( j5.26| >). 
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B. The functions Zi(C, t), Z2(C, t) satisfy the Muskat equation 

d t z (C t) = [ sin(zi(C,Q - z 1 (w,t))[d c z fJ ,((,t) - d w z f _ l (w,t)} ^ 
* M ' Jwer+ cosh(z 2 ((,t) - z 2 (w,t)) - cos(zi((,t) - Zi{w,t)) 

for fj, = 1, 2, C G £l(t) domre , t £ [to — 5', t Q ], with initial conditions 

^(C,*o) = ^(0 

/or /u = 1,2 and ( £ Q(t ) closure '. 

C. The constants 8' , C , cf above may be taken to depend only on the following: 

• 8 (see 1.); 

• The function h(x,t); 

• Bounds for the —norms o/#i (C) — C and Z^iQ on r+(to)>" 

• The constant c in the cord-arc condition 2; and 

• A positive lower bound for the Rayleigh-Taylor function in 3. 



Proof: In section 5.2 we showed that the Muskat equation on a time-varying domain £!(£) is 



equivalent to a generalized Muskat equation on a time-varying strip. In section 5.3 and |5.4| we 



proved a local existence theorem for generalized Muskat. The present result follows at once. 

5.6 Continuing Back to Negative Time 

Let z = (z 1 (C, t), z 2 (C, t)), and A, r, A, k, h(x, t), H(x, t), Sl(t), T±(i) be as in section 3. Recall 
that H(x, t 2 ) < h(x, T 2 ) for all x. 

We suppose we are given initial data z°(() = (zi((), ^(0) on ^(T) closure , holomorphic in 
O(r), and real-valued for real (. 

We assume that the functions z"(£) — z^((,t) (jjl = 1, 2) are 2n— periodic, and that their 
i? 4 (i7(r))— norms are less than 10~ 2 A. 

Starting from the initial data z° at time t = r, we solve the Muskat problem backwards in 



time, until we reach a time £i east at which either the assumptions of theorem 4.3 break down, or 
else we reach ti east = r 2 . We can produce such a Muskat solution, thanks to a routine argument 



that combines the local existence result theorem 5.5 with the a-priori estimates given by theorem 
For t £ [ti east , t] C [t 2 , t], we have 



~\\z(-,t)-z(-,t)\\ 2 H < m) )> -n.DA-. 



Since also 

IW^-zO^ll^^lO^A 2 , 

it follows that 

- z{-, W)||^ ( n( tleM )) < HT 4 A 2 + C(A)rA 2 < 1(T 3 A 2 



Consequently, the assumptions of theorem 4.3 cannot break down at t = ti east , so we have 

^least = T 2 . 

Thus, we have solved the Muskat equation back to time t = t 2 . Our Muskat solution t) 
satisfies: 
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1 . It is holomorphic on 

n(f) = {|$*CI < fc(RC,*)} 

for each fixed t £ [r 2 , r] . 

2. We have 

IK,*) - z(;t)\\ B ^n(t)) < A 

for each t £ [t 2 , t] . 

3. In particular, for time t = r 2 , the Muskat solution z(C, t 2 ) is holomorphic on the smaller 
domain 

n(t) = {|9fd < W,*)} 

for t = t 2 (not the same as 1, due to a defect in our notation, noted in section^. 

4. Estimate 2 for the domain 1 implies easily that 

\\z(-,T 2 ) - z(-,T 2 )|| ff 4( n(T 2)) < ClA 

for the domain 3, as one sees from the Cauchy integral formula and the L 2 — boundedness 
of the Hilbert transform. Here, C\ is a universal constant. 

For the next few paragraphs, we define D,(t) by 3. 

Now, starting from the initial data z(-,t 2 ) at time t = t 2 , we solve the Muskat equation 
backwar ds in time, until we reach a time t\ east £ [— t 2 , t 2 } at which either the ass ump tions of 
break down, or else we reach ti east = — t 2 . (When we apply theorem 



4.3 



4.3 



we use 



theorem 
2C\ \ in place of A; see 4.) 

Again, a routine argument using theorem |5.5| and theorem |4.3| allows us to produce such a 
Muskat solution. 

For t £ [tieast, t 2 ], we learn from theorem 4.3 i 



that 



d 

— ||^-, L) - <-;iiH 4 (n(t)) 



z(; t)- z(;t)\\%Hn ( t)) >-C{A)r- l X 2 . 



Hence, by 4, we have 

lk(-,W)-2(-,ileas,)||^ (a(tlciisl)) < (%\ 2 + C{A)T-\Cl\) 2 T* <2C\\ 2 . 



(Notice that the large constant C(A)t^ 1 X 2 in the conclusion of theorem 4.3 is compensated by 
the short time interval of length 2r 2 ). 

As before, it now follows easily that the assumptions of theorem 4.3 (with 2C\X in place of 
A) cannot break down at ti east . Thus, we get all the way back to time ti east — — r 2 . 

We have therefore proven the following result. 

Lemma 5.6 Let z°(C) = (z® (C), ^§(0) be given, with z°(C) — t) Itt— periodic, holomor- 
phic in £1(t), and having norm at most X in 7J 4 (i7(r)). 

Then there exists a Muskat solution z((,t) = z 2 (C,£)) defined for t £ [— r 2 , t], 

holomorphic in Q,(t)for each t £ [— t 2 , t], equal to z°(£) at t = t, and satisfying 

\M;t) -^(•,*)||j?4 (n(t)) < CX 

for each t £ [— t 2 , t], where C is a universal constant. 
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5.7 Proof of the Main Theorem O 

In the previous sections, we worked with fixed parameters A, r, A, k > satisfying the conditions 
imposed in section[3] In this section, we will let k —> + . We write h n (x, t) to denote the function 
called h(x, t) before; and we write Vt K (t) to denote the region called Vt(t) in the previous sections. 
Let f K be the analytic function on Sl K (r) defined by setting 

«/.(«) - log (I) + (I)) - J. £ log (s to * (|) + sintf (S)) *, 
f K (x)dx =0, 

for x G T. One sees easily that can be taken to be analytic outside the union of the slits 

{5RC = 27rm, ±3C > «} (m G Z); 

thus, is analytic on K (r) as claimed. 

Moreover, the norm of / K (C) in ^ 4 (^ K (''")) is bounded by a universal constant C, since the 
function log \ x\ belongs to L 2 ([— it, it}). 

Lemma 5.6 from the preceding section gives us a Muskat solution 

2w(C>*) = (21,k(C,*)>32,k(C)*)) 

defined for < G [— r 2 , r] and £ G fi K (i), with initial condition 

^(C,r) = (z 1 (C,r) + A/ re (C),z 2 (C,r)). 
Moreover, z K (C, £) is holomorphic on f2 K (t) for fixed t, and we have 

||z K (-,t) - z(;t)\\ H 4(n K yfi < CX 

for t G [-T 2 , t]. 

A routine limiting argument now produces a Muskat solution 

*(C,t) = (*i(C,t),**(C,*)), ie[-r 2 ,r], (en K>0 fif ure (t) 

with the following properties 

1. Zi(C,t) =2 1 (C,r) + A/ (C), z 2 (C,t) = z 2 (C,t), where / (C) = lim K _ y0 + /«(£); 

2. z(£, i) is holomorphic in 

n(*) = n«> n«(t) 

for each fixed i G [— t 2 , t); 

3. | |z(-, i) - i)l| ff*(n(t)) < CX for t G [-r 2 , r]. 

Since 9 4 /o(£) nas a logarithmic singularity at the origin, it follows from 1 that the curve x n- 
z(x, t) is not C 4 — smooth. 

Recalling that d^z 2 (0, r) ^ 0, we see that there is no way to reparameterize the curve x H> 
z(x,t) to make it C 4 smooth. However, for t G [— r 2 , r), we see from 2 that the curve x H> 
2(2;, t) is real-analytic. 

Finally, for t — — r 2 , we learn from 3 that 

\d xZl (x, -r 2 ) - d x z x (x, -r 2 )| < CX, (5.159) 

for all x G T. 

Since, from ( |3.7| i, d x z^(x, — r 2 ) > C20T 2 for all real x, we learn form ( |5.159| l that 

d xZl (x,-T 2 ) > 0, (5.160) 

for all ,x G T, since A is smaller than cr 2 . Thus, at time t = — r 2 , our Muskat curve is a graph. 
The proof of our main theorem ( |1.1[ ) is complete. 
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